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ABSTRACT 


Of all large space structural systems, space stations present a 
unique challenge and requirement to advanced control technology. Their 
operations require control system stability over an extremely broad range 
of parameter changes and high level of disturbances. During shuttle 
docking the system mass may suddenly increase by more than 100% and during 
station assembly the mass may vary even more drastically. These coupled 
with the inherent dynamic model uncertainties associated with large space 
structural systems require highly sophisticated control systems that can 
grow as the stations evolve and cope with the uncertainties and time- 
varying elements to maintain the stability and pointing of the space 
stations. 

This report first deals with the aspects of space station operational 
properties including configurations, dynamic models, shuttle docking 
contact dynamics, solar panel interaction and load reduction to yield a 
set of system models and conditions. A model reference adaptive control 
algorithm along with the inner-loop plant augmentation design for 
controlling the space stations under severe operational conditions of 
shuttle docking, excessive model parameter errors, and model truncation 
are then investigated. The instability problem caused by the zero- 
frequency rigid body modes and a proposed solution using plant augmentation 
are addressed. Two sets of sufficient conditions which guarantee the 
globally asymptotic stability for the space station systems are obtained. 

The performance of this adaptive control system on space stations 
is analyzed through extensive simulations. Asymptotic stability, high 


iii 


rate of convergence, and robustness of the system are observed under the 
above-mentioned severe conditions and constraints induced by control hard- 
ware saturation. It is also found that further actuation level reductions 
can be achieved by using model switching and disturbance modeling 
techniques. 
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CHAPTER I 


INTRODUCTION 

1.1 Adaptive Control and Large Space Structures 

Most of the well-developed control theory, either in the 
frequency domain or in the time domain, deals with systems whose 
mathematical representations are completely known. However, in many 
practical situations, the parameters of the systems are either poorly 
known or time-varying. In such situations the usual fixed-gain 
controller will be inadequate to achieve satisfactory performance in 
the entire range over which the characteristics of the system may 
vary. Hence, some type of monitoring of the system's behavior, 
followed by the adjustment of the control input, is needed and is 
referred to as adaptive control. In other words, while a 
conventional control system is oriented toward the elimination of the 
effects of state perturbations, the adaptive control system is 
oriented toward the elimination of the effects of structural 
perturbations upon the performance of the control system. 

Model Reference Adaptive Control (MRAC) and Self Tuning 
Regulators (STR) are two principal approaches to the adaptive control 
problems. In MRAC, the objective is to make the output of an unknown 
plant asymptotically approach that of a given reference model which 
specifies the desired performance of the plant. With STR, a controller 
for a plant with assumed parameters is first chosen and the control 
gains are then updated with the recursively estimated parameters of the 
unknown plant. 
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MRAC can be classified into the following two types based on the 
adaptation method used: 

(i) Indirect Control, in which the plant parameters are 
estimated and these in turn are used to adjust the 
controller parameters to meet the performance requirement 
dictated by the reference model. 

(ii) Direct Control, in which no effort is made to identify the 
plant parameters, the control parameters are adjusted 
directly so that, the error between the plant outputs and 
the model outputs approaches zero asymptotically. 

Since the adaptive control systems are highly nonlinear closed- 
loop feedback systems, there is a distinct possibility that such 
systems can become unstable. In fact, the proof of global stability 
of adaptive control schemes had been a long standing problem for two 
decades and was not resolved until around 1979-1980. One ideal 
application area for adaptive control is in large space structures 
(LSS). The purposes of applying adaptive control to LSS are to 
reduce the structural and parameter sensitivities of the controller. 
This is due to the fundamental difficulties pf obtaining an accurate 
model for a distributed parameter system. One often has to 
deal with linearized reduced-order models; hence, a great deal 
of uncertainties in the mathematical model describing the dynamics of 
LSS exist. In addition, time variation in the parameter values are 
quite common in LSS environments. Slow time-varying effects may be 
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caused by structural settling, thermal distortions, or reorientations of, 
system or subsystems. Spontaneous changes can also happen, 
especially for space stations. 

1.2 Objectives and Motivations 

After the space shuttle, the next major space endeavor will be a 
permanent manned space station. The launching of an initial space 
station is planned for the early 1990's. By virtue of its mission and 
function, the space station is a large space structure with a very unique 
operational environment. As such, it suffers the same drawbacks as other 
large space structures. These are related to its large size, flexibility, 
and the way it is built and deployed. The size and flexibility prevent it 
from comprehensive ground measurement and test, which implies that pre-flight 
knwoledge of the spacecraft dynamics will be far from precise. In-flight 
system identification will enhance out knowledge on flight 
dynamics but it cannot totally eliminate the model parameter 
uncertainites. Structural flexibility implies infinite dimensionality; 
hence, model truncation is inevitable. With current technology, 
only a relatively small number of states can be handled in 

control design and state estimation. Previous studies of control 

of large space antennas, for instance, have concluded that 
destabilization can occur when the parameters of a design model 

deviate from those of the actual plant by a significant amount 

[1,2]. In addition to model errors, dynamic disturbances of many 
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orders of magnitude greater than those of the conventional 
spacecraft will be routine for space stations. Shuttle docking can 
cause an Instantaneous change of mass of more than 100% accompanied 
by a high intensity shock load. Station assembly, payload 
articulation, crew motion, launching and retrieving of satellites, 
etc., will all contribute to disturbance and model parameter 
variations. The adaptive control problems for space stations are 

summarized in Fig. 1. Although fixed-gain robust control designs can 
desensitize the system performance, it is only effective to a limited 
extent, i.e., for cases where parameter values deviate slightly from 
their nominal values [3]. All these have motivated us to investigate 
the feasibility of applying adaptive control techniques to space 
stations to maintain the stability and pointing. 

1.3 Literature Review 

Indirect adaptive control was first proposed by Kalman [4] in 
1958. The control task may be divided into two parts: a 

prestructured controller and a recursive plant identification 
scheme. The parameter estimates are used by the controller to 
compute appropriate control actions. Feldbaum [5] called this dual 
control because of the two steps involved in finding the control. In 
Ref. 6 Iserman et al. compared six on-line identification and 
parameter estimation methods that can be used for indirect control. 
Ljung [7,8] has proposed a general technique for analyzing the 
convergence of discrete-time stochastic adaptive algorithms, yet the 
problem of boundedness has not been resolved. More recently, Landau 
[9] and Mart in-Sanchez [10] have proven the stability of certain 
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Figure 1 Adaptive Control Problems for Space Stations 



types of indirect adaptive controllers. These results had been 
extended by Johnstone [11] in 1979. Narendra and Valavani [12] 

derived, in 1979, some indirect adaptive control laws by employing a 
specific controller structure and the concept of positive realness 
and showed that these laws are identical to those obtained in the case 
of direct control. Goodwin [13] has utilized a projection algorithm to 
obtain a class of globally convergent adaptive algorithms in 1980 and 
established global convergence of a stochastic adaptive control 
algorithm for discrete-time linear systems [ 14] in 1981. 

Direct adaptive control was first designed using the performance 
index minimization method proposed by Whitaker [15] of the MIT 
Instrumentation Laboratory in 1961 and has since then been referred 
to as the MIT design rule. The performance index used is the 
integral square of the response error. An Improved design rule with 
respect to the speed of response had been proposed in 1963 by 
Donalson et al. [16], who used a more general performance index than 
that of Whitaker. Winsor [17] had also modified the MIT rule in 1968 
to reduce the response sensitivity to the loop gain, at the expense 
of additional Instrumentation. Although some progress was made then, 
none of the design rules mentioned so far are globally stable. From 
then on, stability has become a major concern for subsequent studies. 

The most common application of stability theory to direct MRAC 
had been Lyapunov’s second method. The adaptive rule is obtained by 
selecting the design equations to satisfy conditions derived from 
Lyapunov's second method. Butchart and Shackcloth [18] first 
suggested the use of a quadratic Lyapunov function, which was 
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employed later in 1966 by Parks [19] to redesign a system formerly 
designed by the MIT rule. The use of a different Lyapunov function 
by Phillipson [20] and Gilbart et al. [21] has resulted in the 
introduction of feedforward loops that would improve the damping of 
adaptive response^ Unfortunately, all these algorithms are difficult 
to realize in practice because of the requirement of measuring the 
entire state vector, which is often impossible. 

Landau [22] was the first one to apply Popov's hyperstability 
criterion [23] to multi-input multi-output continuous MRAC systems 
subject to perfect model following conditions. He also used the same 
technique to treat the discrete-time MRAC problems [24,25], as did 
Bethoux et al. [26]. Anderson has given a lucid proof of Popov's 
hyperstability criterion in Ref. 27. An important contribution was 
made by Monopoli [28] (1974), who proposed an ingenious control 
scheme for continuous single-input single-output systems involving an 
auxiliary signal fed into the reference model and a corresponding 
augmented error between the model and the plant outputs, so that the 
use of pure differentiators in the algorithm can be avoided. 
However, as pointed out in Ref. 29 (1978), the arguments given in Ref. 

28 concerning stability are incomplete. Following the augmented 
error signal concept, Narendra [30] (1978) and [31] (1980) and Morse [32] 
(1978) and [33] (1980) succeeded in designing globally stable, asymptotic 
output tracking algorithms for continuous single-input single-output 
systems. Both Narendra and Morse have assumed that the relative 
degrees of the plant transfer function are known. Besides, Morse's 
algorithm is much too complex for use in practical applications. The 
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application of augmented error technique to discrete-time single- 
input single-output systems was made by Ionescu [ 34] , Narendra [35] 
and Suzuki [36]. Johnstone et al. [37] have extended Suzuki' s 
technique to solve some simple non-minimum phase problems by 
optimizing an augmented optimization criterion. Goodwin [ 13] took a 
different approach. Instead of relying upon the use of augmented 
errors or auxiliary inputs, he used the projection theorem to 
establish the global convergence of a class of adaptive control 
algorithms for discrete-time deterministic linear multi-input multi- 
output systems. 

Narendra and Valavani have proved that the direct and indirect 
control would arrive at the same result [38]. Using a typical error 
model, they also found that when all the signals in the plant are 
uniformly bounded, hyperstability and asymptotic hyperstability are 
achieved under exactly the same conditions as stability and 
asymptotic stability in the sense of Lyapunov [39]. 

Most stability proofs for direct model reference adaptive 
control systems have been restricted to single-input single-output or 
at best, multi-input single-output systems [40] (1975). Results 
pertaining to direct MRAC for multi-input multi-output continuous 
systems which do not satisfy the perfect model following conditions 
are limited. Also, the assumption made by the above algorithms that 
the relative degrees (difference between the number of poles and 
zeros) of the plant are known is too restrictive from the engineering 
point of view. 
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One particular adaptive control algorithm applied to multi-input 
multi-output (MIMO) systems was proposed by Sobel et al. [41] (1979) and 
[42] (1982). Using the Command Generator Tracker (CGT) law developed 
by Broussard [43] and a direct Lyapunov stability approach, they 
designed a direct MRAC algorithm that, without the heed for parameter 
identification, forced the error between the outputs of plant and 
model (which need not be of the same order as the plant) to approach 
zero. Like all other MIMO adaptive control algorithms, this one also 
requires that the number of controls equals the number of outputs and 
the plant input-output transfer function matrix is strictly positive 
real for some feedback gain matrix. 

As far as applications of adaptive control to large space 
structures is concerned, direct control is superior to indirect 
control. Since large space structures are infinite dimensional, 
the identification of all or a large number of parameters of the 
plant is clearly impossible or unfeasible. As described in Ref. 44, 
the adaptive controller must then be based on a reduced-order model 
(ROM) whose order is substantially lower than that of the plant. 
However, when the adaptive controller operates in closed-loop with 
the plant, it interacts with the unmodeled residual subsystems, this 
may cause great difficulties or disastrous problems [45]. 

The literature which deals with the application of adaptive 
control to large space structures is very limited. This is due to 
the difficulties associated with LSS in both the model truncations 
and parameter uncertainties. 
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Rohrs et al. [46] (1982) developed a method of analyzing 
stability ~ and robustness properties of a wide class of adaptive 
control algorithms for systems that have unmodeled dynamics and 
output disturbances. According to their investigation, none of the 
algorithms they tested including those of the widely recognized 
researchers are stable under these conditions. 

Bar-Kana et al. [47] (1983) applied and extended the algorithm 
proposed by Sobel et al. [42] to the control of large space 
structural systems, in which they have treated the problems of 
unmodeled dynamics and other model uncertainties. This work also 
suffers from some drawbacks, e.g., it cannot handle rigid body modes 
which always play an important role in large space sructural systems; 
the design and analysis are specific to the control of a simply 
supported beam while LSS are usually much more complex; the 
sufficient conditions for stability are too restrictive, etc. All 
these necessitate further investigation for the application of 
adaptive control to large space structures in general, and to space 
stations in particular. 

1.4 Outline of This Report 

In Chapter II, two space station configurations and their mass 
properties are described. Finite-element, dynamic models for both the 
two-panel and four-panel station configurations are developed in Chapter 
III. In Chapter IV, the adaptive control problem is formulated. The 
control structure is addressed in Chapter V. In Chapter VI, a direct 
model reference adaptive control algorithm together with the plant 
augmentation design is presented and two sets of sufficient 
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conditions for asymptotic stability of the system are derived. In 
Chapter VII, extensive performance analyses through simulation are 
discussed, and conclusions are summarized in Chapter VIII. 
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CHAPTER II 


CONFIGURATIONS AND MASS PROPERTIES OF SPACE STATIONS 

During the last few years, many space station configurations have 
been proposed and studied. The configuration development is indeed an 
evolutionary process, since there are so many factors that need to be 
considered and assessed against various configuration concepts [48,49]. 

Two configurations have drawn particular attention — the NASA 
Space Station Task Force Initial Operation Center (IOC) Baseline Configu- 
ration and the Concept Development Group (CDG) Split-Module Planar Configu- 
ration [49]. They are also referred to as the two-panel baseline 
configuration (or 6 degree-of-f reedom model) and four-panel planar 
configuration (or 19 degree-of-freedom model) , respectively, hereafter. 

2.1 Two-Panel Baseline Configuration 

Figure 2 shows the Task Force IOC Baseline Configuration which 
was developed by the Space Station Task Force with the Jet Propulsion 
Laboratory. The system dynamics of this configuration is dominated 
by the two very large solar panels measuring 250 ft x 40 ft and 
weighing 4000 lbs each on the ground. It also has a 50 ft x 10 ft 
radiator panel, a central bus structure consisting of a resource 
module, habitat module, logistics module, laboratory module, berthing 
truss, payloads, etc. The entire station weighs 134,000 lbs. Note 
that the lighter weight of this configuration compared with that of 
the four-panel configuration is due to the fact that fewer modules 
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Figure 2 Two-Panel IOC Baseline Configuration 


were considered for this configuration rather than the structural 
differences between the two concepts. 

The moments of inertia for this configuration are I xx ■ 
8.75xl0 6 slug-ft 2 , Iyy - 1.58xl0 6 slug-ft 2 , and I z2 - 8.60xl0 6 
slug-ft 2 . Due to the asymmetric design, the products of inertia are 
quite high, I X y ** -9.57x10** slug-ft 2 , Iy 2 ■ -4.89x10** slug-ft 2 , and 
T XZ " 5 .18x10** slug-ft 2 . With the selection of the reference 
coordinates as shown in Fig. 2, the center of mass has a high bias of 
X - 27 ft, Y - -2.3 ft, and Z - 5 ft. 

Figure 3 shows a similar configuration with additional modules 
and payloads attached. The 6 degree-of-f reedom (DOF) dynamic model 
that is developed in Chapter III and used extensively for feasibility 
analysis is based on the mass properties of the configuration shown 
in Fig. 2. As far as dynamic properties are concerned, the 6-DOF 
model is closer to that illustrated in Fig. 4. 

2.2 Four-Panel Planar Configuration 148,49] 

Figure 5 shows a version of the CDG Split-Module Planar 
Configuration. The basic difference between this configuration and 
that of Fig. 2 is that this configuration consists of 4 smaller solar 
panels (100 ft x 50 ft) and 4 larger radiator panels (70 ft x 20 ft) 
and, in addition, this configuration is a dynamically balanced design 
with structural symmetry. Since the solar panels are of much smaller 
size and a more reasonable aspect ratio, the structural strength 
improves and the fundamental modal frequencies increase. 

The main structure of this configuration measures 280 ft in 
length and it supports two resource modules, several pressurized 
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Figure 3 IOC Baseline Configuration with Payloads 
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Figure 5 Four-Panel CDG Split-Module Planar Configuration 


modules, a 30-ft service truss, and payloads* The pressurized 
modules are sized 22 ft x 14 ft diameter determined by the space 
shuttle payload bay size. The solar panels are hinged to rotate 
about the axes parallel to roll (X) and pitch (Y) axes, respectively, 
for solar inertial pointing The radiators are also hinged for 
articulation, and the core or the bus of the station is pointed to 
the nadir direction. Again due to its large size and flexibility, 
the solar panels are the dominant factor for the flexible body 
dynamics. 

Referring to Fig. 5, Table 1 lists the dimensions and masses of 
the major components. Using the body coordinates with origin placed 
at the geometric center of the station, the center of mass is. 


X 


cm 


Vi 



- 1.235 ft 


Y 

cm 


Z 

cm 


IVi 


i m i z i 

pT 


o 


o 


The self moments of inertia of each component are computed using 
the mass and dimension data shown in Table 1. Each of the components 
falls in one of the three basic structures — a rectangular plate, a 
rectangular cube, or a right cylinder. For simplicity, all the 
component masses are assumed uniformly distributed throughout the 
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respective component structure. Table 2 shows the self moments of 
Inertia for each component and Table 3 shows the component moments of 
inertia with respect to the center of the reference frame which is 
also the center of mass if we ignore the small offset along the X-axis. 
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Component Dimension and Mass for Four-Panel Configuration 
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Component Self Moment of Inertia for Four-Panel Configuration 
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CHAPTER III 


DYNAMIC MODELS FOR SPACE STATIONS 

The development of space station control technology requires 
results of the following nature: 

(1) First - order assessment of problems that will have 
significant effects on the space station designs. In order 
to guide the concept formulation and design, these types of 
results must be generated at a high rate of return. 

(2) Analysis of problems that involve multiple interactive 
elements which require somewhat complex implementation. 
Because of the complexity of the problem, more time will be 
required for development and implementation. 

For the first purpose stated above, a 6-DOF finite-element 
dynamic model is developed for the two-panel space station. Because 
of its tractablllty and reasonable simulation turnaround time and 
cost , it Is used extensively to evaluate adaptive control problems 
and performance. For the second purpose, a 19-DOF finite-element 
dynamic model is developed for the four-panel space station. 
However, the analysis in this report is, in general, 

configuration Independent. 

3.1 Finite-Element Model for the Two-Panel Station Configuration 

3.1.1 Dynamic Variables, Coordinates, and Parameters 

The 6-DOF finite-element model is shown in Fig. 6, in which only 
plane motions are considered. The motions of interest are the 
rotations about the X-axis which are tightly coupled with the 

displacements along the Z-axls. Hence, this model has six degrees of 
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Figure 6 6-DOF Finite-Element Model for the Two-Panel Configuration 



freedom — three linear displacements Z^, Z2, Z3 and three bending 
angles 8j, @2, 83 of the central bus and the two outer tips of the 
solar panels. The solar panels are modeled as two uniform beams with 
length L, linear mass density p, and flexural rigidity El. The bus 
and modules are modeled as a rigid core body with mass M2 and moment 
of inertia I2 located at the center of the structure. 

3.1.2 The Stiffness Matrix 

A finite-element technique is used to derive the stiffness and 
mass matrices for the model [50]. To obtain the stiffness matrix by 
using this technique, one starts by dividing the structure into a 
finite number of elements, the properties of each element are then 
determined. The properties of the entire structure are obtained by 
superimposing those of the elements at the associated nodes. There 
are two nodal points associated with each element, and two degrees of 
freedom at each node if only transverse plane displacements are 
considered. The deflected shape of the beam element may be described 
by a set of cubic Hermitian polynomials, as follows , when unit 
displacement is applied at the left end. 


*1 <*) - 1 - 3 (t) 4 2 (t) 

♦ 3 (x) - x (» - if 


( 3 . 1 ) 

( 3 . 2 ) 
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The shape functions for unit displacement at the right end are, 


<■> ■ 3 (if - 2 (if 
*4 (l<) "Hi • *) 


(3.3) 


(3.4) 


The general shape expression is the superposition of these functions, 


v (x) 


*1 (x> v a + ’I' 2 v b + ^3 ^ ®a + ^4 (x) 0 b 


(3.5) 


The stiffness coefficient associated with the beam flexure is 




El (x) (x) *V (x) dx 



(3.6) 


For a uniform beam segment, the stiffness matrix becomes 


F 

a 


6 

-6 

3L 

3L 

F b 

2EI 

-6 

6 

-3L 

-3L 

T 

a 

L 3 

3L 

-3L 

CM 

L‘ 

T b 


3L 

-3L 

L 2 

2L : 


(3.7) 


where F a , F b , T a , and T^ are the nodal forces and torques at nodes a 
and b, respectively, and v fl , v^, 0 g , and 0^ are t * ie corresponding 
nodal displacements — translation and rotation. The stiffness of 
the complete structure is obtained by adding the element stiffness 
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appropriately. For Instance, the structural stiffness of node i at 
which elements 1, m, and n are attached is 


k = k (1) + ; (ra) 
li K ii K ii 


k (n) 

K ii 


(3.8) 


where the hat implies that all the variables are expressed in a 
common global coordinate system, or they have been transformed to a 
common system from their local systems. 

Now we shall use Eqs. (3.7) and (3.8) to derive the stiffness 
matrix for the 6-DOF model. Consider the uniform beam element in Fig. 7 
If the two nodal points taken are located at its ends, there will be 
a total of four DOF — one translation and one rotation at each node. 
Through the finite-element method the stiffness matrix is 


" F l“ 


' Z l“ 


~6 3L -6 3L~ 


“ z r 

T 1 

- k A 

6 ! 

„ 2EI 

3L 2L 2 -3L L 2 


6 ! 

F 2 


z 2 

L 3 

-6 -3L 6 -3L 


z 2 

_ T 2_ 


_°2_ 


_3L L 2 -3L 2L 2 


- 0 2_ 


(3.9) 


where Fj and ?2 are the forces, Tj and T 2 are the torques applied at 
the nodes. For a similar beam element (see Fig. 7), the stiffness 
matrix is 
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Zi Z2 Z2 Z3 

El, L *2^ *2^ El, L ^3^ 



Figure 7 Uniform Bean Elements 
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(3.10) 


By using Eq. (3.8) to combine Eqs. (3.9) and (3.10), i.e. , adding the 
element stiffness at the joining point, the stiffness matrix K for 
the combined uniform beam elements (Fig. 8) is obtained as 
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Figure 8 Combined Uniform Beam Elements 
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3.1.3 The Consistent-Mass Matrix 


The consistent-mass matrix is the mass matrix for the 
distributed mass of the flexible structure. For a beam element of 
length L, mass density p(x), the mass influence coefficient m^j can 
be expressed as 


m. . 

ij 



p(x) 4 n(x) <>j(x) dx 


m. . 

Ji 


(3.12) 


where t|>^(x) and ipj Cx) are the shape functions. The term "consistent" 
signifies that this mass matrix is obtained using the same shape 
functions \J> # (x) as those used for deriving the stiffness matrix. The 
cubic Hermltlan polynomials are generally used for straight beams. 
Therefore, for the straight uniformly distributed beam element shown in 
Fig. 7, the consistent-mass matrix is 



• • 
Z, 



PL 

420 


156 

22L 

54 

-13L 

22L 

4L 2 

13L 

-3L 2 

54 

-13L 

156 

-22L 

-13L 

-3L 2 

-22L 

4L 2 



(3.13) 


Similarly, the consistent-mass matrix Mg for the adjoining beam 
element is 
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( 3 . 14 ) 

Hence, the combined consistent-mass matrix Mq for the 6-DOF model can 
be obtained by using an approach similar to that for obtaining the 
stiffness matrix as, 



( 3 . 15 ) 
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3.1.4 The Lumped-Mass Matrix and System Mass Matrix 

The consistent-mass matrix accounts for the effect of 
distributed mass only, the effect of concentrated mass is accounted 
for by the lumped-mass matrix Mg, 

M d = diag (0, 0, M 2 , I 2 » 0, 0) (3.16) 

The total system mass matrix is, then 

M - M c + M d (3.17) 

3.1.5 Equations of Motion 

The equations of motion for the 6-DOF model can then be written 
as 

(3.18) 

(3.19) 

where 

z p " l z pl e pl z p2 e p2 z p3 0 p 3J T 
■ [Zj 8j Z 2 0 2 Z 3 0j]^ 

F - [Fj Tj F 2 T 2 F 3 T 3 ] T 

Up ■ M-dimensional plant control Input vector 
Yp - M-dimensional plant output vector 


M Z + K Z » F ■ B u_ 
P P P 

y p ” c(° z p + z p ) 
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B =» 6 x M control Influence matrix 
C «= M x 6 measurement distribution matrix 

a is the weighting factor of the position vs. rate measurement 
3.1.6 Modal Coordinates and Modal Properties 

To carry out the analysis of this dissertation, modal 
coordinates are used. The modal model is obtained by setting F = 0 
in Eq. (3.18) and solving the eigenvalue problem. Let <J> be the 
normalized eigenvector matrix which is selected such that 


* T M* = I (3.20) 

4> T K$ « A (3.21) 


where A is the diagonal eigenvalue matrix. Let n be the modal 
amplitude vector, and Z p ** 4>h. Substitute this into Eq, (3.18) and 
premultiply both sides by <fr T , one has the following dynamical 
equation in modal form, 

n + An = $^® u p (3.22) 
After adding damping terms, Eq. (3.22) becomes 

n + diag (2? 1 a) 1 ,...2C6<Jg)h + diag(<o^,... ,a>^)n»$ T Bu p (3.23) 

The corresponding damped dynamical equation in physical coordinates 
can be obtained through transformation. Let D be the damping factor 
matrix, one has 
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(3.24) 


D ® <}>“ T diag ( 2? i 40 ! » • • • » 2 ?6“6^ 1 

Hence, the system equations of motion in physical and modal 
coordinates are, respectively 

j MZp + DZp + KZ p - Bu p 

( y p - C (az p + Z p ) 

and 

I n + diag (2CJUJ,..., 2C 6 ui 6 )n + diag (w^,..., u|)n 
y p » C(a<()n + <frn) 

The modal properties of the 6-DOF model are obtained by using 
the parameters given in Fig. 6. The modal frequencies are, 

■ 0 
u>2 = 0 

a> 3 - 0.04 Hz (3.27a) 

u A - 0.0637 Hz 
ui 5 - 0.3885 Hz 
w 6 - 0.3947 Hz 


(3.25a) 

(3.25b) 

= ^BUp (3.26a) 
(3.26b) 


and the mode shape matrix (eigenvector matrix) $ is 
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.872E-1 

.194E-1 

-.127 

-.961F.-1 

.178 
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-.352F-3 
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.720E-3 
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-. 545E-2 
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-.693E-3 

.155E-1 

. 301E-2 

-.360E-7 

.151E-2 

-.174E-5 

-.352E-3 

-. 157E-4 

.727E-8 

-.450E-3 

.126E-6 

-.117E-3 

-.886E-1 

.116E-1 

-.127 

.961E-1 

.179 

-.117 

-.352E-3 

-. 157E-4 

-.720E-3 

.719E-3 

.546E-2 

-.548E-2 


(3.27b) 

Figure 9 shows the mode shapes corresponding to the above mode 
shape matrix. Of the six modes, there are two zero-frequency rigid 
body rotational and translational modes, two first bending modes — 
symmetric and antisymmetric, and two second bending modes. 

3.2 Finite-Element Model for the Four-Panel Station Configuration 
3.2.1 Dynamic Variables, Coordinates, and Parameters 

Referring to Fig. 10, the main or backbone structure is modeled 
as two flexible beams rigidly attached to the core body (the bus); 
and the solar panels are also modeled as flexible beams with two 
beams jointed together and attached to the ends of the main 
structure. The two payloads, assumed rigid here, are hinge connected 
to the core body to form a balanced structure. To keep the model to 
a tractable size, the beams are assumed torsionally stiff, and hence, 
only bending angles and the associated deflections are modeled here. 

This model consists of 19 dynamic variables, 7 translations and 
12 rotations, as indicated in Fig. 10. Since the beams are assumed 
torsionally stiff, $ 2 represents the roll angles for the entire 
length of the beams associated with the south (north) panels. For 
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Figure 9 Modal Properties for the 6-DOF Two-Panel Configuration Model 
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Figure 10 19-DOF Finite-Element Model for the Four-Panel Planar Configuration 


the same reason, 64 is the pitch angle for the entire main structure, 
therefore, the following constraints apply: 

e 2 « e 6 - e 4 (3.28) 

l t • t 

The payload motions are modeled by the hinge angles, Yg x , Yg y » Yg x , 
Yg y , where in Fig. 10 , Yg x , Yg y , Yg x , Yg y are the corresponding 
inertial attitude angles. The translations at the c.m. of the 

payloads can be computed using the hinge angles.. Therefore, the 
translation variables are Zj, Z 2 , Z3, Z4, Z5, Z 6 , Z 7 , and the 
rotation variables are 0 lt <\> 2 , 83, 64, <t> 4, 65, • 67. Yg x * Yg y , Y 9x , 

I 

and Yg y . 

The model parameters are basically the flexural rigidity (El), 
the length, and the mass density for the beam elements in the model; 
the mass, moment of inertia and length for the payloads and the core 
station. The values for the mass -related parameters are computed 
based on information discussed in Chapter II. 

3 . 2.2 The Stiffness Matrix 

First, the 19 -DOF model is divided into three parts, the solar 
panels — south, the solar panels — north, and the main structure. 
The payloads will not be considered until later. The finite-element 
technique used is the same as that for the 6-DOF model. 

Consider Fig. 11 , the two "south” panels are modeled as two 
uniform beams, each 115 ft in length. Referring to Eq. ( 3 . 11 ), 
one ha 8 
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(3.29) 


Similarly, referring to Fig. 12, the stiffness matrix for the 
north panels is obtained as 


< 



(3.30) 


The elastic model for the main structure is shown In Fig. 13. 
The main structure is also modeled as two flexible beams, but with 
higher flexural rigidity than those of the solar panels. These beams 
are uniform and have a length of 140 ft each. The stiffness matrix 
for this structure is shown in the following equation. 
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Figure 11 Elastic Model for Solar Panels — South 



Figure 12 Elastic Model for Solar Panels — North 
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(3.31) 


Applying the same principle of Eq. (3.8), the stiffness matrices 
of Eqs. (3.29), (3.30) and (3.31) can be combined to yield the system 
stiffness matrix. Reorder the force and displacement vectors and the 
corresponding rows and columns of Eq. (3.29) as follows. 



(3.32) 


The purpose of this reordering process is to move T 20 and 83 to ttie 
top for later use and F 2 and Z 2 to the bottom so that the "south" 
panels can be combined with the main structure. 
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Reorder the rows and columns of Eq. (3.31) so that F 2 and Z 2 are 
placed on the top and Fg and Zg are on the bottom. 



(3.33) 


Reorder Eq. (3.30) so that Fg and Zg sit on the top and Tgg and 
6g at the bottom, 



(3.34) 
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(3.35) 


combine Eqs. (3.32), (3.33) and (3.34) one has 
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Now we shall use the constraints of Eq. (3.28) and 


t 29 = t 60 " t 40 

to combine T 60 and T 2e and rename it as T^, and 9 6 and 8 2 and rename 
it as 0 ^ and rearrange the rows and columns so that they appear 
between F 4 and and between Z 4 and $ 4 , respectively. In addition, 
we also exchange the orders of and F 0 , and and Z 0 . Thus if 
we let F 8 , Z 8 , and Kg be the force vector, displacement vector, and 
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the stiffness matrix, respectively, Eq. (3.36) becomes 


F s ~ * s 


(3.38) 


where 


F g = (Fjr T 10 F 3 T 30 ? 2 T 2 * F 4 t 46 t 4<(. f 6 t 6 <{> f 5 T 50 F 7 t 70> 
Z g = (Zj 0 j Z3 83 Z 2 <^2 Z4 64 $4 Z 0 Z 5 0 3 Z 7 ©7)^ 
and K is shown in Eq. (3.41). 
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where a 


and 6 


2(EI) S 



8 


2(EI) e 



e 


A quick check of symmetry, K g is symmetric as it is supposed to be. 
3.2.3 The Consistent-Mass Matrix 

Consider Figs. 14, 15, and 16 which show the distributed mass 
models for the south and north solar panels and the main structure, 
respectively. By using the same finite-element technique, the 
consistent-mass matrices for the component structures are, 

a. Solar Panel — South 
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Figure 14 Distributed Mass Model for Solar Panels 



Figure 15 Distributed Mass Model for Solar Panels 
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Employing a similar approach for obtaining the stiffness matrix, 
the following equation can be obtained: 


F s - M sc Z s 


(3.A5) 
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where 
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where 


420 


and b 


420 


( 3 . 47 ) 
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3.2.4 The Luraped-Mass Matrix and System Mass Matrix 

Mgc accounts for the distributed mass for the flexible structure 
but not the lumped mass associated with the rigid bodies* Let Mgp be 
the lumped-mass matrix for the station excluding the payloads, then 

M sd = diag (0,0,0,n,0,0,M 4 ,I 4yy ,I 4xx ,0,0,0,0,0,0) (3.48) ‘ 

where M 4 , I 4xx , and l4yy are the mass and moments of inertia of the 
core station defined in Fig. 10. 

The total mass matrix, excluding payloads is, 

M S “ M S C + M SD (3.49) 

The corresponding dynamic equation due to mass and inertia is 


F s - M S *Z S (3.50) 

3.2.5 Payload Dynamics and Hinge Torque Model 

The dynamic model for the payloads (identified as bodies 8 and 9) 
and the hinge coordinates are shown in Fig. 17. 

To include the payload dynamics and the dynamic interactions 
between the payloads and the station; the following expressions are 
obtained using the Lagranglan approach: 
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7gx * 7 8X - « HINGE ANGLE FOR "PAYLOAD 8" ABOUT X-AXIS 

7 8Y * 7 8Y " °4 = H,NGE ANGLE FOR "PAYLOAD 8" ABOUT Y-AXIS 
7 9X * 7 9X * *4 " HINGE ANGLE FOR "PAYLOAD 9" ABOUT X-AXIS 
7 9Y * 7 9Y " " HINGE ANGLE FOR "PAYLOAD 9" ABOUT Y-AXIS 


Figure 17 Payload Dynamics and Hinge Model 
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*4 ” (M 8 +M 9 )Z 4 + (MgLg-M^L^)^ + M8^8b^8y~^9^9b^9y (3-51) 

T 4 q + (MgLg-M 9 L 9 )Z* 4 - (l8y8 + ^8^8 + ^9ys + ^9^9^^4 (3.52) 

" ( I 8ys + M 8 L 8a L 8b + M 8 L 8b )Y *8y _ (I 9ys +M 9 L 9a L 9b +M 9 L 9b^9y 

44) - (Ig xs + ^9x8^4 “ *8xs*^8x “ *9xs Y *9x (3.53) 

F.qs. (3.51), (3.52) and (3.53) are used to replace F 4 , T 4 g, and T 4 ( j, 
In Eq. (3.50). 

The torques applied at the payload hinges are, 

X 8x “ X 8xs^4 + X 8xs^8x (3.54) 

T 9x ° x 9 X s*^4 + X 9xs^9x (3.55) 

T 8y “ ~ M 8 L 8b* z *4 + ( x 8y8 +M 8 L 8a L 8b +M 8 L 8b^ e *4 + (I 8ys +M 8 L 8b^8y (3.56) 
T 9y “ MgLgt,’^ + (Igys+MqLqgL^+MgL^)^ + ( x 9y 8 +M 9 L 9b^9y (3.57) 

t ? « » 

where Yg x » Ygyt Yg x t and Yg y are the hinge angles for payloads 8 and 
9 about the X- and Y-axes. 

The above equations are derived using the Lagrangian approach 
with more general assumptions and then linearized for small angles. 
The detailed derivation is shown in Appendix A. 

3.2.6 Equations of Motion 

l* e t Fp “ (Tg x , T 9x , Tgy, T 9 y) and Zp • (Yg x » Y9x , Y8y» ^9y) be 
the payload forcing and displacement vectors, the corresponding 
vectors for the system can be partitioned as follows. 
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F = 


(3.58) 
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The system mass matrix becomes, 


M ® M c + M d 


(3.59) 
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(3-64) 


The system 

stiffness 

matrix 


r k ' 

0 . _ . 


8 1 

15x4 

K - 

1 . 

— 


°, ,, 1 
4x15 . 

°4x4 


where K g is defined in Eq. (3.41). 


The equation of motion is 
MZ + KZ - F 


(3.65) 


(3.66) 
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3.2.7 Modal Coordinates and Modal Properties 

Let n(t), A, and $ be the modal amplitude vector, eigenvalue 
matrix, and normalized eigenvector matrix, respectively. Let Z(t) = 
<f>n(t), substitute this into Eq. (3.66) and premultiply Eq. (3.66) by 
<J> T , then * t M* = I and 4> T K4> = A, one has the following dynamical 
equation in modal form, 

n + An = 4> t F (3.67) 

where A = diag (u> | , . . . ,a)j 9 ) . Adding damping terms, Eq. (3.67) 
becomes, 

h + diag (2CJUJ,... + diag(wj[, . . . ,u| 9 )n “ <t» T F (3.68) 

The corresponding damped dynamical equation in physical coordinates 
can be obtained through transformation. Let D be the damping factor 
matrix, one has 

D = <T t diag ( 25 jOJj , . . . .aCjqWjgH” 1 (3.69) 

and the equation of motion becomes, 

MZ + DZ + KZ - F (3.70) 

For the purpose of control, let B and C be the control influence 
matrix and measurement distribution matrix, respectively, the system 
equation in physical and modal coordinates are, respectively, 


57 


(3.71) 


/ MZ + DZ + KZ = Bu p 
V y p - C (az p + Z p ) (3.72) 


and 

( n + diag (2C 1 a) 1> ...,25i9Wj 9 )n + diag (tdj , . . . ,ti>j 9 )Ti * 4> T Bu p (3.73) 

( y p - c(a$n + $n) (3.74) 

To obtain the modal properties, i.e., to determine the 
eigenvalues and eigenvectors, for the open-loop system, one can 
either free the hinges for the payloads or clamp them. For the 

latter case, a 15-DOF system results with 12 flexible modes and 3 
rigid body modes. For the former case, however, a 19-DOF system 
results since the payloads are considered rigid bodies and the hinges 
are freed, it ends up with 4 additional rigid or zero frequency 
modes. Since this does not yield additional information, only the 
clamped-hinge case is considered in this dissertation. 

The modal frequencies and mode shapes for the four-panel 
configuration with clamped-hinge case are shown in Fig. 18* These 
modes are divided into three groups. The first bending group 
consists of 6 modes with frequencies ranging from 0.115 Hz to 0.302 
Hz. These modes are formed with the first symmetric or antisymmetric 
bending of the three major structures, i.e., the two solar panel 
pairs and the main structure. The second bending group is caused by 
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Figure 18 Modal Properties for the Four-Panel Configuration Model 






Second Bending Group 
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Figure 18 Continued 


the second symmetric or the antisymmetric bending of the three major 
structures. The frequencies for this group are much higher than 
those of the first group, ranging from 1.67 Hz to 2.34 Hz. The third 
group consists of three rigid body modes with zero frequency. 

The structural and mass parameters used for generating these 
modes are shown in Fig. 10. The flexural rigidity (EI) S *» 9.48xl0 6 
lb-ft 2 has been used for the solar panels and a value of an order of 
magnitude higher has been used for the main structure. 

3.3 Frequency Characterization of Space Station Dynamical Systems 

With the availability of these space station models, the 
frequency characteristics of the various dynamical systems in the 
space station environment are identified as shown in Fig. 19. 

For a nominal orbital altitude of 400 km, the orbital period is 
92.61 minutes or 1.8x10“** Hz rate. For an altitude close to 400 km, the 
orbital rate will be inside the shaded narrow region in Fig. 19. The 
solar panel libration frequency for quasi-solar-inertial pointing 
[51] will be twice the orbital rate, as shown in Fig. 19. A low 
bandwidth attitude control system for the space station will have a 
bandwidth in the range of 0.001 Hz to 0.005 Hz. The two-panel low 
DOF model and the four-panel finite-element model are shown in Fig. 

19 with their modeled frequencies identified by vertical lines. The 
dashed regions extending the modeled modes represent the modal 
spectra that are not Included in the models. The payload attitude 
control systems for a range of applications will have bandwidth in a 
range centered at 1 Hz. The core body including the pressurized 
modules should have structural frequencies above 9 Hz. The figure 
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Figure 19 Frequency Characteristics of Space Station Dynamical Systems 






indicates that the spectral separations of the orbital rate, the 
attitude controllers, and the low frequency modes of the station 
structure are reasonable. However, the same cannot be said about the 
structural modes and the payload controls. For instance, the payload 
controller bandwidth falls between the modes of the first and the second 
bending groups. This result strongly suggests that decoupling control 
for the payload is required. 
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CHAPTER IV 
PROBLEM FORMULATION 

The space station, or the controlled plant can be represented by 
the following state space model: 

( 4 . 1 ) 

( 4 . 2 ) 

where x p (t) is the Np-dimensional plant states, u p (t) is the M- 
dimensional plant control inputs, and y p (t) Is the M-dimensional 
plant outputs. Here it has been implicitly assumed that the number 
of inputs equals the number of outputs. Physically, control inputs 
are the forces and torques generated by actuators, such as thrusters, 
proof masses. Control Moment Gyros, and momentum wheels. The outputs 
are measurements such as linear and angular displacements and rates 
measured by accelerometers, rate and integration gyros, etc. A p , B p , 
C p are the state, control influence, and measurement distribution 
matrices of appropriate dimensions. (A p ,B p ) is controllable and (Ap, 
C p ) is observable. 

A reference model that serves as a performance measure is 
required. This model is chosen to be asymptotically 
stable and can be represented by the following equations: 

( *m (t > " ’SB**) + B m u m (t > (*‘ 3 > 

( " Cm 


x p (t) " Ap x p< c ) + B p u p<0 
y p (t) = C p x p (t) 
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where x^t) is the l^-dimensional model states, u^t) is the M- 
dimensional model commands, and y m (t) is the H-dimensional model 
outputs. A^j, B m , C m are matrices of appropriate dimensions. 

Since the space station is infinitely dimensional while the 
dimension of the reference model has to be reasonably small for 
practical implementation, the following condition will be necessary 
for any adaptive algorithms for LSS (large space structures) or space 
stations : 

N p » N m (4.5) 

Define the output error between the plant and the model as 

e y (t) - y m (t) - y p (t) (4.6) 

Since the reference model specifies the desired performance of the 
plant (space station), the objective is then, without assuming the 
complete knowledge of the plant, to design an adaptation mechanism to 
generate a suitable plant control input u p (t), so that the plant 
output tracks the model output asymptotically, i.e., 

lim ey(t) ■ 0 (4.7) 

t-*» 
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CHAPTER V 


CONTROL ARCHITECTURE 

5.1 Control Architecture for the Two-Panel Configuration 

Referring to Fig. 6, the most rigid location on the station is 
its core on which inertial sensors, accelerometers and actuators are 
located. Control Moment Gyros (CMG's) are assumed and they are 
effective only for antisymmetric modes; the symmetric modes are 
controlled by the thrusters or proof masses. To gain controllability 
and compensate for vibrations of the large flexible panel structure, 
reaction wheels at the tips of the panels are postulated. 
Accelerometers and a target set of vibration sensors for relative attitude 
and rate measurements are also placed at the panel tips. Although the 
panel tips are far from ideal for locating hardware components, the choice 
is nil. For translational control, force actuators are required at 
the bus. 
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In order to apply adaptive control techniques to this 6-DOF 
model, EqSi (3.26a) and (3.26b), i.e., the plant equations, are 
written in state space form. Let be the modal amplitude vector, 
and define the plant state vector Xp as, 


x 


P 


we have 



(5.3) 



the corresponding Ap, Bp, and Cp matrices are 


6x6 


6x6 


-to 


PU 


-25 .<*> , 
Pi Pi 


- to 


p6 


0 - 2C p6“p6 



(5.6) 


(5.7) 
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(5.8) 


C 

P 


{aCi 1 C<|> ] 
P 1 P 


where $ p is the mode shape (eigenvector) matrix for the plant, to p k 
and are the modal frequencies and damping ratios, respectively. 

Rased on the control input vector u p given in Eq. (5.1), the 
control influence matrix B is given by 


0 0 0 0 
10 0 0 
0 10 0 
0 0 10 
0 0 0 0 


0 0 0 1J 


(5.9) 


Since the sensors 


and actuators are colocated, we have, in Eqs. 


(5.7) and (5.8) 


C 


B 1 


(5.10) 


5.2 Control Architecture for the Four-Panel Configuration 

With a similar but more complicated control architectural design 
than that used for the two-panel configuration, the control input and 
output vectors for the four-panel configuration are defined as 
follows: 
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u pi 
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Y-torque at left south panel tip 

u p2 


t 30 


Y-torque at right south panel tip 

u p3 


f 2 


Z-force at the root of south panel 
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X-torque (bending) at the root of south panel 

U P 5 


F 4 


Z-force at central bus 
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X-torque (bending) at the root of north panel 
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Y-torque at left north panel tip 
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Y-torque at right north panel tip 
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( 5 . 11 ) 


( 5 . 12 ) 


11 sets of sensors and actuators are used at 7 locations. The 
problem becomes more involved now because of the more complicated 
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structure and the three-dimensional motion instead of plane motion 
only. 

Again, for the application of adaptive control, state-space 
representations of Eqs. (3.73) and (3.74) are obtained. Let fi 
be the modal amplitude vector, and define the plant state vector Xp 
as, 
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P 



we have 
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the corresponding A p , B p , and C p matrices are 
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■ 2c pi5°pi5 


(5.16) 
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B, 


T 

<t> B 

P 

[aC* 


<V 


(5.17) 

(5.18) 


where 4>p, and Cp^ are the eigenvector matrix, modal frequencies 

and damping ratios, respectively, and the control Influence matrix B 


Is 


B - 


00000000000 
10000000000 
00000000000 
01000000000 
00 100000000 
000 10 0 00000 
0000 1 00 0 0 00 
00000100000 
0000001 0 000 
00000001000 
00 000000 100 
00000000000 
00000000010 
00000000000 
00000000001 


(5.19) 


Again, the sensors and actuators are colocated, hence 


B 


(5.20) 


in Eqs. (5.17) and (5.18). 


72 



CHAPTER VI 

ADAPTIVE CONTROL ALGORITHM 

The equations of the plant (space station) and the reference 
model are described in Eqs. (4.1)-(4.4), and are repeated below for 
convenience : 


plant ( x p (t) = Ap x p (t) + B p u p (t) (6.1) 

{ y p (t) - C p x p (t) (6.2) 

modei j i m (t) = A,,, x m (t) + \ u n (t) (6.3) 

y»(t) - C m x ra (t) (6.4) 


The output error is 

e y( fc ) - y m (t) - y p (t) ( 6 -5) 

The objective is to find Up such that 

lim e^(t) = 0 (6.6) 

t+« 

with the requirement that dimension Np » N„j and Ap, Bp, Cp are not 
completely known. 

The adaptive control algorithm under consideration here is an 
extension of that developed in Ref. 42. Since the Command Generator 
Tracker theory [43] formed the basis of that algorithm, it is 
introduced first in the following. 
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6.1 The Command Generator Tracker Theory 

The Command Generator Tracker (CGT) concept developed by 
Broussard Is a type of model reference control. What makes it 
attractive is that it allows perfect model output tracking even 
though the plant and the model are not of the same order. The basic 
idea of CGT is to generate a plant input Up(t) as a linear 
combination of the model states x^t), model inputs u m (t), and the 
output tracking error e y (t), i.e., 

u p (t) - S 2 ix m (t) + S22“m(t) + Kg e y (t) ( 6 . 7 ) 


so that the output tracking error asymptotically approaches zero. 
The gain matrices S 2 i and S22 are referred to as the CGT gains, and 
K g is selected to stabilize (Ap-B p K e Cp). 

When y p (t) *» y m (t) (i»e., perfect tracking occurs), the 
corresponding plant state and control trajectories are called "ideal 
trajectories" and are denoted as x p (t) and u p (t), respectively. When 
^(t) is constant* x*(t) and u*(t) are assumed to be linearly related 
to the model command Ujjj and model state x^t) as shown in the 
following equations: 



( 6 . 8 ) 


where Sjj, Sj 2 » S21 and S22 are matrices with appropriate 
dimensions. By definition, the ideal trajectory x*(t) is such that 


"hfor the more general case, refer to Ref. 54 . 
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y*(t) = C p x*(t) = C m x m (t) = y m (t) 


(6.9) 


and 




A p x J( fc > 


+ B P u J(t) 


( 6 . 10 ) 


To find the CGT gains S 2 i and S 22 » we first substitute Eq. (6.8) into 
Eq. (6.9). Thus, 

C p x J (t) - V s ll V c > + s 12 ««,> - c m x m (6-U) 


Hence, 


( 6 . 12 ) 

(6.13) 

Differentiating Eq. (6.8), 


c p S 11 “ ^ 

C p s 12 - 0 


x J (t) " S ll x m (t) + S 12 " S 11 *■<*> " S ll ( Vo (t) . + Vm> 

(6.14) 


and substituting Eq. (6.8) into Eq. (6.10) 


x p (t) " V S H X m (t) + S 12 %> + V s 21 **<«> + S 22 V 

" (Ap S 11 + B p s 21> Xm ( fc > + (Ap s 12 + B p s 22> u m ( 6 * 15 ) 


and comparing Eqs. (6.14) and (6.15), we have 
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S 11 An “ A p S 11 + B p S 21 


( 6 . 16 ) 


S 11 B ra “ Ap s 12 + B p S 22 


( 6 . 17 ) 


Rewriting Eqs. (6.12), (6.13), (6.16) and (6.17) in matrix form, we 
obtain 


Ap B p S 11 S 12 s 1 

C p 0 s 12 s 22 _ 


and then define 


n il fi 12 

n 2 i n 2 2 


S 12 s llAn ^11 B m 


S B P 

C P ° 


The resulting GGT equations become 


( 6 . 18 ) 


( 6 . 19 ) 


5 11 " n ll S llAra + n 12 C m 

5 12 “ ft ll s llAn 

5 21 “ n 21 S llAm + fi 22 C m 

5 22 " n 21 s U B m 


(6.20a) 


(6.20b) 


(6.20c) 


(6.20d) 


The so-called perfect model following (PMF) conditions are a 
special case of the CGT when the state vector is available and it is 
assumed that x£(t) “ ^(t). Since x£(t) - Sj jx^t) + S 12 u m , the PMF 
conditions imply that ■ I and Sj 2 * 0. Also, since 


*£ Ct) " V*) 


( 6 . 21 ) 
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• ) 


Apxj(t) + B p u p (t) 

■ Ap + ®p (S 21 x m (t) + S 22 u m ) 

■ < A p + B p s 21> V 1 ^ + B p S 22 “m 

" + B m “m 


( 6 . 22 ) 


From Eq. (6.22), it is found that 

(6.23a) 
(6.23b) 

If the CGT gains S 2 j and S 22 which satisfy the PMF conditions 
(Eq. (6.23)) exist, a valid PMF controller then becomes 

uj(t) + K(x m (t) - Xp(t)) (6.24) 

where K is a stabilizing feedback gain. 

Back to the general CGT, we define the state error as e x (t) • 
x*(t) - Xp(t) and seek a controller which guarantees that e x (t) ♦ 0 
as t ♦ •. Because when x p (t) - x£(t), we have C p x p (t) - Cp x*(t), 
hence 

y p Ct) 4 C„ x p (t> - Cp (t) 4 C m x Tn ( t ) 4 y n (t) (6.25) 

and that is exactly the desired result we require. 


v° - 


B p S 21 “ Am " Ap 

B p S 22 " ®m 
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With the error defined as e x (t) = x p(t) ~ x p (t), *-^ e error 
dynamics becomes 


4 x (t) «* ApxJ(t) + B p u*(t) - Ap x p (t) - B p u p ( t ) 
■ Ap e x (t) + Bp [u p (t) - Up(t)] 


(6.26) 


Note that 

e y (t) - y m (t) - y p (t) - C p (xp(t) - x p (t)) - C p e x (t) (6.27) 

Equation (6.26) will be used in the subsequent stability analysis. 

The CGT concept stated above is, however, non-adaptive. The 
purpose of adaptive control is to eliminate the need for a priori 
knowledge of the plant that is required in CGT. 

6.2 Direct Model Reference Adaptive Control 

Based on the CGT concept, Sobel et al. [42] developed the 
following adaptive control law: 

u p (t) - Kg(t) ey(t) + K^t) x m (t) + K„(t) u^t) (6.28) 

or, let r(t) be a 2M + Njj vector defined as 

r T (t) - [e y T (t) x^t) u m T (t)] (6.29) 
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and K(t) be the M x (2M+N ffl ) adaptive gain matrix defined as 
K(t) - [ K e ( t ) K^t) K^t)] 

Then 

u p (t) ° K(t) r(t) 

where K(t) is a combination of two types of gains, i.e., 

K(t) - K p (t) + K x (t) (6 .32) 

and Kp(t) is the direct gain, Kj(t) is the integral gain defined as 
follows: 

K p (t) = e y (t) r T (t)T (6.33) 

K x (t) = e y (t) r T ( t )T (6.34) 

with 

Kj(0) « K I0 (6.35) 

where T and T are the (2M+N m ) x (2*1+^) gain weighting matrices. Kjq 
is the initial integral gain. Note that u p (t) is highly nonlinear 
and its values are, in part, proportional to the cube of the output 
error e y (t) e y T (t) e y (t). 

The stability conditions of the adaptive controller are 
established via the Lyapunov direct method by considering a positive 
definite function given by 

V(e x» K I> " e x T (t) P e x (t) + Tr[S(K I -K)T-kK I -K) T S T ] (6.36) 


(6.30) 

(6.31) 
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where 


P is a Np x Np positive definite symmetric matrix 
K is a M x (2M+N,,,) dummy gain matrix 
S is a M x M nonsingular matrix 

a* 

Note that dummy gain matrix K does not appear in the control 

A* A* AS A* 

algorithm. It can be partitioned as K = lK g , 1^, K y ] so that 

* * *<*> - *e C p e x (t) + ^ xjt) + ^ u m (6.37) 

where K g , K x and K u are, like K, dummy gains. 

Introducing the control algorithm into the error equation given 
by Eq. (6.26), using Eqs. (6.31) to (6.33) and recalling from Eq. (6.8) 

that uj(t) = S n x m (t) + S 22 u m , gives: 

“ Ap + B p f s 2 1 Vt) + s 22 “ra 

- K x (t)r(t) - C p e x (t) r T (t) T r(t)J (6.38) 

The time dependence of the variables is omitted in the sequel for 
brevity. Thus the adaptive system is described by 

e x 

and 

Kl 


- Ape x + BpfSzjX^ + S2 2 Uni - Kjr - Cpe x r T Tr] (6.39) 

- C p e x r T T (6.40) 
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For the positive definite function given in Eq. (6.36), its time 
derivative is given by 

V = e x T Pe x + e x T Pe x + 2 Tr[S(Kj - K)T -1 K i t S t ] (6.41) 

Substituting Eqs. (6.39) and (6.40) into Eq. (6.41), we have 

V = e x T P[A p e x + B p S 21 x m + BPS 22 U™ - B p K x r - B p C p e x r T Tr] 

+ [Apex + B p s 21 x m + ®p s 22 u m “ B p K I r " B p c p e x rT ^ ]Tpe x 
+ 2 e x T C p T S T S(Kj - K)r 

- e x T (PAp + Ap T P)e x - e x T (PB p C p + C p T B p T P)e x r T Tr 

- 2 e x TpB P K I r + 2e x TpB p<S 2 l^ + S 22 %> 

+ 2 e x T C p T S T S(K X - K)r 

- «, T < PA p + V P)e x - *x T< PB P C P + C p T B p Tp) e x r T Tr 

ip ip ip 

- 2 e x [C p S T s - PB p ]K x r - 2 e x T C p T S T S K r 

+ 2 p B p<S 2 l*m + s 22%> < 6 * 42 > 

Choosing C p such that C p T S ■ P B p , Eq. (6.42) becomes 

V - ^(PAp + Ap T P)e x - 2 e x T PB p (S T S)‘ 1 B p T P e x r T T r 

- 2 e x T C p T S T S K r + 2 e x T PB p (S 21 x m + S^) (6.43) 
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Substituting -Eq. (6.37) into Eq. (6.43), we obtain 


V = e x T [P(Ap - R p K e Cp) + (A p - B p K e C p ) T P]e x 

- 2 e x T PB p (S T S)~ 1 B p T P e x r T T r 

+ 2 e x T PB p [(S 21 -K x )x m + (S 22 - K u )u m J (6.44) 

with the choice ^ = S 2 j and K u = S 2 2* Eq. (6.44) reduces to 

V - e x T [P(A p - BpKgCp) + (Ap - B p K e C p ) T P]e x 

- 2 e x T PBp(S T S)“ 1 B p T P e x r T T r 
Let Q - - P(Ap - BpK e C p ) - (Ap - B p K e C p ) T P 

Eq. (6.45) becomes 

V - - e, 1 Q e x - 2 e x T P B p (s T S)" 1 B p T P e x r T T r (6.47) 

If T is chosen to be positive semideflnite, the second term in 

Eq. (6.47) will be negative semideflnite in e x# Then we choose P 

~ 

such that for some K e , Q is positive definite. Consequently, V is 
negative definite in e x , and V is a Lyapunov function for establish- 
ing the asymptotic stability of the zero state of Eqs. (6.39) and 
(6.40). 


(6.45) 

(6.46) 
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C ~ 



To summarize, a sufficient condition for asymptotic stability is 

(1) T is positive definite 

(2) T is positive semidefinite 

(3) P B p = C p T (S T S) 

(4) P is chosen such that P(A p -B p K g C p ) + (Ap-B p K g C p ) T P 
is negative definite for some K g . 

Under the above conditions, the plant output will asymptotically 
track the model output. Furthermore, since the derivative of V is 
negative semidefinite in the augmented state (e x (t), Kj(t)], the 

adaptive gains will be bounded. 

It should be noted that conditions (3) and (4) together are 
equivalent to requiring that the plant transfer matrix 

Z(s) = C p (sl-Ap + B p K e C p ) B p be strictly positive real for some 

feedback gain matrix Kg. For the definition of positive realness and 
strictly positive realness of matrices, see Appendix B. 


6.3 Instability Problem Caused by Rigid Body Modes 

It is well-known that the zero frequency rigid body modes are 
unstable modes. Simulation results show that this adaptive control 
algorithm fails to stabilize these modes and yield stable states or 
outputs. These problems may be verified analytically. They are 
treated as critical cases of stability problems for autonomous 
differential equations developed in Ref. 52. Let 

x = Ax + g(x) (6.48) 
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be the given differential equation, where Ax represents the linear 
part and g(x) the nonlinear part. The real parts of the eigenvalues 
of A are assumed to be nonpositive. The critical case refers to the 
situation where some of the eigenvalues have zero real parts. If the 
matrix A has a double zero eigenvalue and if the reduced equation can 
be written in the form 

a, l =u 2 + (6.49) 
U) 2 =:aa) l +a,aJ l + ^ + ' * • _hl) 2^ tJl> l + ^ ,UJ l + ^ + * * (6.50) 


(where g is at least of second order, Y and 6 are positive integers), 
then we have the following sufficient conditions for instability due 
to Lyapunov [53], 

The equilibrium is unstable, if one of the following conditions 
is satisfied: 

(1) Y even; 

(2) y odd, a>0; 

(3) y odd, a<0, 6 even, Y > 5 + 1, b>0; 

(4) y odd, a<0, 6 odd, Y > 26 + 2; 

(5) y odd, a<0, 6 odd, Y - 26 + 1, b 2 +4(5+l)Y > 0; 

(6) the equation for u >2 contains only the term h(&>i,u> 2 )* 
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To apply the above criteria to our case, the system equations 
are rewritten in the form of Eqs. (6.49) and (6.50) by defining 


J - [C,,, -C p ] 


(6.51) 


a I x m 


(6.52) 


Then 


since 


C m x m - C p X P " Jx 


(6.53) 


r =1 x 


(6.54) 


xm 

0 T 


(6.55) 


Hence 


r Tr “ ey ey + *m ^ xn + ^um “ra 


x TjT T Jx + x TO " T xm x,,, + C (constant) 


(6.56) 


Therefore the system equations (6.1) and (6.3) can be rewritten as. 
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r. ■ 

X 

m 


A 

ra 

0 

A 

1 

e 

3 b 

m 

i 

• 

X 

P 


0 

A 

P 

A 

X + 

B Jx (x T J T T Jx + x T T X + C) 
p ey m xm m 


+ 


BLr 
P I 


(6.57) 


T _ 
e r T 

y 


(6.58) 


As stated in Chapter VII, in the present study, we concentrate 
on the space station attitude hold only, i.e., u,,, ** 0. In this 
special case, the corresponding oij's are x pi (modal amplitudes) and 
t^'s a*® Xpi (modal amplitude rates), respectively. For a certain 
mode, its corresponding block in matrix A is 

' 0 *1 
; 25 i“i "ij 

Hence for rigid body modes, it becomes [° °] which is exactly the 

form represented by Eqs. (6.49) and (6.50). Taking the 6-D0F model 
as an example, it is found through tedious expansion that y ■ 3, and 
a < 0; hence, the possibilities of satisfying condition 1 or 2 are 
excluded. One of the remaining conditions may still he satisfied 
although the expansion will be extremely tedious and is not done 
here. Note that the above conditions are sufficient conditions; 
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therefore, the equilibrium state of the system may still be unstable 
even though none of them are satisfied. 

If u m is equal to some constant instead of zero, the analysis of 
stability will be different. Since the term BpJxC in Eq. (6.57) is 
linear in Xp, it should be added to Ap. This changes the 
corresponding matrix A in Eq. (6.48) and its eigenvalues. The zero 
frequency modes would disappear, and the stability of the system will 
be entirely determined by the eigenvalues of matrix A. If any of 
them has a positive real part, instability of the system can be 
concluded immediately. 

6.4 Plant Augmentation 

To solve this unstable rigid body modes problem, a method 
referred to as plant augmentation is proposed. The plant 
augmentation is accomplished by introducing an inner control loop in 
the plant. 

Consider the equation of motion (Eq. (3.18) or (3.66)) before 
the damping terra is added, , 

MZp + KZ p - Bu p (6.59) 

After the inner loop is Introduced, the above equation becomes 

MZp + KZ p - Bu p - K il Zp (6.60) 

where K IL is the inner loop control gain matrix. By rewriting Eq. 
(6.60) as follows, 
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MZ p + (K + K IL )Z p - BUp 


(6.61) 


one can see that the modal characteristics of the plant have been 
altered due to By choosing the values and structure of Kj l , the 
rigid body modes will no longer have zero frequencies. As a result 
of this plant characteristic change, a stable adaptive control system 
can be realized. It is important to note that to design such an 
inner loop, one does not require accurate knowledge of the plant. 
This is because the inner loop controller can be made very robust by 
choosing the loop only at the location where the controllability . is 
the highest for the rigid body modes. Furthermore, the exact values 
of the augmented rigid mode frequencies are not important and what is 
important is that they are different from zero. Looking from another 
point of view, the stability of the adaptive system has been improved 
by the highly robust inner control loop. 

Taking the two-panel configuration as an example, the inner loop 
gain Kjl is chosen to have the following form: 

K il - diag(0, 0, K z2 , K 92 , 0, 0) (6.62) 

This selection is based on the fact that the rigid body modes are 
largely determined by the core body of the space station. To 
determine K z2 , we refer to the following block diagram: 
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Since the transfer function H can be written as 


H 



(6.63) 


the natural frequency for the rigid body translational mode can 

then be estimated as. 



(6.64) 


Similarly, the natural frequency 0 )q 2 for the rigid body rotational 
mode is 



Hence, the selection of ^2 and u 02 determines the values of K£ 2 and 

K 02 . 

Extending the above approach to the four-panel configuration is 
straightforward. is now of the form: 

■ diag(0, 0, 0, 0, 0, 0, Kg4» ^94* ^$4» 

( 6 . 66 ) 
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Including the effects of payloads, the natural frequencies of 
the rigid body modes are estimated as follows. 


0) 


Z4 


I ^ 

= J M~+M 


Z4 

„+M. 

8 9 


(6.67) 


(JL) 


04 


: / *94 


M 9 L 9 +I 8ys +I 9ys 


( 6 . 68 ) 


0) 


I 


4>4 


I 4xx +I 8xs +I 9xs 


(6.69) 


Again the selection of <dq^ and determines the values of k Z4» 

K 04 , and K^, respectively. 

With the inner loop introduced, the space station adaptive 
control system block diagram is shown in Fig. 20. 

6.5 Sufficient Conditions for Global Asymptotic Stability 

After solving the problem of zero frequency rigid body modes via 
plant augmentation, we need to find sufficient conditions that will 
make the space station adaptive control system globally 

asymptotically stable. Referring to Section 6.2, we have to select a 
positive definite P satisfying PBp ■ C p ^(S^S), and for which 

A# 

Q - - PCAp-Bp K e C p )-(Ap-Bp K e C p ) T P is positive definite for some K g . 
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REFERENCE MODEL 



Figure 20 Space Station Adaptive Control System Block Diagram 








Before deriving the sufficient conditions for globally 


asymptotic stability, we consider first the characteristic equation 
of a symmetric matrix for determining its eigenvalues. 

Let W be a symmetric matrix having the following partitioned 

form: 



(6.70) 


Let A be an eigenvalue of W. Then 


[i; 


(6.71) 


where [Vj V2]” is an eigenvector of W corresponding to A, i.e.. 


(6.72) 

(6.73) 

From Eq. (6.73) 

A 3 v j ■ ( A I - a 2 > v 2 (6.74) 

Hence 

V 1 " A 3~ ! (XI ” A 2 )v 2 <6 ' 75> 


A 1 W 1 + A 3 V 2 " X V 1 


A_v. + A-V„ - A v„ 
3 1 2 2 2 
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Substituting Eq. (6.75) into Eq. (6.72), we have 


A 


3 


i . e . , 


(XI - A J )v 1 

(xi - Aj)a 3 -1 (ai - a 2 )v 2 


[a 3 - (xi - Aj)a 3 -1 (xi - a 2 )]v 2 - 0 


Define the matrix inside the brackets in Eq. (6.77) as A,i. 

A = A 3 - (XI - A (XI - A 2 ) 

Since A is a diagonal matrix, in order to have a nontrivial 
v 2 , at least one element of A must be zero. 

Expanding Eq. (6.78), we have 


a - x 2 a" 1 - xa 3 _1 (a l +a 2 ) + (a 3 " 1 a 1 a 2 - a 3 ) 


To find the first set of sufficient conditions, let 


I ai 
ai I 


S T S - I 


K - K, K, 
e 11 


(6.76) 

(6.77) 
e. , 

(6.78) 
solution 

(6.79) 

(6.80) 

(6.81) 

(6.82) 
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then 



T T 
a<f> C 
P 


T T 
4> C 
P 


(6.83) 


and the condition PBp - Cp T (S^S) is satisfied. 

Since P is required to be positive definite, all of its 
eigenvalues must be positive. Referring to Eqs. (6.70) and (6.80), 

the corresponding A^, a **d A^ here are 


Aj = I 


a 2 “ I 


Ao ■ ol 


In view of Eq. (6.79), we have 


* “ I “4 (2D + (" I ‘ ®I) 


or 


A » X 2 I - 2X1 + (I - a 2 I) 


(6.84) 

(6.85) 

( 6 . 86 ) 


(6.87) 


( 6 . 88 ) 
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Hence a typical pair of eigenvalues are found by 


A 2 - 2A + (1 - a 2 ) « 0 


The roots of Eq. (6.89) are given by 


1,2 


2 ± A-4(l-a 2 ) 

■ s— “ 1 * ot 


Note that a > 0, hence 


Aj"l' + o>0 


For X2<-l-a>0,we require 


a < 1 


Referring to Eq. (5.6) or (5.16), we define 


w 2 0 



and assume the same damping for all modes, 
C - ?1 ** ?2 “ * ’ • “ ^Np 


(6.89) 


(6.90) 


(6.91) 


(6.92) 


(6.93) 


(6.94) 
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Then the matrix Ap can be written as 


0 I 


-A -2?A 


1/2 


( 6 . 


Substituting P, A p , B p , C p and K g into the equation for Q, 


Q - -P(Ap-B p K e C p )-(Ap-B p K e C p ) T P 



I 

al 


0 I 


0 

s - 

al 

I 

( 

-A -2?A 1/2 

• 

T T 
4> C 

Ti 


L -1 

r 


_ - 

1 

L ” J 


r ' 

\ -A -2?A 1/2 


T T 
<|> C 

LP J 


K i T|C i [ oC V c \fj 

t r 

) I al 
al I 


I 

al 

( 

0 I 

„ .1/2 

- 

0 

. T T T 

«■ 

0 

T T T _ , 

al 

I 


-A -2CA ' 


a$ C K. K,C<fr 
L p lip 

♦p c K i K i c *pJ 


( 


0 I 


0 

0 

T r 

V 

-A -2CA 1/2 


T T T 

_°*p C K L Kl C* p 

T T T 

♦p C K 1 K l% 

) 


I al 


al I 


I al 
al I 


-A-a<f> T C T K.^K.C$ -2CA 1/2 -«J» T C T K. T K.C* 
P lip P lip 


m T T T " 

0 -A-a<|> ‘C K, K,w 


I al" 

P IIP 



.1/2 T T T 



. 1 C K 1 *1% 


al I 
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2 T T T 1/2 T T T 

-aA-a <fr ‘C K. K.C* I- 20 CA 1 ' -a$ *C K. K.C* 

p lip p lip 

TTT 1/2 TTT 

-A-a<*> *C K, R.C* aI- 2 ?A ' -$ *C K. K.C* 

P 1 1 P P 1 1 P . 


-aA-a 2 <j> T C T K. T K.GJ> -A-o$ T C T K. T K.C4> 

P lip P lip 

1/2 TTT 1/2TTT 

I-2acA ' -a* C K, K.C* OI-2CA 1 ' -<J> C K, K.C$ 

pllp pllp. 


2 t T T 1/2 TTT 

ctA+a $ A C K, K,C<j> -I+2acA i/ %a<J> l C X K.C* 

pllp pllp 

T T T 1/2 TTT 

A+a<J> C K. K.C* -aI+2?A +4> CL K.C<f> 

Pllp pllp 


+ 


aA+a 2 <J> T C T K 1 T K C<{> 

P IIP 

1/2 TTT 
-I+2acA ' +a<{. C K, K,C<j> 

P IIP 


A+a<J> T C T K. T K C$ 

pllp 

1/2 TTT 
-aI+2?A ' +♦ C K, K,C<fr 
P IIP. 



2aA 

A-I+2o?A 1 ^ 2 


0 rp fp *P *P 

cT* T C\\C6 a4 l c\ l KC* 

p lip p lip 

a 



+ 2 



A-I+2acA^ 2 

-2aI+4cA 1/2 


•p fp ip *P *P *P 

«♦, c * P c K i K i c M 



2aA 

A-I+2acA 1/2 


’♦V 0 " 

T 

° K 1 

a 


+ 2 

P 


A-I^a^A 1 ^ 2 

-2aI+4£A 1/2 


T T 

0 4> C 

P J 

L k i J 



(6.96) 
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The second term in Eq. (6.96) is apparently positive semidefinite. 
Hence in order to make Q positive definite, the first term must be 
positive definite. Let 



2aA 

A-I+2o?A 1 ^ 2 

A-I+2acA 1 ^ 2 

-2aI+4?A 1/2 


(6.97) 


referring to Eq. (6.70), the corresponding A^, A 2 and A 3 here are 


Aj « 2aA 

(6.98) 

A 2 «* - 2aI+4cA 1//2 

(6.99) 

A 3 ■ A-I+2asA^2 

( 6 . 100 ) 


Substituting them into Eq. (6.79), we have 
A - A 2 (A-I+2a£A 1/2 ) -1 - A( A-I+2acA 1 / 2 ) _ 1 (2aA-2aI+4cA 1/2 ) 

+ (A-I+2otcA 1 / 2 ) - 1 (2aA) (-2aI+4cA^ 2 ) - (A-I+2acA 1/2 ) (6.101) 

In order that (A— I+ 2 aCA^ 2 )— 1 

|A-I+2acA 1/2 | * 0 (6.102) 

i.e., ta 2 -l+2aCa) * 0 , for i - 1, ... N (6.103) 

11 p 
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» 


for 1 = 1 


• • • • * 


(6.104) 


1 • € • , 


a * 



N 

P 


Expanding Eq. (6.101), we have 



r 1 o i 


1 

A 


2tmi 2 -2 o+4Cu>, 0 

-w 2 

A • X Z 

2 0 
-1+2 aCu^ 

0 * 

- X 

Wj^-l+ZaCw^ 

X 

Q 

u 

X 

X 

' l 


X 

A 

X 

X. 

X 

X 

^Np-^^Np 


%- 1+2ot x 



“Np- 1+2aC “Np. 


VI 


“2 " 

2aw (-2cH-4t« ) 


— 

2 

-l+2atw^ 

X 

X 

2 0 
-l+2atw^ 

* 


X 

„ ' 2o “np ( - 2, *«V 


X 

W Np” 2+2a ^ W Np 

“„p- 1+2t <“Np 



l.e., the characteristic equations are. 


1 


2 2 2 2 
2 2001 ^ -20+450)^ 2au> 1 (-2a+4c« )-(a>^ -l+2acu> 1 ) 


2 ” 2 
ai 1 -1+20^0)^ a) i -l+2a?u 1 


co^ -l+2a5u) t 


for i - 1, ..., N, 


(6 


or 


2 2 4 2222 22 3 

X -(200)^ -20+45^) X-(u i + l+4a z -2^ -405^+40 -4 acw^,) 


for 1 - 1, . .., Np 

Note that for the quadratic equation: 
aX 2 + bX + C = 0 


(6 


(6 


.105) 

-0 

.106) 

0 

.107) 

.108) 
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(6.109) 


the roots are, 


•v -bWb 2 -4ac 

X 1 2a 



Hence, In order to have Xj > 0 and *2 > 0 we must have 

b < 0 
ac > 0 


( 6 . 110 ) 


( 6 . 111 ) 

( 6 . 112 ) 


Referring to Eq. (6.107), the conditions for positive X^'s are then, 


- o(l-co|) > 0 , 


for i - 1,..., N p (6.113) 


1-4ocoj 1 -(2-4o 2 -4o 2 C 2 )u» 1 2 -4oc« 1 3 -hd 1 4 <0 , for i - 1, ..., N p (6.114) 

To summarize, the first set of sufficient conditions for 
globally asymptotic stability is given by 

T Is positive definite 
T is positive semidef inlte 


a < 1 

(6.115) 

2 

1-w 

ft + , 

(6.116) 

2^ 

24^ - o( l-w i ^) >0 

(6.117) 

l-4a5wj-(2-4a^-4a252j (1J ^2 - 4 a ^ (J j^3 +u ^4 < q 

(6.118) 


for 1 - 1, ...» Np. 
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To find the second set of sufficient conditions, we select 


P « 


“ 

81 al 
al I 


(6.119) 


S T S = I 


( 6 . 120 ) 


k i Tk i 


( 6 . 121 ) 


Then 


PB 

P 


81 al 


0 


r t t 

a<f> C 





P 




a 


al I 


T T 


, T T 


* c 


* c 



L P J 


P J 



( 6 . 122 ) 


and the condition PBp ** Cp^(S T S) is satisfied. 

Since P is required to be positive definite, we have, first of 

all, 

8 > 0 

and the eigenvalues of P must be positive, 
the corresponding Aj t A£ and A-j here are 


A 1 - 

81 

(6.124) 

A 2 ■> 

I 

(6.125) 

A ^ » 

al 

(6.126) 


(6.123) 

Referring to Eq. (6.70), 
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Substituting them into Eq. (6.79), it becomes 
2 

* - 1 - 1 < 6I+I > ♦ <1 1 - «> 

or 

A ® X 2 I - (8+1)1 + (8-a 2 )I 

Hence, a typical pair of eigenvalues are determined by 


(6.127) 


(6.128) 


X 2 - (8+l)X + (6-a 2 ) - 0 


Referring to Eqs. (6.111) and (6.112), 
Xj's are 


(6.129) 

the conditions for positive 


8 > -1 (automatically satisfied) 

and 

0 > a 2 


(6.130) 


(6.131) 


Substituting P, Ap, Bp, C p and K g into the equation for Q, we have 


0 ~ - P(A p - B p K e c p ) " (V B pVp )Tp 


81 ol 


al I 


T T T 1/2 T T T 

-A-a4 A C K. \ C* -2CA 4/ -♦ *C K. K.C4> 
pllp pllp. 
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_ TTT 

0 -A-a4> C K, K,C<t> 

pllp 


' BI al ’ 

.1/2 T T T 



I . -2CA ' -* p C Kj KjC^ 


1 

P 

M 

M 


r 7 T T T 1/2 T T T 

-aA-ct <j> l C R, K,C| Bl-2a?A w -a<J> 'c K, K,C| 

pllp pllp 

TTT 1/2TTT 

-A-a<j> L c K, K, C4> aI-2?A ' C K. K,C$ 

, P lip P llpj 

-aA-al^^ Kj Kj C<j> p -A-a^C Kj Kj C<J> p 

BI-2aCA 1/2 -a<j> T C T R T K « aI-2cA 1/2 -<j» T C T R T R C* 

pllp pllp 

aA+a 2 $ T C T K , T K C<|> -Bl+2aCA 1/2 +o* T C T K T K C4» 

pllp pllp 

A+a<J> T C T R. T K C<f> -aI+2?A 1/2 +(j) T C T K T K.C<J> 

pllp pllp_ 

aA+a 2 * T C T K T K Oj> 

P lip 

1/2 T T T 
-BrrtaCA 1 '*-*** l C R, R.Cf 

L p i i p 

2 TTT 

2aA+2a ♦ ‘c L LC* 

P lip 

1/2 TTT 

A-BI+20CA 1 ' +2a* C K. K.C* 

P lip 


A+a$ Vk/r.C* 

P IIP 

1/2 T T T 
-al+2th' +4 C K, R.C(j> 

P 1 1 PJ 

1/2 T T T 
A-8I+2acA 1/4 +2a$ C K, K.C* 
p ii; 

1/2 T T T 
-2aI+4tA ' +2$ C R. K.C* 

P IIP 
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Again the second tern In Eq. (6.132) Is apparently positive 
semidef lnite. Hence In order to make Q positive definite, the first 
term must be positive definite. Let 



referring to Eq. (6.70), the corresponding Aj, 
Aj - 2aA 

A 2 - - 2aI+4cA 1/Z 
A 3 - A-SH-2a£A 1/2 


A 2 and Aj here are 

(6.134) 

(6.135) 

(6.136) 
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Substituting them into Eq. (6.79), we have 


A - A 2 (A-0I+2acA 1/2 ) _1 -A(A-ei+2acA 1/2 )" 1 (2aA-2oI+4?A 1/2 ) 


+ (A-Sl+2acA 1/2 ) -1 (2aA)(-2aI+4?A 1/2 )-(A-ei+2o?A 1/2 ) (6.137) 


In order that (A-8I+2acA*/2)~l exists, we must have 


|A-6l+2ocA 1/2 | * o 


i.e., id -8 + 2a?(i) * 0 for i » 1, ..., N 

i i p 


(6.138) 


i.e. , a * 


8-u i 

2^T 


for i - I, . . . , N 


(6.139) 


Expanding Eq. (6.137), we have 



r 1 n 1 


r 1 n 1 


2a«,^-2a+4;w_ 0 

7 - 

2 0 
-0+2oCu^ 

V 

X 

X 

0 1 

- X 

2 ® 
Uj -8+2otu^ 

X 

X 

X 

0 1 


1 1 
N 

X 

X 

X 

0 2^-2«f4c U|Ip 


*V 6+2at “Hp_ 


*V 8+JatU N P _ 



2oBj ) 




'' 2«« (-2a+4c«) 

0 — 5e_1_ £E_ 


“hp • (42 % 


•, -e+2at», 


0 


(6.140) 
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Thus, the characteristic equations are, 


oj^ -8+2 


2 2 2 2 
* 2au. -2a+4?ai 2aa» (-2a+4cu> )-(w -B+2acw ) 

1 X+ i i — -0 

-8+2050^ -8+20^0)^ 


for i - 1, ..., N p (6.141) 


i.e. 


» 


2 2 4 2 2 2 2 

X -(200)^ -2o+4Cu i )X-(io i +8 +4a c 


Referring to Eqs. (6.111) and (6.112) 
Xj's are then, 

2 Coij - ad-oij 2 ) > 0, 

8 2 -4a8Co) i -(28-4a 2 -4o 2 C 2 )a) i 2 - 4a?a) i 3 


2 2 2 3 

-2Bu)^ -4a8Cu^+4a -4acu>^ )=>0 

for i - 1, ...» N p (6.142) 

, the conditions for positive 

for i - 1, ...» N p (6.143) 

+a) i 4 < 0, 

for i - 1, N p (6.144) 
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In summary, the second set of sufficient conditions for globally 


asymptotic stability is given by 


T is positive definite 

T is positive semidefinite 

6 > a 2 

8-u 2 
“ * 2 ^ 

2Co) i - aO-ujj 2 ) > 0 

8 2 - 4ag£;u>^ - ( 28-4a 2 -4a2£2) {1) ^2 _ 4a£ (1 j^3 +a j^4 < q 

for i “ 1 N p 


( 6 . 145 ) 

( 6 . 146 ) 

( 6 . 147 ) 

( 6 . 148 ) 


As we compare the above two sets of sufficient conditions, it is 
apparent that the second set gives us more freedom for selecting the 
system parameters. 
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CHAPTER VII 


PERFORMANCE ANALYSIS AND PRACTICAL CONSIDERATIONS 
The purpose of this chapter is to investigate the dynamic 
responses of the space stations under the severe disturbances as the 
adaptive controller described above is applied without assuming the 
complete knowledge of the system parameters. In addition, from the 
point of view of practical implementation, means of reducing the 
control efforts and their effects on the performance of the system 
are studied. 

Crew motion, reboost, and vehicle docking are the major 
disturbance sources. They will also cause changes of mass 

property. Crew motion will shift the center of mass, reboost 
will result in gradual mass reduction, and vehicle docking will 
spontaneously Increase mass and inertia of the system. From the 

point of view of time-varying effect and the level of disturbances, 
space shuttle docking is by far the most significant source of 

disturbance. Since most of the cases studied here are 
related to shuttle docking, the docking characteristics and 

devices are described first in the following sections. 

The simulation programs are written in Advanced Continuous 
Simulation Language (ACSL). The program and numerical results of the 
simulation of adaptive control during shuttle docking for the 19-DOF 
four-panel space station are listed in Appendices D and E, 
respectively. 
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7.1 Shuttle Reaction Control Subsystem Residual Rates 

The shuttle Reaction Control Subsystem (RCS) consists of two 
major parts, the primary (PRCS) and the vernier (VRCS) subsystems. 
There are a total of 44 thrusters, 38 of them are associated with the 
PRCS, each has a nominal thrust level of 870 lbs; and the other 6 are 
associated with the VRCS with a thrust level of 24 lbs each. Phase 
plane control laws are employed to determine when actuations are 
needed and jet select logics are used to determine what thrusters are 
to turn on. The states of the system are estimated by a two-stage 
state estimator with a dual cycle time of 80 ms and 160 ms. 

PRCS is normally employed for AV change, attitude maneuvers, and 
coarse attitude control; and the VRCS is for fine attitude control. 
Since shuttle docking requires maneuvers, PRCS must be used. Due to 
the high thrust level of the PRCS, and with more than one jet used at 
the same time to maintain attitude and approach rate while 
maneuvering, large residual rates result. The best achievable 
(minimum) residual rates, i.e., rates obtained under ideal 

conditions, are AV * 0.05 ft/sec and Am ® 0.20 deg/sec. However, 
these minimum rates are difficult to realize under nominal 

operational conditions and much higher rates are expected. These 
expected rates are on the order of AV « 0.50 ft/sec and Aw - 1.00 
deg/sec. Figure 21 shows the shuttle control system block diagram 
and residual rates. 

For the purpose of performance analysis, the following 

assumptions are made regarding the shuttle and space station docking: 
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• PRIMARY REACTION CONTROL SUBSYSTEM: FOR AV CHANGES, ATTITUDE 



o 



no 


Figure 21 Shuttle Reaction Control Subsystem (RCS) Residual Rates 


(1) Throughout the docking period, the space station attitude 
control system will maintain operational on attitude hold 
mode. 

(2) Just prior to the contact, the shuttle RCS is set at 
passive mode, i.e., no thrusters are allowed to fire. 

(3) Once contact is made, latching is assumed, i.e., no 
separation is allowed. 

7.2 Design of Shuttle Docking Devices 

Shuttle hard docking is a rather idealized condition. Under 
this condition, the shuttle momentum is transferred to the space 
station for a short period of time, At. At the end of At, the 
station and the shuttle are latched together as one integrated body. 

The initial momentum of the shuttle is determined by the shuttle 
mass, M s «* 7.81 x 10 3 slugs (2.52 x 10^ lbs) and inertia, I B » 7.54 x 
10^ slug-ft^, and the shuttle residual rate AV and Au> (see Fig. 21). 

The final velocities are, of course, determined by the combined system 
mass and inertia. 

The concept and design of a shuttle docking device which can 
simulate both hard and soft docking is shown in Fig. 22. The space 
station and shuttle are considered as two separate bodies coupled by 
a set of angular and rectilinear springs and dampers. Let M g and I g 
be the shuttle mass and inertia, and M 2 and I 2 the mass and inertia 
of the station. The values of the spring constants and damping 
factors can be computed using the following equations: 
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Figure 22 Concept and Design of Shuttle Docking Device 



(7.1) 



(7.2) 


(7.3) 


(7.4) 


where K , D , to., £ are the spring constant, damping factor, natural 

A A A A 

frequency and damping ratio for the set of rotational spring and damper, 

, D , a) , C T are lor the set of rectilinear 

Li Li Li 

spring and damper. The detailed derivation of Eqs. (7.1) - (7.4) is 
given in Appendix C. 

With this design, shuttle hard docking can be simulated by using 
extremely stiff springs and dampers while for soft docking, much 
weaker springs and dampers are employed* 

7.3 Performance of Adaptive Control on the Two-Panel Space Station 
As stated in Chapter III, the 6-DOF two-panel space station is 
used extensively to evaluate adaptive control problems and 

performance because of the associated reasonable simulation 
turnaround time and cost. These results provide valuable guidance to 
the understanding of the performance of space station adaptive 
control systems. 


respectively. Similarly, 
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7.3.1 Augmented Plant Modal Properties 

Following the approach described in Section 6.4, the following inner- 
loop control gains are chosen: 

= 45662 ft-lb/rad 
K z2 = 4464 lb/ft 

This plant augmentation has caused changes of modal properties from 
those shown in Eq. (3.27) of the unaugmented plant. The modal frequencies 
and mode shape matrix for this augmented plant are. 


“pl 

° 

0.01163 Hz 

“p2 

S3 

0.039 

Hz 

u p3 

- 

0.0656 

Hz 

“ P 4 

SB 

0.1684 

Hz 

“p5 


0.3892 

Hz 

“ P 6 


0.3947 

Hz 


and 


-.921E-1 

.128 

-.922E-1 

.144E-1 

.178 

.178 

.382E-3 

-.704E-3 

.705E-3 

-.264E-3 

-.544E-2 

-.549E-2 

.382E-9 

.177E-3 

.641E-8 

-.158E-1 

.185E-2 

•226E-5 

.332E-3 

-.745E-8 

-.465E-3 

-.959E-9 

•135E-6 

-.117E-3 

•921E-1 

.128 

.922E-1 

.144E-1 

.178 

-.177 

.382E-3 

.704E-3 

.705E-3 

.264E-3 

. 546E-2 

-.548E-2 



(7.6) 


The corresponding mode shapes are plotted In Fig. 23. As can be 
seen, the four lower frequency modes have changed substantially, and 
most Importantly, the two evolved rigid body modes no longer have 
zero frequencies. 
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Figure 23 Augmented Plant Modal Properties for the 6-OOF Two-Panel Configuration Model 





With the nodal dampings assumed to be ■ 0.5% for all modes, 
the Ap, Bp, and C p matrices in Eqs. (5.6), (5.7) and (5.8) can be 
readily determined. 

7.3.2 The Selection of the Reference Model 


To evaluate the performance of the adaptive controller, the 
reference model is selected to have lower order, significantly 
different model parameters, and high damping. It consists of 4 modes 


(corresponding to the 4 low frequency modes of the plant) or 8 states 


x,,,, 4 inputs u m , and 4 outputs y m which are defined as follows, 



u 


m 


U m3 

U m4 


m 


1 


' “Vi + ?«1 ' 

y m2 


aZ m2 + Z m2 

y m3 


a6 m2 + %2 

— i 


. ae m3 + 9 »3 . 


(7.7) 


(7.8) 


(7.9) 


the corresponding system matrices A^j, B,,, and are 


°4x4 

! *4x4 

. i _ 

| 

-<») 2 

nl X 2 
"“m4 

| " 2C ml“ml^ 

| ‘ 2? m4 a5 m4 . 


(7.10) 
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(7.11) 


(7.12) 


where B and C are shown in Eqs. (5.9) and (5.10), respectively. The 
modal frequencies for the reference model are, 


<o i ** 0.02 Hz 
mi 

o > m 2 “ 0.03 Hz 
“m3 ° 0.04 Hz 
“m4 ” °* 06 Hz 


(7.13) 


The modal dampings are = 0*707 for all modes. The position to rate 
measurement weighting factor a - 0.2 and the gain weighting matrices 
T - T - diag (2.5 x 10 9 , 2.5 x 10 9 , 2.5 x 10 11 , 2.5 x 10 9 , 1000, 
1000, 1000, 1000, 1000, 1000, 1000, 1000, 400, 400, 400, 400) are 
used. 

7.3.3 Adaptive Regulator Control 

7. 3. 3.1 Controller Performance with High Initial Transient 

The purpose here is to evaluate the convergence property of the 
adaptive controller for the attitude hold and vibration suppression 
under very large initial transient conditions. The initial conditions 
for the plant are: 
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pl 

- 3.699 ft 

Z P1 " - 

0.877 ft/sec 

pi 

0.860 deg 

®pl " 

0.336 deg/sec 

p2 “ 

0.345 ft 

Z p2 " 

0.035 ft/sec 

P 2 " 

0.937 deg 

0 p2 “ 

0.037 deg/sec 

P3 " 

4.071 ft 

Z p3 = 

1.045 ft/sec 

p3 “ 

0.723 deg 

0 p3 " 

0.387 deg/sec 


These conditions were taken 10 seconds after the shuttle docked to 
the space station with the shuttle initial approaching rates of 0.2 
deg/sec and 0.05 ft/sec. The corresponding values in the modal 
coordinates are obtained through the following transformation. 



(7.15) 


The initial conditions for the reference model are 


^mi - °* 9 V 
%i " °* 9 ^>pi 


(7.16) 


The simulation results for this esse are shown in Figs. 24-29. 
Figure 24 shows the outputs of the plant and the model. The model 
outputs damp out quickly because of its high damping, and the plant 


The remaining figures are at the end of this report. 
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outputs track them asymptotically and converge within 100 seconds 

from the transient starts. The same behavior is observed in Figs. 25 

and 26 for the plant and model physical states. In Fig. 25, Z 2 p 

follows Z 2m better than Zj p to Z lm or Zj p to Z 3 m because a force 

actuator is employed at the central bus; thus, at this position, 

linear displacement can be controlled more effectively. Figure 27 

indicates that the four lower frequency modes of the plant also 

closely follow the corresponding modes of the model although mode 2 

damps out at a slower pace. Mode 1, the rigid body rotational mode, 

is apparently the dominant mode of the system. The two high 

frequency modes shown in Fig. 28 also converge to zero within 100 

seconds. The required adaptive control inputs, u pi , are plotted in 

Fig. 29. The demand on bus control torque is quite high, almost 

3,000 ft-lb under these high initial transient conditions. 

7. 3. 3. 2 Controller Performance with High Initial Transient and 
Measurement Noise 

The measurement noises of the sensors are assumed to be zero- 
mean gausslan white processes. The level of the noises are taken to 
be two orders of magnitude lower than those of the peak outputs. Under 
these assumptions, the measurement uncertainties for the 
accelerometers are 0.012 inch and for the gyros are 20.63 sec* The 
qualities for the sensors in existence today are much better than 
these. For instance, the accuracy of the gyro used on High Energy 
Astronomy Observatory - 1 (HEA0-1) is 0.5 sec". 

The plant and model outputs are shown in Fig. 30. Apparently, 
the measurement noises have no significant effects on the high rate 
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of convergence and stability of the system although yp^ has some 
noticeable fluctuations throughout the time history. Figure 31 shows 
the adaptive control inputs. The bus control force Up 2 fluctuates 
highly during the first 25 seconds which is quite different from that 
of the case without measurement noises. After the initial transient, 
all the control inputs are associated with sustaining or even 
increasing fluctuations. Control energy is wasted by reacting to 
noises after the steady state is reached. One way to solve this 
problem is to employ a threshold at the actuator input. Control is 
applied to the system only when the input exceeds the threshold 
levels. Nevertheless, controller robustness is observed even with a 
high degree of measurement noises applied to the system. 

7. 3. A Adaptive Control During Shuttle Docking 

The following simulations are designed to test the controller 
performance and stability using shuttle docking dynamics. To 
simulate these cases, a docking disturbance term B 4 U 4 is added to the 
right-hand, side of Eq. (3.18), where 


r 0 01 



Lo 0 


and is the docking force, Tj is the docking torque applied to the 
space station. and T<j are determined by the following equations: 
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° d L ^shuttle “ z p2^ + ^shuttle “ Z p2^ 


(7.18) 


T d " d A ^shuttle “ ®p2> + K A < Shuttle ~ 8 P 2> (7.19) 

and the equations of motion of the shuttle are characterized by 
• • 

^shuttle * ” *d^s (7.20) 

^shuttle = " ^d^s (7.21) 

where M g is the mass and I 8 Is the moment of Inertia of the shuttle. 
The shuttle residual rates (rates prior to docking) are: 

^shuttle “ 0,05 ft /<*ec 

Shuttle (0) " °* 2 de g/sec 

7. 3. 4.1 Shuttle Hard Docking 

The control objective here is to stabilize the space station so that 
it returns to its prior docking condition after docking occurs. Thus 
all of the initial conditions for the plant and the model are set to 
zero. For regulator control, the input command u^ to the model is 
zero. Since there are no model disturbances employed under the 
"unscheduled event," the reference model stays at its quiescent state 
throughout the simulation period. Thus the plant, while under the 
influence of the docking disturbances, is commanded to follow the zero 
output. 


(7.22) 

(7.23) 
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The following docking parameters are used in this case: 


D l - 2.313 x 10 3 lb/ft/sec 
K l - 1.028 x 10 3 lb/ft 

, (7.24) 

D a = 1.86 x 10 ft-lb/rad/ sec 

K a - 8.25 x 10 7 ft-lb/rad 

the corresponding damping ratios and natural frequencies are 


linear: ■ 0.707 

o>L ® 0.1 Hz 


(7.25) 


angular: C A ■ 0.707 

<u A » 1.0 Hz 


(7.26) 


This case is termed "simulated hard docking** because the springs 
and dampers used are extremely stiff. The simulation results are 
plotted in Figs. 32-39. Figure 32 shows that all of the plant 
outputs converge to zero within 100 seconds after docking begins. 
The results here are surprisingly good since , in addition to the high 
docking disturbances and the parameter and truncation errors, we also 
deal with a sudden increase of plant mass and inertia by more than 
100%. The same performance is observed for the plant and model 
physical states in Figs. 33 and 34. In Fig. 35, all but mode 2 are 
well damped. Mode 2 shows a very lightly damped oscillation, 
however, since the amplitude of this mode is very small and it 
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presents no noticeable effects on any of the physical outputs. At 
higher amplitude, the damping is also higher, the oscillation will be 
decayed at a faster rate and it will pose no real problem to the 
system. The two high frequency modes in Fig. 36 also show good 
convergence rates. Figure 37 illustrates the time history of shuttle 
angular position and rate. It is not surprising to note that they 
appear almost exactly the same as the central bus angular position 
and rate (Figs. 34(b) and (d)) because after the hard docking, the 
shuttle and the station are supposed to combine and become one 
body. Figure 38 indicates that the peak of the docking torque is 
63,158 ft-lb. 

To achieve the above performance, a peak torque as high as 7500 
ft-lb has to be generated by the bus torque actuator (u p 3 > as shown 
in Fig. 39. This is extremely high and presents an implementation 
problem. Of course, hard docking will not be a realistic docking 
option due to its high transient loads to the space station. Under 
soft docking conditions, the torque should be much less. 

7. 3. 4. 2 Shuttle Soft Docking 

The conditions for the shuttle soft docking simulation are the 
same as those of shuttle hard docking except that much weaker springs 
and dampers are employed here. Specifically, the following docking 
Interface parameters are used: 
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I>L = 693.8 lb/ft/sec 


K l - 92. 49 lb/ft 

D A = 5.57 x 10 4 f t-lb/rad/sec 

K a - 742 ft-lb/rad 


(7.27) 


the corresponding damping ratios and natural frequencies are 


linear: Cl “ 0.707 

ii>L = 0.03 Hz 


(7.28) 


angular: C A = 0.707 

u A - 0.003 Hz 


(7.29) 


Because of the very low system natural frequencies yielded under 
soft docking, it will take a much longer time for the space station 
to reach its quiescent state. Hence, the simulation time of 400 
seconds is used in this case. The simulation results are shown in 
Figs. 40-49. Referring to Fig. 40, all plant outputs except y p 2 
converge slower than those of the hard docking case, especially yp 3 - In 
Figs. 41 and 42, Z p j, Z p 3 and 0 p2 also show a significantly slower 
convergence rate. The same behavior is observed for modes 1, 2, 3 
and 6 in Figs. 43 and 44. In contrast to the hard docking case as 
expected, the shuttle angular position and rate as shown in Fig. 45 
have experienced high excursions. However, with the adaptive control 
system on the space station, the angular position and rate of the central 
bus have deviated very little from their nominal values (Figs. 42 (b) and 
(d)). Figure 46 indicates that the peak of the docking torque is 190 ft-lb. 
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substantially lower than the 63,158 ft- lb of the hard docking — a 
factor of 332. 

What makes the soft docking a desirable' docking option Is, in 
addition to the reduced shock load, that the demand of control 
efforts drops drastically. For instance, the peak of Up^ drops to 300 
ft-lb from 7,500 ft-lb of hard docking (Fig. 47). The comparison of 
relative panel tip displacement and acceleration by employing 
adaptive control vs. low gain bus control are shown in Figs. 48 and 

49, respectively. The dynamic load and panel deflection are markedly 
reduced by using adaptive control. For example, with the adaptive 
control, the peak panel deflection drops to 0.41 ft. from 6.6 ft. 

7. 3. 4. 3 Shuttle Hard Docking with Actuator Saturation 

One way to relieve the high actuation demand on the control 
hardware is to employ actuator saturation through gain limiting. ■ To 
study the effect of actuator saturation, shuttle hard docking is 
again employed here (for soft docking will take too long to 
simulate). The following limits for the control inputs are used: 


1 “pi 1 

- 50 ft-lb 

1 °p2 | 

- 12 lbs 

max | » pJ j 

- 1000 ft-lb 

max | u . 1 

- 50 ft-lb 

1 P* 1 



(7.30) 
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Figures 50-54 indicate that the system performance has degraded 
compared with those without saturation (Figs. 32-36), as expected. 

The peaks of the responses are higher and it takes longer time 
periods for them to converge. It is also seen that the responses are 
more jittering due to the bang-bang effect. Larger angular excursion 
of the shuttle is also observed in Fig. 55. The docking force and 
torque remain unchanged (Fig. 56). Figure 57 shows that the control 
inputs Upi and Up^ have not yet totally stayed within their linear 
operation region at the end of 150 sec. 

As the gain limit becomes more severe (one half the values shown 
in Eq. (7.30)), the system performance degrades further as indicated in 
Fig. 58. Though it is believed that the outputs will still converge, 
it will take a much longer time than the above cases. Moreover, only 
Up 2 stayed within the linear operation during the simulation period 
(Fig. 59). 

Regardless of the severe actuation saturations applied, the 

results show no sign of threatening the system stability. 

7. 3. 4. 4 Shuttle Hard Docking with Model Switching and Disturbance 
Modeling 

Intuitively it is true that the control effort can be reduced if 
after the shuttle docking contact is made with the space station, one 
uses a new reference model in which the shuttle mass and Inertia are 
Incorporated and a simulated disturbance similar to the docking 
disturbance is injected. The plant outputs will now follow the new 
model transient Instead of the zero model outputs. The reason for 
the control inputs to be smaller is because of the expected smaller output 
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errors. This has motivated us to study the effect of model switching 
and disturbance modeling. 

Since the plant is assumed to be poorly known, after 
incorporating the shuttle mass properties, the mode shapes and modal 
frequencies of the 4 modes of the model are deliberately selected 
with an error up to 20%. The disturbance force and torque inputs to 
the reference model are taken to be 2.25 sec-pulse of 173.78 lbs and 
0.225 sec-pulse of 116953 ft- lb amplitude, respectively. 

The simulation results are shown in Figs. 60-65. The plant 
outputs follow the reference model outputs reasonably well. By 
comparing Figs. 60-64 with Figs. 32-36, it is found that the 
convergence rates are about the same for the two cases. As far as 
the peak control inputs are concerned, u p3 drops significantly from 
7500 ft-lb to 2880 ft-lb, and u p 2 drops from 120 lbs to 75 lbs as 
shown in Fig. 65. However, u pl and u p ^ rise from 230 ft-lb to about 
800 ft-lb. These results indicate that the concept of model 
switching is more involved than one may think. More study in the 
area of peak control effort reduction will be required. In general, 
one should not expect to reduce all the control efforts, but only the 
critical ones. 

7.4 Performance of Adaptive Control on the Four-Panel Space Station 

The purpose of this section is to demonstrate the performance of 
adaptive control on the much more complex 19-DOF four-panel space 
station model. 
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7.4.1 Augmented Plant Modal Properties 

Plant augmentation has altered the modal properties from those 
of the original plant. The new modal frequencies and mode shapes are 
shown in Fig. 66. The three evolved rigid body modes no longer have 
zero frequencies, and the evolved modes 1 and 2 become pure rotation 
about the X and Y axes, respectively. Note that the frequencies of the 
modes in the first bending group are also pushed higher while those 
of the modes in the. second bending group remain unchanged due to the 
fact that the inner-loop controller is a low-bandwidth system. 

With the modal dampings assumed to be Cp^ ■ 0.5% for all modes, 
the A p , Bp and C p matrices in Eqs. (5.16), (5.17) and (5.18) can be 
readily determined. 

7.4.2 The Selection of the Reference Model 

Again, to evaluate the performance of the adaptive controller on 
the four-panel space station, the reference model is selected to be a 
lower order system with high damping and significantly different 
parameters from those of the plant. It consists of 9 modes 
(corresponding to the rigid body modes and the modes in the first 
bending group of the plant) or 18 states 11 inputs and 11 

outputs y m which are defined as follows 



(7.31) 
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1 



y ra - 


“ml 

u m2 

“m3 

“m4 

“m5 

“ra6 

“m7 

“ra8 

“m9 

“mlO 

“rail 


y«nl 


a0 ml + 0 ml 

y m 2 


a0 ra3 + 0 ra3 

y m 3 


aZ m2 + Ki2 

ym4 


• 

°*m2 + *m2 

y m 5 


a ^m4 + ^m4 

y m 6 

- 

a0 m4 + 0 m4 

y m 7 


a< Nn4 + ^m4 

y m 8 


aZ ra6 + Z m 6 

y m 9 


• 

a ^m6 + ^m6 

ymio 


o0 m5 + 0 ra5 

_ y m ii_ 


_ a0 ra7 + 0 m7_ 


(7.32) 


(7.33) 
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The corresponding system matrices A m , B m and C m are, 


A 

m 


9x9 


^9x9 


-U) 


ml 




- 0 ) 


m9 


- 25 „ 9 > 


B 

m 




(7.34) 


(7.35) 



(7.36) 


where B and C are shown in Eqs. (5.19) and (5.20), respectively. To 

/ 

show the breadth of the controller performance, the modal frequencies 
are selected to have a 40% error and mode shapes a 30% error. The 
high modal damping ratios ■ 0.707 for all modes and position to 
rate measurement weighting factor a ■ 0.2 are employed. 

7.4.3 Adaptive Regulator Control with High Initial Transient 

The purpose here is to evaluate the convergence property of the 
adaptive controller for the attitude hold and vibration suppression 
under very large initial transient conditions. 
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The initial conditions for the plant are: 


Z P1 - 

- 1.649 ft 

Z P1 

“ 

0.283 ft/sec 

6 P1 - 

0.845 deg 

9 pl 

a 

- 0.179 deg/sec 

Z p3 3 

1.617 ft 

Z p3 

a 

- 0.284 ft/sec 

0 p3 

0.830 deg 

0 p3 

a 

- 0.180 deg/sec 

Z P 2 

- 0.004 ft 

Z P 2 

a 

0.002 ft/sec 

♦ P 2 “ 

- 0.004 deg 

♦ P 2 

a 

0.017 deg/sec 

v- 

- 0.001 ft 

z p4 

a 

0.0002 ft/sec 

V- 

0.722 deg 

0 p4 

- 

- 0.029 deg/sec 

♦p4 3 

0.012 deg 

• 

♦p4 

a 

- 0.035 deg/sec 

Z P6“ 

0.004 ft 

Z p6 

- 

- 0.001 ft/sec 

♦ P 6 3 

- 0.003 deg 

• 

♦p6 

a 

0.017 deg/sec 

Z P5" 

- 1.618 ft 

Z p5 

- 

0.283 ft/sec 

0 P5 3 

0.830 deg 

S3 

m 

- 0.180 deg/sec 

Z P7 3 

1.649 ft 

Z P7 

a 

- 0.284 ft/sec 

0 P7 3 

0.846 deg 

0 p7 

- 

- 0.179 deg/sec 


(7.37) 


These dynamical conditions were taken 10 seconds after the shuttle 
docked to the space station with the following Initial approaching 
rates: 

^shuttle “0.05 ft/sec 

^shuttle 3 0,2 de g/ 8 ec (7.38) 

♦shuttle 3 °* 2 deg/sec 
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The corresponding initial values in the modal coordinates are 
obtained through the following transformation, 



The initial conditions for the reference model are. 


(7.39) 


^rai - °*« V 
n ml = 0.8 n pl 


(7.40) 


The gain weighting matrices T - T * diag (2.5 x 10®, 2.5 x 10®, 

2.5 x 10 6 , 2.5 x 10 6 , 2.5 x 10 6 , 2.5 x 10 9 , 2.5 x 10®, 2.5 x 10 6 , 

2.5 x 10 6 , 2.5 x 10®, 2.5 x 10®, 100, 100, 100, 100, 100, 100, 100, 
100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 400, 400, 400, 
400, 400, 400, 400, 400, 400, 400, 400) are used. 

The adaptive system simulation for the four-panel configuration 
is implemented in a fashion similar to that for the two-panel 
configuration except that the dimension and complexity are much 
greater and the execution time and cost are high. Because of this, 
relatively fewer cases have been studied here. However, the results 
are very encouraging as shown in Figs. 67-73. The simulation period 
is 75 seconds. Figure 67 indicates that the plant outputs closely 
track the model outputs; and with the exception of y p 4 » y p 7 and y p 9 > 
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all outputs converge within 75 seconds. Figures 68 and 69 also show 

excellent physical state responses. The three bending angles $ ^ , $ 4 , 

• • • 

$6 and their rates $ 2 , $ 4 , converge slower than other states and 
rates. The modal responses are shown in Figs. 70-72. All except 
mode 5 show relatively high damping rates. Mode 2, the rigid Y- 
rotational mode, strongly dominates the system dynamics and has a 
very high rate of damping. This mode is excited due to the fact that 
the simulation assumes a huge shuttle angular momentum about the 
pitch axis. The demand for control efforts shown in Fig. 73 is 
substantially lower than that required for the two-panel 
configuration. The peak of control torque is 400 ft-lb at the 
central bus and 350 ft-lb at the panel tips as compared with 3,000 
ft-lb and 1,000 ft-lb, respectively, for the two-panel configuration. 
This is because the four-panel configuration is much more stiff and 
massive. 

7.4.4 Adaptive Control During Shuttle Hard Docking 

Unlike the hard docking case for the two-panel configuration, it 
is assumed here that the residual angular momenta of the shuttle are, 
about its roll and pitch axes. The docking disturbance term B^u^ is 
then defined as 
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0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


0 0 0 
1 0 0 


0 0 1 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 



( 7 . 41 ) 


where F d is the docking force, T d is the docking torque about the 
pitch axis (Y-axis), and P d is the docking torque about the roll axis 
(X-axis). Again the shuttle is assumed to dock with the space station 
at the central bus (the core). F d , and P d are determined by the 
following equations: 


F d - Dp (Z 8huttle - Zp 4 ) + K F (Z 8huttle - Zp 4 ) (7.42) 
T d - D t (Shuttle " ®p4> + k T ^shuttle “ 0 p4> (7.43) 
**d “ °P (^shuttle “ ^p4^ + ^p ^shuttle “ $p4^ (7.44) 
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and the equations of motion of the shuttle are characterized by 


^shuttle = “ *d^s 

(7.45) 

^shuttle = “ ^d/*9 

(7.46) 

• • 

^shuttle “ ” ^d^6 

(7.47) 

where M g = 7820 slugs is the mass; Iq = 7.54 x 10^ slug-ft 2 is the 

moment of inertia about the pitch axis; and 1^ “ 1.0 x 10^ 
is the moment of inertia about the roll axis of the shuttle. 
The shuttle residual rates (rates prior to docking) are: 

slug-ft 2 

^shuttle C°) “ 0.05 ft/sec 

(7.48) 

^shuttle <°> “ °- 2 deg/sec 

(7.49) 

• 

♦shuttle a 0,2 deg/sec 

(7.50) 


Again all of the initial conditions for the plant and the model 
are set to zero, and since there are no reference model disturbances, 
the reference model stays at its quiescent state throughout the' 
simulation period. Thus the plant is commanded to follow the zero 
output at the presence of docking disturbances. 

The docking parameters used and the corresponding damping ratios 
and natural frequencies are listed as follows: 


135 


linear: D p = 3.06 x 10 3 lb/ft/sec 

K p = 1.36 x 10 3 lb/ft 
C F - 0.707 
a F = 0.1 Hz 

angular: D<j> * 1.46 x 10 7 ft-lb/rad/sec 
(pitch) Kj. = 6.48 x 10 7 ft-lb/rad 

s 0.707 
oij = 1.0 Hz 


angular : . 

Dp “ 7.08 x 10 8 ft-lb/rad/sec 

(roll) 

K p = 3.15 x 10 7 ft-lb/rad 


Cp ■ 0.707 
o) p - 1.0 Hz 


(7.51) 


(7.52) 


(7.53) 


The gain weighting matrices are T - T ■ diag (2.5 x 10 6 , 2.5 x 10 8 , 

2.5 x 10 6 , 2.5 x 10 6 , 2.5 x 10 6 , 2.5 x 10 8 , 2.5 x 10 8 , 2.5 x 10 6 , 

2.5 x 10 6 , 2.5 x 10 6 , 2.5 x 10 6 , 1000, 1000, 1000, 1000, 1000, 1000, 

1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 1000 , 

1000, 400, 400, 400, 400, 400, 400, 400, 400, 400, 400, 400). 

The simulation results are shown in Figs. 74-82. 

Figure 74 shows that all of the plant outputs converge to 
equilibrium state within 300 seconds after docking begins. 

The same performance is observed for the plant states and their 
rates in Figs. 75 and 76. In Figs. 77-79, all but mode 5 are well damped. 
However, the amplitude of mode 5 is very small and it is believed 
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that it will converge after some longer period of time. Mode 2 
strongly dominates the system dynamics due to the huge shuttle 
angular momentum applied about the pitch axis. Due to the higher 
structural stiffness, the required control efforts (Fig. 80) drop 
drastically compared with those in the hard docking case for the two- 
panel configuration while the docking force and torque maintain 
unchanged (see Figs. 38 and 82). The linear and angular positions of 
the shuttle shown in Fig. 81 are practically the same as those of the 
station central hus (Figs. 75(d), (j) and (n)), as expected. 
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CHAPTER VIII 


CONCLUSIONS 

The feature of a large space deployable structure is its complex 
flexible dynamics. Flexible dynamics are characterized by extremely 
high system dimensions and parameter uncertainties. Model truncation 
plays an important role in spacecraft control due to the limited 
computer capability and related hardware available today or in the 
near future. Control systems that can adequately deal with truncated 
dynamics and model parameter errors are necessary for large space 
structural systems. Space stations, among all large space structural 
systems, have stringent operational requirements and present a unique 
challenge for control engineers and researchers. For space 
stations, in addition to the parameter uncertainties and 
truncation errors, the control system has also to deal with 
the growth, time-varying dynamics, and high-intensity environmental 
disturbances. Adaptive control provides a potential solution to 


these problems. 

A 

direct 

model 

reference adaptive 

control 

algorithm for 

the 

control 

of 

space stations is investigated. 

This algorithm 

along with 

the 

proposed inner-loop 

plant 


augmentation technique (for the unstable rigid body dynamics) 
form a potentially robust control system for space stations. Two 
sufficient conditions have been derived to assure the globally 
asymptotic stability. Extensive control simulations and analyses 
have been conducted for two space station configurations with 
emphasis on generic properties and practical implementation issues. 
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The adaptive system 

developed 

here 

has 

exhibited 

high 

performance and 

robustness. 

However, 

in 

common 

with all 

other 

adaptive control 

algorithms, 

this one is 

also 

nonlinear, complex 

, and 


high gain. The demand on control effort, especially during stringent 
high disturbance operations, is so high that it has exceeded the 
capability of the realistic hardware. To cope with this problem, a 
number of solutions have been proposed and investigated. These 
solutions are gain limiting, reference model switching, disturbance 
modeling, and disturbance load reduction (e.g., soft space shuttle 
docking). With these methods integrated into the controller, the 
required hardware capability has been drastically reduced yet 
stability and robust performance of the system are still observed. 
Further research in this area is required and fruitful results are 
expected. 

Specific conclusions are summarized as follows: 

(1) The study results show promising potential applica- 
tion of adaptive control techniques to space 
stations. 

(2) The proposed inner-loop plant augmentation method as 
part of the adaptive system has improved the system 
convergence significantly and stabilized the rigid 
body modes. 

(3) High rates of convergence and robustness have been 
observed throughout the simulated cases. 
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Specifically, 

(i) The system is robust a) in the presence of unmodeled 
dynamics — model truncation, b) with poorly known 
plant dynamics, and c) in the event of instant change 
of system mass property by more than 100%. 

(ii) It shows a good convergence property even under severe 
dynamic conditions including a) high initial elastic 
deformation and initial attitude errors, b) strong 
shuttle docking disturbances and dynamic interactions. 

(4) Shuttle hard docking is an example of the harsh space 
station operational environment. The high gain requirement 
associated with the adaptive controller will far exceed 
hardware limitations. Two methods of reducing the control 
demand have been investigated. Gain limitation can be 
applied to set a practical limit to the control effort and 
maintain stability of the system at the expense of higher 
transient and longer settling time. Model switching 
together with disturbance modeling provides a means for 
output error reduction and hence, the reduced control 
demand. However, since our knowledge about the plant is 
incomplete, tuning through simulation is required to 
achieve good results. 

(5) It is believed that a combination of load reduction (e.g., 
shuttle soft docking), gain limiting, model switching, and 
disturbance modeling will result in satisfactory system 
performance with realistic hardware implementation. 
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Figure 24 Adaptive Regulator Control for Two-Panel Configuration 





Adaptive Regulator Control for Two-Panel Configuration 




Figure 26 Adaptive Regulator Control for Two-Panel Configuration 



Adaptive Regulator Control for Two-Panel Configvirat 
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Figure 28 Adaptive Regulator Control for Two-Panel Configuration 



146 


Figure 29 Adaptive Regulator Control for Two-Panel Configuration 



Figure 30 Adaptive Regulator Control with Measurement Noise for Two-Panel Configuration 



PRECEDING PAGE BLANK NOT FILMED 


149 


Figure 32 Adaptive Control During Shuttle Hard Docking for Two-Panel Configuration 




Figure 33 Adaptive Control During Shuttle Hard Docking for Two-Panel Configuration 
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Figure 34 Adaptive Control During Shuttle Hard Docking for Two-Panel Configuration 
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Figure 36 Adaptive Control During Shuttle Hard Docking for Two-Panel Configuration 
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Figure 37 Adaptive Control During Shuttle Hard Docking for Two-Panel Configuration 
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Figure 39 Adaptive Control During Shuttle Hard Docking for Two-Panel Configuration 




Figure 40 Adaptive Control During Shuttle Soft Docking for Two-Panel Configuration 
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Figure 41 Adaptive Control During Shuttle Soft Docking for Two-Panel Configuration 



Adaptive Control During Shuttle Soft Docking for Two-Panel Configuration 
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Figure 43 Adaptive Control During Shuttle Soft Docking for Two-Panel Configuration 
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Figure 44 Adaptive Control During Shuttle Soft Docking for Two-Panel Configuration 
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Figure 45 Adaptive Control During Shuttle Soft Docking for Two-Panel Configuration 
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Figure 46 Adaptive Control During Shuttle Soft Docking for Two-Panel Configuration 



With Adaptive Control With Low Cain Bus Control 
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Relative Panel Tip Displacements 
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Relative Panel Tip Acceleration 




Adaptive Control During Shuttle Hard Docking with Actuator Saturation for Two-Panel 
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Figure 52 Adaptive Control During Shuttle Hard Docking with Actuator Saturation for Two-Panel 





Adaptive Control During Shuttle Hard Docking with Actuator Saturation for Two-Panel 
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Figure 54 Adaptive Control During Shuttle Hard Docking with Actuator Saturation for Two-Panel 





Configuration — Adaptive Control Inputs 




Adaptive Control During Shuttle Hard Docking with More Severe Actuator Saturation 







Adaptive Control During Shuttle Hard Docking with More Severe Arto.-t-nr Saturation 
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Figure 60 Adaptive Control During Shuttle Hard docking with Model Switching and Disturbance Modeling 



Adaptive Control During Shuttle Hard Docking with Model Switching and Disturbance Modeling 
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Figure 62 Adaptive Control During Shuttle Hard Docking with Model Switching and Disturbance Modeling 
for Two-Panel Configuration — Plant and Model Rotational State Responses 
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Figure 63 Adaptive Control During Shuttle Hard Docking with Model Switching and Disturbance Modeling 
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Figure 6 4 Adaptive Control During Shuttle Hard Docking with Model Switching and Disturbance Modeling 
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Figure 65 Adaptive Control During Shuttle Hard Docking with Model Switching and Disturbance Modeling 
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Evolved Second Bending Hroup 




r-j 
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Figure fifi Continued 
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Adaptive Regulator Control for Four-Panel Configuration 
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68 Adaptive Regulator Control for Four-Panel Configuration 
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Figure 6B Continued 










Lgure 69 Adaptive Regulator Control for Four-Panel Conf iguratlon 
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Figure f>9 Continued 










Figure 71 Adaptive Regulator Control for Four-Panel Configuration 




Figure 72 Adaptive Regulator Control for Four-Panel Configuration 
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Figure 72 Continued 
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Figure 73 Adaptive Regulator Control for Four-Panel Configuration 



Figure 73 Continued 





Figure 73 Continued 
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Figure 74 Adaptive Control During Shuttle Hard Docking for Four-Panel Configuration 
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Figure 74 Continued 
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Figure 75 Adaptive Control During Shuttle Hard Docking for Four-Panel Configuration 















Figure 7S Continued 
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Figure 76 Continued 




Figure 7 f> Continued 
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Figure 76 Continued 
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Figure 79 Adaptive Control During Shuttle Hard Docking for Four-Panel Configuration 




Figure 79 Continued 




Figure 80 Adaptive Control During Shuttle Hard Docking for Four-Panel Configuration 





Figure 80 Continued 




Figure 80 Continued 




Shuttle T, Inear and Angular Position 




Figure R2 Adaptive Control During Shuttle Hard Docking for Four-Panel Configuration 


APPENDIX A 


DEVELOPMENT OF HINGED PAYLOAD DYNAMICS USING 
LAGRANGIAN APPROACH FOR FOUR-PANEL CONFIGURATION 



Lagrange's equation: 



(A.l) 


Starting from taking position vectors for points A, B and C 


then 


—A " Z 4 —A 


-B “ Z 4 *A + L 8a -B 


" Z 4 *A + L 8a -B + L 8b 


* A " Z 4^A 


-B * Z 4 —A L 8a -B 
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Z 4 —A * L 8. (9 4 i x S»> 


Z, e + L d 6. el 
4 -A 8a 4 -B 


e 

— c 


-C " Z 4 -A + L 8a 6 4 -B + L 8b (Y 8Y - x -C ) 


Z 4 -A L 8a 6 4 -B + L 8b Y 8Y 


Similarly, 


— D 


Z. e. + L n 6. e' 
4 —A 9a 4 -D 


— E " — A * ^9a 9 4 — D + ^*9b * 9Y — E 


• v 2 ■ V • V 

•• v c -c -c 


K 9 *4 * L 8b ■'EY * 2 *4 L 8. *4 ( Sa • Sj> 


‘ 2 2 4 L Sb 'BY ( — A * 9 2 L S. L 8b 9 4 'BY <4 • si’ 


* 2 4 9 L 8 a 9 4 9 L 8b x 8Y- 2i 4 L 8. 9 4 £OB9 4 


" 2 z 4 L 8b t 8Y C0S V 8Y + 2 L 8a ^8b ®4 Y 8Y C08 Y 8Y 


Similarly, 


v: - 


V • V 
-E -E 


2 4 9 L 9a 9 4 9 L 9b '9Y 9 2 2 4 L 9« *4 «* 9 4 


+ 2 Z 4 L 9b Y 9Y COS Y 9Y + 2 L 9a L 9b 6 4 Y 9Y COS Y 9Y 
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Hence 


2 ®4 Z 4 + 2 m 8 V C + 2 “9 V E + I *4YY 6 4 + I I 8YS Y i 


1 * 2 ' 1 *2 1*2 1 *2 
2 X 9YS Y 9Y + 2 X 4XX *4 + 2 X 8XS Y 8X 4 2 *9XS Y 9X 


1 *2 1 *2 2*2 2 * • 9 

2 "4 2 4 * 1 "8 Z 4 * L L * I 4 «* * ’sv) 


‘ 2 Z 4 L 8a ®4 C0S ®4 “ 2 Z 4 L 8b (y 8Y 4 ®4 ) COS ( ®4 4 Y 8Y > 


+ . 2 L 8a L 8b ®4 (y 8Y + ®4 } C0S Y 8Y 


By 4 2 “9 ^4 4 L 9a ®i 


4 L 9b (0 4 + y 9Y )2 4 2 Z 4 L 9a ®4 cos ®4 


4 2 Z 4 L 9b (y 9Y 4 ®4 ) C0S ( ®4 + Sy> 


4 2 L 9a L 9b ®4 (y 9Y + ®4 ) C08 Y 9Y 4 2 Z 4YY ®4 


• l.i 

9YJ 2 


2 J 8YS ( ®4 + Y 8Y ) * 2 J 9YS ( ®4 * Y 9Y )2 * \ *4XX *4 


4 2 J 8XS ( *4 4 Y 8X) 4 2 J 9XS ( *4 4 Y 9X )2 
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V V - 0 


L ■ T - V - T 

• • 

There are 7 independent variables: Z^, 6^, Tgy» tgy, 

• t 

Ygx> y 9x* one equation for each of them. The generalized forces are 



similarly 
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d 

dc 


»4 Z 4 * "8 Z 4 ' "8 L 8. , 4 c °‘ *4 * " 8 L 8b ( »8Y * «4 ) co ' <e 4 * " r 8Y ) 


* “ 9 Z 4 * m 9 L 9. 6 4 C0B e 4 * m i L tb Sy * V cos 


(, 4 


Expand it, we have. 


"4 Z 4 4 “8 Z 4 ' m 8 L 8. e 4 co " 9 4 * "8 L 8« e 4 ,in 6 4 

- "e L 8b <t 8Y 4 V cos <e 4 4 ’by’ 4 Vsb (, 8Y 4 *4* 2 ‘ in <T 8Y* V 

•* •• *2 

♦ » 9 Z A ♦ m 9 L 9a®4 cos ®4 “ B 9 L 9a®4 sin ® 4 

* n 9 L 9b Sy 4 V COS <Y 9Y + V 

“ ®9 L 9b <t 9Y * ®4 )2 8in <Y 9Y + ®4 ) " F 4 (A,2) 


230 


Equation for 8^: 


d |3L \ 3 L 


i«6* | 


d_ 

dt 


”8 L 8a 8 4 * "s L 8b (8 4 * ' , SY ) ‘ m 8 Z 4 I '8a c0 ‘ *4 


* ”8 Z 4 L 8b cos <e 4 4 1r 8Y > * "8 l 8« l 8b T 8Y eos ’ 8Y 


* 2 ”8 L 8a L 8b e 4 co * T 8Y * “ l 9 I '9a # 4 * “9 L 9b <8 4 * T 9Y ) 


* “9 2 4 L 9a co ' 8 4 * V4 L 9b co8 (8 4 'Si 1 


+ m 9 L 9a L 9b Y 9Y COS y 9Y * 2 n 9 L 9a L 9b e 4 C0S y 9Y + hw* 4 


+ I 8YS <6 4 + ^Y* + X 9YS 


<8 4 * 1, 9Y ) ] 


•[ 


"8 2 4 I -8a 8 4 8i " *4 4 ”e 2 4 L 8b <T 8Y * V ,l ° ( V 1 8Y > 


B 9 2 4 l 9. 8 4 ,in 8 4 - "9 2 4 L 9b Sy * V * ln 


<8 4 4 ''9Y ) ]' 


l 4e 
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Expand it, we have, 

n, 8 L 8a e 4 4 “8 L 8b 6 4 4 m 8 L 8b Y 8Y " ®8 Z 4 L 8a C0S 6 4 


" m 8 Z 4 L 8b cos (6 4 4 y 8Y } 4 m 8 L 8a L 8b Y 8Y COS Y 8Y 


m 8 L 8a L 8b Y 8Y sin Y 8Y + 2 "8 L 8a L 8b e 4 C0S Y 8Y 


• • j •• 2** 2** 

2 m 8 L 8a L 8b e 4 Y 8Y Sin Y 8Y 4 m 9 L 9a e 4 4 m 9 L 9b e 4 + m 9 L 9b Y 9Y 


4 m 9 Z 4 L 9a cos 6 4 4 nl 9 Z 4 L 9b COS (0 4 4 Y 9Y } 


4 ID 9 L 9a L 9b Y 9Y C0S Y 9Y ‘ m 9 L 9a L 9b Y 9Y 610 Y 9Y 


4 2 ID 9 L 9a L 9b e 4 COS Y 9Y ~ 2 V^aSb® 4 Y 9Y Sln Y 9Y 


4 I 4YY 6 4 4 I 8YS 6 4 + I 8YS Y 8Y 4 I 9YS®4 + Z 9YS Y 9Y “ T 46 


Equation for y i 


8Y' 


d_ 

dt 



9L _ _ 

■■ B 7* 

81 8 Y 8V 


z .e. , 


dt |™8 L 8b <6 4 4 y 8Y 


> - V4 L 8b COS (6 4 4 Y W ) 4 ro 8 L 8a L 8b®4 COS Y 8Y 


♦ I 


8YS 


<6 4 4 Y Jv ) j“ [ m 8 Z 


4 L 8b (y 8Y 4 V 6in <6 4 4 Y 8Y ) 


m 8 L 8a L 8b 6 


8 4 <»»Y * ‘A 8 Sin Y 8 yJ 


8Y 
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Expand it, we have, 


m 8 L 8b 6 4 + m 8 L 8b Y 8Y " m 8 Z 4 L 8b C0S (6 4 + Y 8Y } * m 8 L 8a L 8b 6 4 COS 


+ I 8YS 6 4 + I 8YS Y 8Y + ir, 8 L 8a L 8b 6 4 sinY 8Y = T 8Y 


Equation for Y g : 


i.e. , 


d_ 

dt 



3 L 
SY 


9Y 


9Y 


d_ 

dt 


m 9 L 9b <6 4 + Y 9Y> + m 9 Z 4 L 9b COS <6 4 + Y 9Y> 


• • 


+ m 9 L 9a L 9b 6 4 COS Y 9Y + *9YS (6 4 + Y 9Y> 


' m 9 Z 4 L 9b (y 9Y + ®4 ) Sin (0 4 + Y 9Y } 


“ m 9 L 9a L 9b 6 4 (y 9Y + V sin Y 9Y “ T 9Y 


.1 

9Y 


Expand it, we have, 


V>b*« * ”9 L 9b> 9Y * V* L 9b cos ( ®4 * 'W 4 “9 L 9a L 9b®4 cos ' 


4 1 9YS® 4 4 I 9YS r 9Y 4 m 9 L 9a L 9b 9 4 * ln T 9Y " T 9Y 


8Y 

(A. 4) 


9Y 

(A. 5) 
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After linearization , Eqs. (A. 2), (A. 3), (A. 4), and (A. 5) become: 
Equation for Z^: 
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I 6 * 1 v ' * T 

9XS*4 a 9XS' 9X a 9X 
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APPENDIX B 


DEFINITION OF POSITIVE REALNESS AND STRICTLY POSITIVE 
REALNESS OF MATRICES [55] 


Definition B.l 

An m x m matrix H(s) of real rational functions i9 positive 
real if 

(1) All elements of H(s) are analytic in the open right half 
plane Refs] > 0 (i.e., they do not have poles with positive 
real parts). 

(2) The eventual poles of any element of H(s) on the axis Refs] 
= 0 are distinct, and the associated residual matrix of 
H(s) is positive semidefinite Hermltian. 

(3) The matrix H( jui)+H T (-j<j)) is a positive semidefinite 
Herraitian for all real values of w which are not a pole of 
any element of H(s)< 

Definition B.2 

An m x m matrix H(s) of real rational functions is strictly 
positive real if 

(1) All elements of H(s) are analytic in the closed right half 
plane Refs] > 0 (i.e., they do not have poles with non- 
negative real parts). 

(2) The matrix H(jw)+H T (-ju) is a positive definite Hermitian 
for all real to. 
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APPENDIX C 


DERIVATION OF FORMULAS FOR CALCULATING THE 

SPRING CONSTANTS AND DAMPING FACTORS OF DOCKING DEVICE 

Consider the angular springs and dampers first. Referring to 

the following figure, the moment of inertia I2 of the core body of 

the space station is used to represent the space station and I e 

s 

represents the shuttle 



the equations of motion can be written as 

1 2* e *2 + d A<® 2-V + K A( 0 2- e s> “ T 2 (C.l) 

I 2 9 8 + D A <V 0 2> + K A (0 8 -0 2 ) - T s (C.2) 


Assume that docking occurs at time t - 0, and 0 2 (O) «* 0 g (O) “ 0 (this 
will not affect the values of K A and D a ). Take Laplace Transforms, 
Eqs. (C.l) and (C.2) become 


( I 2 ® 2 + Da® ♦ K A > 9 2(®> " ( d A® + K A ) 0 s(®) “ T 2 (8) < c ‘ 3 > 

(V 2 + D a s + K A )0g(s) - (D a s+K a )0 2 (s) - T g (s) (C.4) 


238 


From Eq. (C.4) 


(D.s + K,)9 (s) 

V s > ' . 2 ~ ~ ~ .. + — 


V> 


Is + D.s + K. 
s A A 


I s + D.s + K. 
s A A 


(C.5) 


Substituting Eq. (C.5) into Eq. (C.3), we have, 


2 ( V +K A ) (D a 8+K A ) 

_ LIT \A/_\ ^ n / _ \ IK / _ \ < fi ** 


(1 2 s + V + V e 2 (s > - 2 


<V +D A 8+ V 


e 2 (s) = t 2 (s) + ^ 2 


<V +D A 8 +K A> 


T (s) 
6 


(C.6) 


1 • 0 • I 


(I 2 s 2 +D a s+K a ) (I 8 S 2 +D a s+K a )0 2 (s) - (D a s+K a ) 2 0 2 (s) 

= (I g s 2 +D A s+K A ) T 2 (s) + (D a s+K a ) T g (s) (C.7) 

For the purpose of obtaining the characteristics of the system, the 
right-hand side of Eq. (C.7) is set to zero. Then we have 

[I 2 I 8 s 2 + (I 2 +I 8 )D A s + (I 2 +I 8 )K a ]s 2 = 0 (C.8) 

for s f 0 ( *.’ s 2 - 0 implies two rigid modes at the origin) 

I 2 I s® 2 + (I 2 +i s> d A 8 + ( t 2 +I 8> k A " 0 (C.9) 
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or 



(C.10) 


(C.ll) 


(C.12) 


(C.13) 


(C.14) 


For the linear springs and dampers, 
by using the same approach 



and I>l are obtained as follows 

(C.15) 

(C.16) 


where M2 is the mass of the core body of the space station and M s is 
the mass of the shuttle. 
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ORIGINAL PAGE IS 

OF POOR QUAiiTX 


APPENDIX D 

PROGRAM LISTING FOR THE SIMULATION OF ADAPTIVE CONTROL DURING 
SHUTTLE HARD DOCKING TO FOUR-PANEL SPACE STATION 


1 PROCRAM 

2 

3 ' 

4 : 


APPLICATION OF ADAPTIVE CONTROL ON SPACE STATION (15 DOF) 
CREATED OCT 1984 —BY C.H.C.IH 


S 


6 


7 

INITIAL 

8 


9 


10 

s 

11 

t 

12 

0 

13 


14 

INTEGER 

15 

ARRAY 

16 

ARRAY 

17 

ARRAY 

18 

ARRAY 

19 

ARRAY 

20 

ARRAY 

21 

ARRAY 

22 

ARRAY 

23 

ARRAY 

24 

ARRAY 

25 

ARRAY 

26 

ARRAY 

27 

ARRAY 

28 

ARRAY 

29 = 

ARRAY 

30 

ARRAY 

31 

ARRAY 

32 

ARRAY 

33 


34 

CONSTANT 

3S 



isssssaasssS&sssa INITIAL =a=3s3ss3ss3&=as3sssaaassnssc3ss' 


************%**************** t**t****xx*tttx********x*** ' 

* TYPE AND DIMENSION OF VARIABLES *' 

**XX************X*X*X*X***X****X*XX**XX**X**XX*XX***XX* ' 

N,N2,M,L,L2,Q,II,J,K,NSTPi,NSTP2 

ZP < 15) ,ZPD( IS ) ,ZM(15) , ZMD ( IS), ZP 0(15) , ZPD0(15),ZPDD(15> 

AP (30,30) , BP (30, 11), CP (11, 30) 

AM(18,18) ,BM( 18,11 > ,CM(1 1,18) 

XP ( 30 ) ,XPD(3Q) , XP 0 ( 30 > ,XM ( 18 ) , XMD( 18) ,XM0(18) 

B(15,il) ,C (11 ,15) ,BPF( 15,11 > ,CPF(il, 15) ,BMF( 7,11) ,CMF (11,9) 
YP(ii) ,YM(ii) , UP < 11 ) ,UM (11) 

ETA( 15) ,ETA0(15) ,ETAD(1S> ,ETAD0(1S> ,ETADD(15> 

ETAS( 9 ) ,ETASD(9 > , ETAS0 ( 9) , ETASD0 ( 9 ) 

EVAL ( IS ) , EVEC ( IS , IS ) , EVECT (15,15), EVECI ( 15 , IS ) 

EVECS (15,9), EVECST (9,15), EVECSB (15,9) 

TA(40,40),TB(40 ,40) ,TAS(40> , TBS (40) ,R(40 ) ,RT(1,40) 

DAMPFC ( 15) , DAMP (1S,15),ZETA(9> , DAM (9 , 9 > , U2( 15 , IS ) ,U2M(9,9) 
KI(11,40),KP(11,40),KID(11,40),KI0(11,40),MASS(15,15) 

DUMY4 (1,40), DUMY5 (1,40) ,DUMY6 ( 11),MT(15,15),MTT(1S,15) 

DUMY7 ( 11 ) ,DUMY8(30 ) ,DUMY9<30 ) ,DUMYi 0 ( IB) , DUMYil ( 18) 

EYS(il) ,EY(ii,i),STIFF<15,15) ,K1(15,1S> 

FDDOCK (3) ,B1 ( 15 , 3 ) , BDF ( 15 , 3 ) , BD ( 30 , 3 ) 

DUMY12 (30 ) , DUMY13 (30 ) , SCR1 ( 15) ,SCR2( 15 ) ,FREQ(15) 

N=1S * ' NUMBER OF PHYSICAL COORDINATES' 

'(ALSO NUMBER OF MODAL COORDINATES)' 


36 = CONSTANT N2=30 
37 ; CONSTANT M=il 


30 = 
39 s 


40 

41 

42 

43 

44 

45 

46 


CONSTANT L=9 


CONSTANT 

CONSTANT 

CONSTANT 


L2*18 
UM=li*0 . 
TFINAL“i 0 0 . 


CINTERVAL CINT-0.2 


0 


$ 'NUMBER OF PLANT STATES' 

•'NUMBER OF PLANT INPUTS AND MODEL INPUTS 
'(ALSO NUMBER OF PLANT OUTPUTS AND MODEL 
' OUTPUTS) 

•'NUMBER OF THE MODES RETAINED(TQ FORM 
'THE MODEL) 

•'NUMBER OF MODEL STATES' 

•'MODEL INPUT' 

•'TIME TO STOP' 

•'COMMUNICATION INTERVAL' 


47= ' 

48= ' 

49= ' 


*******************************************************' 
« DEFINE THE PHYSICAL PARAMETERS OF THE SPACE STATION *' 

********************* 8 *********************************' 


50 = 

SI .CONSTANT 

52 = CONSTANT 

53 .CONSTANT 

54 = CONSTANT 

55 CONSTANT 

56 CONSTANT 

57 CONSTANT 

58 CONSTANT 


C»32. 174 
ZETA=9*0 .707 
EIS=9475776. 
LS«115. 
LE-140 . 

RHOS- . 541 
RHOE*l . 048 
M4-4165.35 


•'GRAVITATIONAL ACCELERATION' 

•'DAMPING RATIOS OF THE REFERENCE MODEL' 


S9 CONSTANT M8=994.72 
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60 CONSTANT 

M9-994 . 72 

61 CONSTANT 

I4XX-3 . 869E+6 

62 CONSTANT 

I4YY-1 343E+6 

63 : CONSTANT 

I8XS=2 . 437E+4 

64 ; CONSTANT 

I9XS=2 437E+4 

65 = CONSTANT 

1SYS-S . 637E+4 

66 ; CONSTANT 

19YS-S . 637E+4 

67 ; CONSTANT 

L3A=7 . 

60 : CONSTANT 

£88=11 • 

69 .CONSTANT 

L9A=7. 

70 CONSTANT 

L98=ii . 

71 ; CONSTANT 

FAC=10 . 

72 i CONSTANT 

P I»3 . 14159265 

73: 

L8=L8A+LBB 

74 = 

L9=L9A+L9B 

7S 

£IE=FAC*GIS 

76 : CONSTANT 

TAS=S*2 . SE+6, 2*2. 5E+8 ,4*2. SE+6, 18*1 000 . ,11*40 0 . 

77 : CONSTANT 

TBS=S*2 . 5E+6 ,2*2 . 5E+8 ,4*2 . 5E+6 , 18*1000 . ,11*400 . 

7D CONSTANT 

DAM PFC S 15*0 .005* 'MODAL DAMPING RATIOS OF THE PLANT' 

79 : CONSTANT 

ALPHA- 0.2 

SO : CONSTANT 

NSTP 1-200 

81 = CONSTANT 

NSTP2-50 

82. CONSTANT 

TSW=S. 0 

83: 

DTRCC=PI/180 .0 

04 

Q-M+L2+M 

8S: 

DO 40 11=1, N 

86 : 

DO 30 J=i,N 

87 : 

K1 ( II, J>-0 . 

88: 30. 

•CONTINUE 

89: 40. 

.CONTINUE 

90: 

Kl(7,7)«28184. 

91 : 

K 1 ( 8 , 8 ) -339324 . 

92: 

Ki <9 , 9 >-244308 . 

93- 

DO 60 JI»1,M 

94 : 

DO 50 1=1 ,Q 

9S = 

K10 ( II , J)-0 . 

96: 50. 

.CONTINUE 

97: 60. 

.CONTINUE 

98. 


99 = 

DO 62 11 = 1 ,N 

100: 

DO 61 1=1,3 

101 : 

81(11, J>«0. 

102: 61. 

.CONTINUE 

103: 62. 

.CONTINUE 

104 : 

81(7, 1)-1 . 

105: 

81(8, 2)=1 . 

106: 

Bl(9,3) = l . 

107: 

DO 65 11 = 1 , N 

108: 

DO 63 J-l ,M 

109 : 

B( II , J)-0 . 

110: 63. 

.CONTINUE 

111> 65. 

.CONTINUE 

112: 

B(2,i)=l. 

113: 

B( 4,2 )-l . 

114: 

B(5,3)*l- 

115: 

8(6,4 ) -1 . 

116< 

B(7,S)»i ■ 

117: 

B(8,6)-l . 

118: 

B( 9 ,7 ) — 1 
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119: 


B(10,8)«i. 


120 : 


B( 1 1,9 )*1 . 


121 : 


8(13,10 >»1 . 


122 : 


B< 15,11 )*i . 


123: 


DO 68 11*1, M 


124: 


DO 66 J=i , N 


1 25 : 


CCII,J) = 0 . 


126: 

66. 

.CONTINUE 


127: 

68. 

. CONTINUE 


128: 


C(i,2)=i. 


129: 


C <2 , 4) = 1 . 


130 : 


C < 3 ,5) =1 . 


131 : 


C(4,6)=i . 


132: 


C(5,7)»l. 


133: 


C(6,B)»i . 


134: 


C(7,9)-i. 


135: 


C(8,10>-1. 


136: 


C(9,il)»l . 


137: 


C( 10 ,13)*1 . 


133: 


C(il,iS)-l. 


139: 




140 : ' 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXX 

141 : ' 


* DEFINE THE DOCKING DISTURBANCE PARAMETERS 

X 

142: ' 


1 XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX*XX*XXXXXXX»XX«*XXXXXX 

143: 




144 ; CONSTANT 

ZDAM=3 . 06E+3 $ 'SOFT DOCKING LINEAR DAMPING COEFF' 


14S : CONSTANT 

ZSTI-1 . 3S97E+3 X'SOFT DOCKING LINEAR SPRING CONST' 


146 ; CONSTANT 

RDAMT*1 . 4S93E+7H 'SOFT DOCKING ANGULAR DAMPING COEFF' 


147 = CONSTANT 

RDAMP=7 . 078E+6 


148 : CONSTANT 

RSTIT=6 . 4845E+7X 'SOFT DOCKING ANGULAR SPRING CONST' 


149: CONSTANT 

RSTIP=3 . 14SE+7 


150 CONSTANT 

MSIiUT=25160 0 . X'MASS OF THE SHUTTLE' 


151 : CONSTANT 

ISIIUTT=7 . S4E+6 6 'MOMENT OF INERTIA OF THE SHUTTLE' 


152 ; CONSTANT 

ISHUTP-i. OE+6 


153: 




1S4 : ' 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX ' 

155. ' 


X CALCULATES M AND K MATRICES, THE EIGENVALUES AND 

X' 

156: ' 


* EIGENVECTORS OF THE OPEN LOOP SYSTEM, FORM THE 

X ' 

157. ' 


* DAMPING MATRIX IN MODAL COORDINATES , ASSEMBLE AP 

X ' 

1S8 . ' 


It BP AND CP MATRICES OF THE PLANT 

X ' 

159 : ' 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXtt ' 

160: 




161 : PROCEDURAL (MASS ,STIFF,EVAL , E VEC , DAMP .FREQ , AP , BP ,CP,BD*=. . . 


162: 


EIS, EIE , RHOS , RHOE ,LS,LE, M4 , M8 , M9 , I4XX , 14 YY ,I8XS , I9XS 

I • ' 

163: 


I8YS , I9YS ,L8A ,L9A ,L8B ,L9B,L8,L9,DAMPFC,B1 ) 


164: 


CALL RE ADMK <M AOS, STIFF.EIS.EIE, RHOS, RHOE ,LS ,LE , M4 ,M8 , 


165. 


M9 ,T4XX , T4YY , ISXS , I 9XS , ID YS , I 9YS , L8A , LV A , 


166- 


L8B,L9B,L8,L9> 


167. 


CALL MADD(STIFF, XI .STIFF, N,N) 


168. 


CALL OPLOOP (MASS, STIFF, EVAL.EVEC) 


169. 


DO 69 II«i,N 


170. 


SCR1 (I I )*ABS(EVAL( XI ) ) 


171. 


SCR2(II)-SQRT(SCR1<II> > 


172 = 


FREQ( II )«=SCR2( II )/ (2 . XP I ) 


173: 

69. 

.CONTINUE 


174. 


CALL MPRINT (MASS ,N ,N, 12 , 12H MASS MATRIX) 


175 1 


CALL MPRINT ( STIFF, N,N, 17, 17H STIFFNESS MATRIX) 


176: 


CALL MPRINT (EVAL ,i,N,19, 19H EIGENVALUES (WX*2) > 
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177 = 
173 = 
179 = 
100 = 
181 = 
182 = 
1 63 : 
184: 
185 : 
186: 
107: 
180: 

189 = 

190 = 

191 : 

192 = 
193. 
194: 
195: 
196: 
197: 
198: 

199 = 

200 = 
201 : 
202 = 

203 = 

204 = 

205 = 
206: 
207: 
208: 
209. 
210 = 
211 : 
212 . 
213 = 
214. 
215: 
216 = 
217 = 
218: 

219 = 

220 = 
221 : 
222 = 

223 = 

224 . 

225 = 
226. 

227. 

228. 
229. 
230 = 
231. 

232 = 

233 = 


CALL MPR INT ( EUEC , N , N , 20 ,20H MODAL MATRIX (EUEC)) 

CALL MPR INT( FREQ , 1 ,N , 23, 23H MODAL FREQUENCIES (HZ)) 

DO 80 11=1, N 
DO 70 J-1,N 
W2 ( II , J )*0 . 0 
70. . CONTINUE 

U2 ( II, II )*ABS(EVAL(I I ) ) 

80 . . CONTINUE 

CALL MPRINT(U2,N,N,43, . . . 

43H GENERALIZED STIFFNESS MATRIX OF PLANT (U2)) 
CALL MTRANS<EVEC,EVECT,N,N> 

CALL MMULT ( EUECT ,MASS» MT , N , N , N ) 

CALL MMULT (MT, EUEC ,HTT,N,N,N) 

CALL MPRINT(MTT,N,N,4,4H MTT ) 

DO 90 I 1=1 , N 
DO 90 J=i,N 
DAMP ( 1 1 , J ) = 0 . 0 
90 . CONTINUE 

DO 95 11*1, N 

DAMP(II,II)*2.0#DAMPFC(II)*SQRT(U2(II,II) ) 

95. . CONTINUE 

CALL MCOPY ( EVEC , EVEC I ,N,N> 

CALL INVERS(EVECI,N> 

CALL MPRINT(DAMP,N,N,5,5H DAMP) 

DO 110 11=1, N 
DO 110 J=i,N 
AP ( I I , J )=0 . 

110. . CONTINUE 

DO 120 II* 1 , N 
DO 120 J*1 ,N 
AP ( II , N+ J )=0 . 

IF(II.EQ.J) AP ( 1 1 , N+ J > = 1 . 

120 . . CONTINUE 

DO 130 11*1, N 

DO 130 J=1,N 

AP ( N+T I , J ) *-l . *U2(II,J> 

130. . CONTINUE 

DO 140 11*1, N 
DO 140 J= 1 , N 

AP(N+II,N+J)=-1 *DAMP(II, J) 

140 . . CONTINUE 

CALL MMULT (EVECT,B1 , 8DF,N ,N,3> 

CALL MMULT(EVECT,B,BPF,N,N,M) 

CALL MMULT < C , EVEC , CPF , M ,N ,N> 

DO 143 11=1, N 
DO 142 J*1 ,3 
BD(IX, J)*0 . 

BD(N+II, J)*BDF (II ,J) 

142. . CONTINUE 

143. . CONTINUE 

DO 145 11*1, N 
DO 144 J*1,M 
BP (II ,J)*0 . 0 
BP(NHI,J)-BPF(II,J) 

144. . CONTINUE 

145. . CONTINUE 

DO 148 11*1 ,M 
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234 


00 147 J=i,N 

233. 


CP (II, J)=ALPHAXCPF<IJ, J> 

236 - 


CP(ii,N+j>*cprcii, j> 

237: 

147. 

.CONTINUE 

238: 

148. 

CONTINUE 

237: 



24tl : END 


T- 'END OP PROCEDURAL' 

241 : 



242: ' 


XX*£*£XXXXXXXXX*X***X*XX*XXX*************#***X£SXXX*XX 

243: ' 


* ASSEMBLE AM , BM , CM MATRICES OP THE REFERENCE MODEL * 

244 . ' 


XXXX**XXX*»X**SXX*X»XXX*XXXXXXXX*«»X****XX«XXXXXXXX*XX 

24S : 



24 6: PROCEDUR AL.. ( AM , DM , CH , EVECS , EVECST .DAM , U2M-7ET A , E VEC , U2 , B , r , ALP 1 IA > 

247 = 


DO ISO 1 1 = 1 , L 

24f! . 


DO ISO J=1,L 

249: 


U2M( II , J)=0 . 

250 : 

150. 

.CONTINUE 

251: 


DO 1S2 11=1,6 

252 


W2h ( II , I I ) =W2 ( 6+II ,6+1 I )*0 .6 

253. 

152. 

.CONTINUE 

254: 


DO 1S3 11=1,3 

2SS: 


U2M(6+II ,6dI>=U2(12+II ,12+II>#1 .4 

2S6: 

153. 

.CONTINUE 

257: 


DO 160 11=1, L 

2S8: 


DO 155 J=1,L 

259: 


DAM (II , J ) =0 0 

260 : 

1.55. 

.CONTINUE 

261: 


DAM (II ,11 >=2 . 0*ZETA( II > *SRR T < U2M ( I I ,11 ) > 

262: 

160. 

.CONTINUE 

263: 


DO 170 11*1, L 

264: 


DO 170 J=1 , L 

265: 


AM ( 1 1 , J ) = 0 . 

266 : 

170. 

.CONTINUE 

267: 


DO 180 11 = 1 , L. 

268: 


DO 180 J=i , L 

269: 


Ah ( I I ,L+ J) =0 . 

270. 


IFdI.EO J) AM (II,L+J> = 1 . 

271 : 

180 . 

.CONTINUE 

272: 


DO 190 11=1, L 

273: 


DO 190 J=1 ,L 

274: 


AM (L+I I , I ) =-i . XW2M< II , J > 

27S: 

190 . 

.CONTINUE 

276: 


DO 200 11=1, L 

277: 


DO 200 J-l.L 

278: 


AM (L+I I ,L+J>=-1 . *DAM(II , J> 

279: 

200 . 

.CONTINUE 

280: 


DO 205 11 = 1, N 

281: 


DO 205 J*1 , 9 

282: 


EVECSB< II , J )=EVEC ( I 1,6+ J) 

203 = 

205. 

CONTINUE 

284: 


DO 212 J=1,L 

28S: 


EVECS( 1 , J)=EVECSD(i , J)*0 .7 

286: 


EVECSO, J)=EVECSB<2,J)#0 . 7 

287. 


EVECS (3, J ) =EVECSB<3, J>*1 . 3 

288: 


EVECS<4, J)«EVECSB<4,J)*1 .3 

289: 


EVECS<5,J>=EVECSB<5,J>*1 . 0 

290 = 


EVECS ( 6, J)=EVECSB<6 , J >#0.7 

291 : 


EVECS ( 7, J)=EVECSB(7 , J>#1 . 0 

292 


EVECS (8, J>=EVECSB(8,J>*1 . 0 
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293 = 

EVECS<9, J)*£VECSB<9, J)*l . 0 

294. 

EVECSflO , J>=EUECS8<10, J>*1 . 0 

295 

EVECC(il,J>=EVECSB<il,J>*i .3 

296 = 

EVECS<12,J>=EVECSBU2, J>*0.7 

277; 

EYECS ( 13, J )=EVECSB<1 3 , J >X0 .7 

298: 

EVEr.S<14,J>=EVECSB<14, J)X1 .3 

299: 

EVECS ( 15 , J)=EYECSB(1 5 ,J)XI . 3 

300: 212. 

•CONTINUE 

301 : 

CALL MTRANS<£YECS , EVECST,N, L ) 

302: 

CALL HMULT (EVECST ,B , BMF ,L , N ,h> 

303: 

CALL MMULT (C , EUECS,CMF , N,N,L > 

304: 

DO 215 11*1, L 

305: 

DO 215 J=1,H 

306: 

BhCII , J>*0 . 

307. 

IiM(L+II,J) *BMF ( II, J ) 

308: 215. 

CONTINUE 

309 = 

DO 220 Il-l.h 

310. 

DO 220 2*1, L 

311 : 

CM ( 1 1 , J)=ALPHA*CMF ( II, J) 

312: 

Ch< II ,L*J)*CMF<II,J> 

313= 220. 

.CONTINUE 

314: 


315: END 

•'END OF PROCEDURAL' 

316: 


317= ' 

«******t****t***»t**«***t»****«*x***:**K**x:*****s****X*** 

310: ' 

* SET THE INITIAL CONDITIONS IN PHYSICAL COORDINATES 

319: ' 

XX***************X*«*X**XX*X*XXXXXXXXXl*»*’.***»**«XX;t£*** 

320: 


321 "CONSTANT 

Z10=0 . 

322 "CONSTANT 

Ti 0*0 . 

323 CONSTANT 

230=0. 

324 = CONSTANT 

T30 = 0 . 

325 CONSTANT 

Z20=0. 

326 = CONSTANT 

P20=0 . 

327. CONSTANT 

Z40*0 . 

328 = CONSTANT 

T40-0. 

329: CONSTANT 

P40-0. 

330 CONSTANT 

260*0 

331 "CONSTANT 

P60*0 . 

332 ■ CONSTANT 

zso*o. 

333 ; CONSTANT 

T50-0 . 

334 "CONSTANT 

270=0. 

335 . CONSTANT 

T70*0 . 

336 "CONSTANT 

Z1D0-0. 

337 • CONSTANT 

TIDO-O. 

338 = CONSTANT 

23D0-0 . 

339. CONSTANT 

T3D0-0. 

340 . CONSTANT 

Z2DO-0. 

341 "CONSTANT 

P2D0=0 . 

34 2 "CONSTANT 

Z4D0-0. 

343 "CONSTANT 

T4D0-0. 

344 .CONSTANT 

P4DO-0. 

345 < CONSTANT 

26D0=0 . 

346 CONSTANT 

P6DO*0. 

347. CONSTANT 

25DO=0. 

348 .CONSTANT 

TSD0*0. 

349 = CONSTANT 

Z7D0-0. 
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350 

351 

352 

353 


CONSTANT 


354 

355 


356 

357 
356 

359 

360 


T7DO*0 . 

2P0 (1)=Z10 

ZPO <2)*T lOtDTRCC 

ZPO < 3 )*Z30 

ZP0C4)«T30*DTRCC 

ZP0(5)=Z2O 

ZPO (6)»P20*DTRCC 

ZP 0(7 > =Z40 

ZPO (8)«T40*DTRCC 

ZP0(9)*P40*DTRCC 

ZPO(10)«Z60 


361 


ZPOC 11)*P60*DTRCC 


362: 

363-. 

364. 

365. 

366. 

367. 

366. 

369: 

370. 

371 = 

372: 

373. 

374 = 

37S. 

376 = 

377. 

378 = 

379. 

380: 

381 CONSTANT 
302. CONSTANT 


ZP 0 (12 )=Z50 

ZP0(13>=T50tDTRCC 

ZPO ( 14 >*Z70 

ZP0(1S)-=T70*DTRCC 

ZPDOCi J-Z1D0 

ZPDO<2>«"T1DO*DTRCC 

ZFDQ(3>*Z3D0 

ZPDOC 4)*T3Dfl*DTRCC 

ZPD0(S)*Z2D0 

ZPDO(6)*P2DO#DTRCC 

ZP DO ( 7 )=Z4D0 

ZPD0(B)«=T4D0*DTRCC 

ZPDO(9)=P4DO*DTRCC 

ZPDOC 10) =Z6D0 

ZPDO(il)«P6DO*DTRCC 

ZPDO < 12)"Z5D0 

ZPDO C 13)*T5D0#DTRCC 

ZPDO ( 14 )*Z7D0 

ZPD0(1S)«T7D0*DTRCC 

DSHUi-0.2 

SHU1»0 . 


383 = CONSTANT 

384 . CONSTANT 

385. CONSTANT 
386 CONSTANT 
307. 

308< 

389: 

390. 

3V1: 

392.' 


DSHU2-0.2 

SHU2=0. 

ZDSHUOeO . 05 
ZSHUO-O . 

DSHUT0-DSHU14DTRCC 

SHUT0-SHU1*DTRCC 

DSHUP0«DSHU2*DTRCC 

SHUP0>SHU2*DTRCC 

*txx*t*t*t*****tt********xtx**t*tt****t*t#***********x 


393. ' * FIND THE INITIAL CONDITIONS IN MODAL COORDINATES X 

394. ' ***t**********«*tM*****«**»***t*******S**«»*******»** 

39S: 

396 . PROCEDURAL (XP 0 ,XMO»EUECI ,ZPO ,ZPDO ) 

397. CALL MVMULT CEVECI , ZPO , ETAO.N, N) 

398. CALL MVMULTCEVECI ,ZP DO, ETADO,N,N) 

399. DO 250 IZ»1,N 

400. XPO<II>»ETAO<II> 

401< XP 0 (N+II )*ETADO < II ) 

402. 250. .CONTINUE 

403. DO 258 II«1,L 

404 . ETASO ( II )«=0 .9XETA0 (6+1 1 ) 

405. ETASDO(IX) a O . 9XETAD0 (6*11 ) 

406. 258. .CONTINUE 

407. DO 280 XI-l.L 
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408 



XMOLII )*ETAS0 <11 > 

409 



XM0<L+II>=ETASD0<II> 

410 


200. 

CONTINUE 

411 



CALL HrRINT(XP0,i,N2,4,4H XP0J 

412 



CALL MPRINTLXM0 , 1 ,L2,4 ,4H XM0> 

413 




414 

END 


* 'END OF PROCEDURAL' 

415 




416 

0 


»«»*****X«****«*X*«X*****X«**X*****X*«**X*X***»**X****' 

417 

0 


% DEFINE THE SCALING MATRICES TA AND TB *' 

418 

0 


**«***X****************«***»***»************X*X**«****' 

419 




420 

PROCEDURAL 

-<TA,TB*TAS,TBS> 

421 



DO 290 II=i,Q 

422 



DO 290 J-1,Q 

423 



TALII, J>-0. 0 

424 


290. 

CONTINUE 

42S 



DO 300 11*1,0 

426 



TALII, II >oTAS(II> 

427 


300 . 

CONTINUE 

428 



DO 310 11-1,0 

429 



DO 310 J-l ,Q 

430 



TBLII, J)-0 . 0 

431 


310. 

CONTINUE 

432 



DO 320 1 1-1 , Q 

433 



TBLII, II)-TBS(II) 

434 


320 . 

CONTINUE 

435 




436 

END 


* 'END OF PROCEDURAL' 

437 




438 




439 

END 



440 




441 

•wmm 

X 3C8BS 

iBatacBiuBSM* END OF INITIAL ■■•naB»aBBBBoaasMS*«BSBs ' 

442 




443 




444: 




44S 




446> 




447 

' mmm 

CSSSBI 

ISItBCBBBBBSSSCBBSa DYNAMIC BSBBBBBBBBBBSBB&SailBOBBBSSr SC 

448 




449 

DYNAMIC 


450 



TERMTLT . GT . TFINAL ) 

451 



TERMT <XP < 1 J . CT . 1 . 0E+1 0 ) 

452 



TERMT<XP<2).GT.i.0E+10> 

453 

PROCEDURAL (NSTP-T, TSU ,N5TP 1 .NSTP2 ) 

454- 



IF(T.GT.TSW) GO TO 325 

455 



NSTP-NSTPi 

456 



GO TO 327 

457. 


325.. 

NSTP -N5TP 2 

458> 


327.. 

CONTINUE 

459 

END 


•'END OF PROCEDURAL' 


460: 

461 s DERIVATIVE 

462: 

463: DERIVATIVE 
464: 
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465 

466 

467 

468 

469 

470 

471 

472 

473 

474 

475 

476 

477 

478 
479: 
400 
481 
4B2 

483 

484 

485 

486 

487 

488 

489 

490 

491 

492 
493: 

494 

495 

496 

497 

498 

499 

500 

501 
502: 

503 

504 

505 

506 

507 

508 

509 
510. 
*11. 
S12 
513. 
514 : 

515. 

516. 

517. 
S18> 

519. 

520. 
521 : 
522 


**(********** *******(* ******* ***************** **r £**»**«* 

' * CALCULATE THE SOFT DOCKING DISTURBANCE FORCE & TOROUC * 

' X******************************************************** 

PROCEDURAL (FDOCK ,TQDCKT,TQDCKP .FDDOCK .DDSHtJT , DDCHUP .ZDDSHU- . . . 

R DAMT , RDAMP ,ZDAM .DSHIIT , DSHUP , ZDSHU , ZPD , RSTIT, . . . 

R STIP , ZSTI , SHUT , SHUP ,ZSHU , ZP , ISHUTT , ISHUTP , MSHUT , G ) 
TQDCKT-RDAMT* < DSHUT-ZPD < 8 > > +RSTIT* < SHUT-ZP (8 ) ) 

TRDCKP -RDAMP* < DSHUP -ZPD < 9 > > *RSTIP*< SHUP-ZP < 9 > > 

FDOCK =ZD AM* (ZDSHU-ZPD < 7 ) >*ZDTI* (ZSHU-ZP (7) ) 

DDSMUT =- ( TQDCKT/ISHUTT > 

D DSHUP --(TQDCKP/ISHUTP > 

ZDD5HU— < FDOCK/ <MSHUT/G > ) 

FDDOCK (1) -FDOCK 
FDDOCK (2 J-TQDCKT 
FDDOCK (3>-TQDCKP 

END *'END OF PROCEDURAL' 

' **X*****************************************X*********' 

♦DEFINE THE INTEGRAL GAIN KI AND PROPORTIONAL GAIN KP*' 

' *y*«******x***x*«***X*X*X********:****** *************** - 

PROCEDURAL(KI , KP-YM.YP ,XM ,UM ,TA ,TB ,KI 0 > 

DO 330 II-l , M 
EYS<II)-YM<II)-YP<II> 

330. . CONTINUE 

DO 340 II-i.M 
R(II)-EYS<II) 

340 . . CONTINUE 

DO 350 II-1.L2 
R(M+II)-XM<II> 

350. . CONTINUE 

DO 360 II«1,M 
R < M+L2+ 1 1 > -UN (II) 

360.. CONTINUE 

' DO 365 II=i,H 

EY<II,i)=EYS<II) 

365. . CONTINUE 

CALL UTRANS<R,RT,8> 

CALL HMULT<RT,TA,DUMY4,i,8,Q> 

CALL MMULT<RT,TB,DUMY5,1 ,Q,Q> 

CALL MMULT (EY .DUMY4 , KID ,M , 1 ,Q > 

CALL MMULT ( EY , DUM Y5 , KP , M , 1 , Q ) 

' KI«INTVC(K ID ,K 10 > 

END * 'END OF PROCEDURAL' 

: ' ****************************************************** ' 

' * FIND THE CONTROL UP FOR THE PLANT *' 

>' ************X********************************7X*******' 

PROCEDURAL (UP-K I ,KP,R I 

CALL MVMULT<KI , R .DUMY6 , M , Q) 

CALL MVMULT(KP,R,DUMY7,M,Q) 

CALL VADD<DUNY6,DUMY7.UF,M> 

END • 'END OF PROCEDURAL' 
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523 

524 : ' 

525 = ' 
526 . ' 
527 = 
520 = 

529 = 

530 = 

531 = 

532 = 

533 = 

534 : 

535 = 

536 = 

537 = 

538 = 

539 . 

540 = 

541 = 

542 = 

543 = 

544 = 

545 = 

546 = 

547 = ' 


***************** ************************************* ' 
* INTEGRATE THE EQUATIONS OF THE PLANT AND THE HODEL *' 
******************************************************' 

CALL MVMULTt AP , XP ,DUMY8,N2 ,N2> 

CALL MVMULT ( BP , UP , DUMY9 , N2 , M > 

CALL MVMULT (BD,FDDOCK, DUMY 12 , N2 , 3 ) 

CALL VADD < DUM Y8 , DUMY9 , DUMY 1 3 , N2 ) 

CALL VADD < DUMY13 ,DUMY12,XPD, N2 ) 

XP«INTVC<XPD,XPO> 

CALL MVMULT < CP , XP , YP ,M ,N2 > 

CALL MVMULT < AM,XM , DUMY1 0 ,L2 , L2 > 

CALL MVMULT(BM,UM,DUMYli ,L2,M> 

CALL VADD (DUMY 10 , DUMY1 1 , XMD , L2> 

XM~ INTVC< XhD ,XMO) 

CALL MVMULT ( CM , XM , Yh , M , L2 ) 

DSIIUT* INTEG (DDSHUT , DSHUTO ) 

DSHUP»INTEG < DDSHUP ,DSHUPO > 

SHUT=INT£G (DSHUT , SHUT 0 > 

SHUP“ INTEG (DSHUP ,SHUPO > 

2DSHU= I NTEG ( ZDDSHU , 2DSHU0 > 

ZSHU= INTEG (ZDSHU , 2SHU0 > 

******************^* *********************** *********** ' 


548 : 

S49: 

550 

551 

552 

553 

554 

555 

556 
557 . 
558 = 
559 . 

560 . 

561 = 

562 = 

563 = 

564 

565 

566 

567 
S6Q 

569 

570 

571 

572 . 

573. 

574. 

575. 

576. 

577. 
578 = 

579. 

580. 
581 = 


' * OBTAIN RESULTS IN PHYSICAL COORDINATES *' 

' ****************************************************** ' 

PROCEDURAL < ZP , ZPD.ZPDD, ZM, ZMD=XP , XPD , XM ,EVEC > 

DO 370 11*1, N 
ETA<II>*XPUI) 

ETAD (I I ) *XP <N+ II) 

ETADD<II)*XPD<N+II> 

370. . CONTINUE 

CALL MVMULT (EVEC,ETA,ZP ,N,N> 

CALL MVMULT <EVEC,ET AD, ZPD,N,N> 

CALL MVMULT < EVEC , ETADD,ZPDD ,N ,N> 

DO 380 11*1 , L 
ETAS< 1 I )=XM< II > 

ETASDdI )*XH(L+II> 

380. . CONTINUE 

CALL MVMULT < EVECS , ETAS , ZM , N , L > 

CALL MVMULT (EVECS , ETASD , ZMD , N ,L ) 


END 


END 


•'END OF PROCEDURAL' 


KSBBBBBSBBBBEBaB END OF DERIVATIVE BBBBBBBBBBBBB BBBB bbbbbbb 

****************************************************** ' 

* CHANGE RADIANS TO DEGREES *' 

* *********************** *********************** ******* ' 

ZlP-ZP(l) 

TIP «ZP (2) /DTRCC 
Z3P«ZP<3> 
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58;: = 

T3P=ZP < 4 ) /BTRCC 

583 : 

Z2P=ZP (S ) 

584 : 

P2P=ZP(6)/DTRCC 

505: 

Z4P=ZP (7) 

586: 

T 4P =ZP ( 8 > /BTRCC 

507: 

P4P=ZP(9)/BTRCC 

588: 

Z6P=ZP< 10) 

589: 

P6P=ZP (11 ) /BTRCC 

590: 

Z5P=ZP ( 12) 

591 . 

T5P=ZP < 13) /BTRCC 

S92: 

Z7P-ZP< 14) 

593 

T7P=ZP < 15)/DTRCC 

594 : 

Z1BP“ZPB < 1 ) 

595: 

T1DP»ZPB<2)/DTRCC 

596: 

Z3DP*>ZP0 C 3) 

597 : 

T3DP“ZPB (4 )/DTRCC 

590 : 

Z2I)P«ZPB<5) 

599: 

P2DP=ZPB<6)/BTRCC 

600 : 

Z4DP=ZPB< 7 ) 

601 . 

T4DP=ZP 0(8) /BTRCC 

602. 

P40P=ZPB ( 9 )/DTRCC 

603: 

Z6DP=ZPB< 10 ) 

604: 

P6BP“ZPB< 1 1 ) /BTRCC 

60S: 

2SBP=2PB< 12 ) 

606 

TSDP=ZP0<13)/DTRCC 

607: 

Z7BP-ZPB<14) 

608: 

T7DP*ZPB < IS) /DTRCC 

609: 

Zlh«ZM(i) 

610 : 

T 1M=ZM<2) /DTRCC 

611 : 

Z3M=ZM<3) 

612 

T3M*Zh < 4 )/BTRCC 

613: 

Z2M=ZM<S) 

614: 

P2M=ZM<6)/BTRCC 

615- 

Z4«=ZM<7> 

616< 

T 4M=Zh( 8 ) /BTRCC 

617. 

P4M=ZM<9 ) /BTRCC 

6 1 0 * 

Z6M*ZM< 10) 

619: 

P6M"ZM< 1D/DTRCC 

620: 

Z5H«=ZM< 12) 

621: 

TSM=ZM( 13)/BTRCC 

622: 

Z7M»ZM<14) 

623 : 

T7M*=ZM ( 15 )/BTRCC 

624: 

Z1DM«ZMD<1> 

625: 

T1BH-ZMB(2)/BTRCC 

626: 

Z3Bh=ZMD(3> 

627: 

T3DM-ZMB(4)/BTRCC 

620. 

Z2DM=ZM0<5> 

629: 

P2DH«ZMB<6)/DTRCC 

630 < 

Z4BM=ZMB<7 ) 

631. 

T 4DM=ZMBC 8 )/DTRCC 

632: 

P4DM*ZMB<9)/DTRCC 

633. 

Z6BM=ZMB(10) 

634: 

P6BM-ZHB(il)/BTRCC 

635: 

Z5DM«ZHB(12) 

636: 

T5BM«ZND C 13)/BTRCC 

637: 

Z7BM*ZMB< 14) 

638. 

T7DM«.ZND( IS) /BTRCC 

639: 
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t>40 . 
641 = 
642: 
643 = 
644. 

645 = 

646 = 

647 = 

648 • 

649 = 

650 = 

651 = 
632 = 

653 = 

654 = 

655 = 

656 = 

657 = 

658 = 

659 = 

660 = 
661 = 
662 < 

663 = 

664 = 

665 = 

666 = 

667 = 

668 = 

669 = 

670 = 

671 = 

672 = 

673 = 

674 = 

675 = 

676 = 

677 = 

678 = 

679 = 

680 = 
681 = 
682 = 

683 = 

684 = 

685 = 

686 = 

687 = 

688 = 
689 * 
690 < 
691 < 
692 > 
693 . 

694 = 

695 = 

696 = 


0SHANT-DSHUT /DTRCC 
SHANT-SMUT/DTRCC 
DSHANP-DSHUP/DTRCC 
SHANP-SHUP/BTRCC 
XP1S=XP(1 ) 

XP14-XP (2) 

XP13-XP (3 ) 

XP12-XP (4) 

XPil=XP (5 > 

XPiQ-XP (6) 

XP9=XP (7) 

XP8-XP (8) 

XP7=XP(9> 

XP6-XP(10> 

XPS-XPUl) 

XP4-XP < 12) 

XP3-XP (13) 

XP2-XP ( 14) 

XPi-XP (IS) 

XM9-XM ( 1 ) 

XH8=XM(2) 

XM7«XM(3) 

XM6-XM<4> 

XMS-XM(5) 

XM4=XM(6) 

XM3-XM(7> 

XM2=XM(B> 

XM1-XM(9) 

YP 1-YP < 1 ) 

YP2-YP <2 ) 

YP3-YP (3) 

YP 4-YP (4 ) 

YPS-YP <S) 

YP6-YP (6) 

YP7-YP (7) 

YP8»YP(B) 

YP9»YP(9> 

YPIO-YP(IO) 

YP1 1-YP (11) 
YMi-YM(l) 

YM2*YM(2) 

YM3-YM ( 3 ) 

YM4=YM(4> 

YHS-YM(5) 

YM6-YM (6) 

YH7-YM(7> 

YMB=YM(8) 

YH9-YM(9) 

YM10-YMU0) 

YMil-THdi) 

UP1-=UP(1> 

UP 2-UP (2) 

UP3=UP (3) 

UP 4-UP (4) 

UP5-UP <5> 

UP 6-UP (6) 


697: UP7«UP <7) 

698: UP8«UP(8) 

699: UP9«UP(9> 

700-. UPIO-UP(IO) 

701.- UPll*UP(il> 

702 = 

703 :END 
704 = 

70S: / «SMssBscasBcsrs»BBSRBssa END OP DYNAMIC >B3Sss BasaasasessaaBuss ' 
706: 

707 = 

708: 

709: ' ■sesssccauttsciiBsatsBsiaMi TERMINAL sssssasBiasBaeiasBiaiicsBsasa ' 
710: 

711 TERMINAL 
712: 


713: ' 

714: ' 

715: ' 

716: 

717: 

718: 

719: 

720: 

721: 

722 = 

723 = 
724: 
725: 

726 = 

727 = 
728. 


****************************************************** ' 

* PRINT SYSTEM PARAMETERS *' 

*!M************»***********«******«*****************X*' 

CALL MPRINT<EVECT,N,N,34, . . . 

34 H TRANSPOSE OF MODAL MATRIX (EVECT) ) 

CALL MPRINT(EVCCI ,N,N,32, . . . 

32H INVERSE OF MODAL MATRIX (EVECI ) ) 

CALL MPRINT(EVECS,N,L,3i, . . . 

31 H TRUNCATED MODAL MATRIX (EVECS)) 

CALL MPRINT( DAMP , N , N , 35, . . . 

3SH DAMPING MATRIX OF THE PLANT (DAMP)) 

CALL MPRINT(DAM,L,L,34, . . . 

34H DAMPING MATRIX OF THE MODEL (DAM)) 

CALL MPRINT(U2,N,N,43, . . . 

43H GENERALIZED STIFFNESS MATRIX OF PLANT (U2)> 


729 = 
730. 


CALL MPRINT(U2M,L,L,44 , . . . 

44H GENERALIZED STIFFNESS MATRIX OF MODEL ( L‘2M ) ) 


731: 

732: 

733: 

734: 

735: 

736. 

737 = 

738 = 

739. 

740. 
741: 
742. 

743 = 

744 = END 
745. 

746 : 
747. 
740. 

749 = 
750: END 


CALL MPRINT (AP ,N2 , N2, 1 0 , 1 OH MATRIX AP) 
CALL MPRINT (BP ,N2, M, 1 0 , 10H MATRIX DP) 
CALL MPRINT (CP ,M ,N2, 10 , 1 OH MATRIX CP) 
CALL MPRINT(AM,L2 ,L2 , 1 0 , 10H MATRIX AM) 
CALL MPRINT (BM ,L2 ,M , 1 0 , 1 OH MATRIX BM> 
CALL MPRINT (CM ,M,L2, 1 0 , 1 OH MATRIX CM) 
CALL MPRINT (TA,Q,Q,10,10H MATRIX TA> 
CALL MPRINT(TB,Q,Q,10,10H MATRIX TB) 


•mmmmmmmmwmmmm END OF TERMINAL 


VEND OF PROGRAM' 
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1 : 

2 = 

3 : 

4 = 

5 • 

6 ■ 

7 : 

8 
9 : 

10 : 
11 : 
12 : 
13 : 

14 : 

15 = 
16 : 
17 : 
18 : 
19 : 
20 > 
21 : 
22 : 

23 = 

24 : 

25 : 
26 : 

27 = 

28 = 
29 : 

30 : 

31 = 
32 : 
33 . 
34 = 
35 > 
36 = 
37 . 

38 = 

39 = 

40 : 

41 = 
42 : 
43 : 

44 = 

45 : 

46 = 

47 = 
48 : 
49 = 

so > 
Si: 
S2 = 
53 . 

54 = 

55 = 
56 : 
57 . 
SB: 
59 = 


SUBROUTINE READMK < MASS , STIFF ,EIS , EIE ,RHOS ,RHOE ,LS ,LE , M4 , 
*M8 ,M9 ,I4XX ,I4YY , IOXS , I9XS , ISYS , I9YS.L8A ,L9A,L8B ,L9B , 
*L8,L9> 

REAL MASS < 1 5, 15) .STIFF <15,15 > ,MC< 15, IS) , MIX 15, 15) 

REAL I4XX, I4YY.IQXS, I9XS, I8YS , I9YS ,L8A , LBB , L9A , L9B ,L8 , L9 
REAL ALPHA, BETA, A, B.EIS, EIE, RHOS , RHOE , LS ,LE , M4 , M8 , M9 
ALPHA-2 . *EIS/<LS**3 ) 

BETA-2. *EIE/<LE**3) 

A=RHOS*LS/420 . 

B-RHOE*L£/420 . 

00 10 1=1,15 
DO 10 J=i ,15 
STIFF< I, J)=0 . 

10 CONTINUE 

STIFF<i,l)=6.*ALPHA 
STIFF! 1 ,2 )«3 . tLStALPHA 
STIFF < 1 ,S)=-6 . tALPHA 
STIFF <1 ,8)=3 . tLStALPHA 
STIFF<2,2)=2 . t(LStt2) tALPHA 
STIFF < 2,5 )*~3 .tLStALPHA 
STIFF <2 , 8 )=<LStt2) tALPHA 
STIFF< 3 ,3 ) =6 . tALPHA 
STIFF<3,4)=-3. tLStALPHA 
STIFF<3,5>=-6. tALPHA 
STIFF <3, 8 )*-3 . tLStALPHA 
STIFF <4 ,4 ) =2 .t <LStt2) tALPHA 
STIFF<4 ,5 )=3 . tLStALPHA 
STIFF < 4 , 8 ) = <LStt2 ) tALPHA 
STIFF<S,5)*12. tALPHA+6.fBETA 
STIFF < S,6 )*3 . tLEtBETA 
STIFF(S,7 )=-6 . t£<ETA 
STIFF <5, 9) =3. tLEtBETA 
STIFF(6,6)=2.t<LEtt2)tHETA 
STIFF <6,7)=-3. tLEtBETA 
STIFF (6, 9 )= (LEtt2)tBETA 
STIFF <7,7)=12. tBETA 
STIFF <7, 10 )=-6. tBETA 
STIFF <7, 11 )»3. tLEtBETA 
STIFF<8,8)«B.t<LStt2) tALPHA 
STIFF ( 8, 12 >-3. tLStALPHA 
STIFF <8, 13)*(LStt2)tALPHA 
STIFF (8,14 )=-3 . tLStALPHA 
STIFF<8,15)=<LStt2)tALPHA 
STIFF < 9, 9 > =4 . * < LEttOtBETA 
STIFF<9, 10 )»-3 . tLEtBETA 
STIFF <9,11 )-<LEtt2)tBETA 
STIFF < 10 , 10 >=12 . tALPHA+6 . tBETA 
STIFF < 10 , 1 1 )=-3 . tLEtBETA 
STIFF! 1 0 , 12>«-6 . tALPHA 
STIFF < 10 , 1 3)«-3 . tLStALPHA 
STIFF<10,14>«-6. tALPHA 
STIFF <10 , 1 5>»3 . tLStALPHA 
8TIFF<li , 11 )«2. t(LEtt2)tBETA 
STIFF < 12 , 12>«6 . tALPHA 
STIFF <12, 13)=3. tLStALPHA 
STIFF <13, 13)«>2.t<LStt2> tALPHA 
STIFF < 1 4 , 14 ) »6 . tALPHA 
STIFF <14, 15)»-3. tLStALPHA 
STIFF! 15, 1S>=2. t<LStt2)tALPHA 
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60 = 

DO 20 1-1,14 

61 : 

K-I+i 

62 . 

DO 20 J-K , 15 

63: 

STIFF< J,I)-STIFF(I, J) 

64= 20 

CONTINUE 

65> 


66 : 

DO 30 1-1,15 

67: 

DO 30 J-1,15 

6B: 

MC(I , J>=0 . 

69: 30 

CONTINUE 

70 : 

MC(1 ,1>-1S6.*A 

71 : 

MC(1 ,2>-22.*LS*A 

72: 

MC<1 , 5 >=54 . <A 

73 = 

MC(i,8)»- 13.XLSXA 

74: 

MC(2,2)»4.*(LS**2>*A 

75 : 

MC(2,S>-13.*LS*A 

76 : 

MC(2,B)«- 3 *(LS**2)*A 

77: 

MCC3,3)«156.*A 

78: 

MC(3,4>— 22.XLSXA 

79: 

MC(3,5)-S4. *A 

80 : 

MC(3,Q)-13.XLS*A 

81 : 

HC(4,4)-4.*(LS**2)«A 

82 : 

MC<4,5) — 13.XLSXA 

83: 

MC(4,8)— 3.*(LS**2)*A 

84> 

MC(5 , S >-312 . *A+ 1S6XB 

85. 

HC(S,6)«22.*LE*B 

86 = 

MC(5 , 7 )-54 . *B 

87: 

HC(S,9)— 13.XLEXB 

88: 

MC(6,6)«4.*(LE**2)*B 

89: 

MC(6,7)-13. *LE*B 

90 : 

MC(6,9)=- 3. *(IE**2)*8 

91 = 

MC(7,7)-312.«B 

92: 

MC(7,10)-54.*B 

93: 

MC(7,il>— 13.XLEXB 

94: 

MC(8,8)-16.*(L5«2)*A 

95: 

MC(8,12)— 13.XLSXA 

96: 

MC(B,13>— 3.«(LSM2)*A 

97: 

MC(8,14)-13.*LS*A 

98: 

MC<8,15)— 3.«(LS»*2)*A 

99: 

MC(9,9)-8.»(LE**2>*B 

100 : 

MC(9,10>— 13 . *LE*8 

101: 

MC(9,li>— 3.*(LE**2)*B 

102: 

MC(10 , ID )— 156. XB+312 . SA 

103: 

MC(10,11 )--22.*LE*B 

104 = 

MC(10, 12)^54. *A 

105 = 

HC(10,13>-13.*LS*A 

106: 

NC (10,14)— 54. XA 

107: 

MC ( 1 0 , IS )--13. XLSXA 

108: 

HC(11 ,li)-4.X<LEXX2)XB 

109. 

HC(12,12)-156.«A 

110: 

HC< 12,13 >-22. XLSXA 

111> 

HC (13, 13) -4.X LSXX2 > XA 

112< 

HC( 14 , 14 >-156 . XA 

113: 

MC(14,15>— 22. XLSXA 

114: 

HC(1S,15>-4.X<LSXX2>XA 

115: 

DO 40 1-1,14 

116> 

K-I+1 

117. 

DO 40 J-K, 15 

118: 

HC ( J , 1 ) -MC ( I , J ) 
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119 

120 
121 
122 

123 

124 

125 
12 6 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 


40 CONTINUE 

DO 50 1-1,15 
DO SO J-i , 15 
HD < I , J) - 0 . 

SO CONTINUE 

HD<7,7>=M4+M8>M? 

MD<7,8)=M9*L9-MQ*L8 

MD<8,8> = J4YY+Mn*<l_S**2>*M9*(L9**2> + IBYS«I9YS 
HD(9,?) = I4XX-» I8XS+I9XS 
DO 60 1 = 1 ,14 
K=I*i 

DO 60 J=K , 15 
MD<J,I)=HD(I,J) 

60 CONTINUE 

DO 70 1-1,15 
DO 70 J=1 ,15 

MASS ( I , J ) = NC ( I, J)+MDCI , J> 

70 CONTINUE 
RETURN 
END 


SUBROUTINE OPLOQP (MASS, STIFF, EVAL, EVEC) 


1 : 

2 = 
3C 
A . C 
SC 
6 i C 
7.C 
8 : C 
9.C 
10 C 
11C 
1 2 : C 
13 C 
14:C 
IS i C 
16. C 
17C 
1 3 . C 
19. C 

20 C 

21 C 

22 . C 

23 = C 
24. C 
2S . C 

26. C 

27. C 

28. C 

29. C 

30. C 
31 . 

32. 

33. 

34. 

35 . 
36 . 

37. 

38. 

39. 

40. 
41 . 
42. C 
43 C 
44. C 
4S . 
46< 

47. 

48. 

49. C 
SO . 

51. 

52. 

53. 
54 > 
SS. 
S6C 
57. 
58 . 
59. 

60 .C 

61 C 

62 
63. 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk* 


k k 

k MASS AND STIFF ARE THE INPUT H AND K MATRICES. k 
k SUBROUTINE FIRST FINDS SIMILARITY TRANSFORMATION PHIl SUCH * 
k THAT THE COLUMNS OF PHI1 ARE EIGENVECTORS OF HASS NORMALIZED k 
k SUCH THAT k 
k k 
k (PHIl ) ' (MASS ) (PHIl )=( IDENTITY) (' DENOTES TRANSPOSE) k 
k k 
k THEN FINDS ANOTHER SIMILARITY TRANSFORMATION PHIS NORMALIZED k 
k SUCH THAT * 
k k 
k (PHI2) / (PHI1)MSTIFF><PHI1XPHI2>«DIAG(EVAL) k 
k k 
k WHERE EVAL ARE THE EIGENVALUES OF THE OPEN-LOOP SYSTEM, * 
k I.E., THE SQUARES OF THE EIGENFREQUENCIES . * 
k k 
k FINALLY, THE PRODUCT k 
k » 
k EVEC=(PHIi) (PHI2) k 
k k 
k IS RETURNED - THE J COLUMN IS EIGENVECTOR OF OPEN- LOOP k 
k SYSTEM CORRESPONDING TO J EIGENVALUE. k 
k k 
k (UNTIL THE SECOND CALL TO SYMQRR, PHI2 IS USED FOR k 
k INTERMEDIARY CALCULATIONS) * 
k k 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkAkkk 

PARAMETER N=15 

REAL MASS < N ,N) 

DIMENSION STIFF<N,N> 

DIMENSION EVEC(N,N),EVAL(N) 

DIMENSION PHI1(N,N) ,PHI2(N,N) , PHIl T< N,N> , WORK (500 ) 


CALL MCOPY( MASS, PHIl , N ,N ) 

CALL SYMQRR (*601, PHI 1,N,N, EVAL, UORK) 

RENORMALIZE FIRST SIMILARITY TRANSFORMATION 

DO 100 I«i,N 
DO 100 J=i,N 

PHI1(I,J)=PHI1(I,J ) /SQRT (EVAL ( J ) ) 

100 CONTINUE 

CALL MTRANS(PHI1 ,PHI1T,N,N> 

CALL MMULT<STIFF ,PHI1,PHI2,N,N,N) 

CALL MMULT (PHI 1T,PHI2,EVEC, N,N,N> 

CALL KCOPY (EVEC,PHI2,N,N) 

CALL SYMQRR(*701 ,PHI2,N,N, EVAL, WORK) 

EVAL NOU CONTAINS THE SQUARES OF THE E I GCNFR EQUENC IES 

CALL MMULT (PHIl, PHI2 ,EVEC ,N,N,N) 

EVEC NOU CONTAINS THE CORRESPONDING EIGENVECTORS AS 
COLUMNS 

CO TO 900 
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64 = 

65 = 

66: 601 CONTINUE 

67= WRITE(6,602) 

60: 602 FORMAT (IX, 'CLUTCHED ON FIRST CALL TO SYMQRR ' ) 

69: CO TO 900 

70 = 701 CONTINUE 

71: WRITE(6,702) 

72: 702 FORMAT (IX, 'CLUTCHED ON SECOND CALL TO SYMQRR') 

73: 900 CONTINUE 
74: RETURN 

75= END 


1: SUBROUTINE INUERS(A.N) 

2 

3: DIMENSION A(N,N> .WORK (1000) 

4 . 

5 : C THIS SUBROUTINE WILL CALL AINU FOR AN ACSL PROGRAM 

6 : C AND IT WILL HANDLE THE ERROR TRAPPING 

7. 

0: CALL AINUR (A,N,N,S100 , WORK ) 

9: 

10 > RETURN 

11 . 

12 = 

13 :C HANDLE ERROR BOMB 0UT8 

14: 100 CONTINUE 

15: URITE(6,1000) 

16: 1000 FORMAT ( IX , ' ************* INUERSION PROBLEM ft******#**##**#*') 

17: 

18: RETURN 0 

19: END 
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1. SUBROUTINE MADD < A, B ,C,L ,N> 

2 C THIS SUBROUTINE CALCULATES THE SUM A+B AND STORES THE RESULT IN C 

3 -C L = ROUS OF A AND B , N » COLS OF A AND B 

4-- DIMENSION A(L,N),B<L,N),C<L,N) 

S: DO 100 1*1 ,L 

6: DO 100 J=i,N 

7. C<I,J> - A<I, J)+B(I, J> 

8 100 CONTINUE 
?: RETURN 

10: END 


1= SUBROUTINE MMULT1A , B,C,L ,M,N> 

2 ; C THIS SUBROUTINE CALCULATES THE PRODUCT AB AND STORES THE RESULT IN C 

3 * C L * ROUS OF A , M * COLS OF A AND ROWS OF B, N = COLS OF B 
4- DIMENSION A (L ,M ) , 8< M , N ) , C<L , N) 

5: DO 100 1*1, L 

6: DO 100 J=i,N 

7 C<I,J> = 0 . 0 

8: DO 100 K-i ,M 

9- CC * A<I,K>*B(K ,J) 

10. Cd.J) - C(I,J) + CC 

11:100 CONTINUE 
12: RETURN 

13: END 


1: SUBROUTINE MVMULT (A,B,C,M,N> 

2C THIS SUBROUTINE CALCULATES THE PRODUCT OF MATRIX A AND 
3 : C VECTOR B AND STORES THE RESULT IN VECTOR C 
4: DIMENSION AIM, N), B(N) ,C<H> 

S: DO 100 1*1, H 

6: C(I)-0.0 

7: DO 100 J*i,N 

8. CC*A(I, 

9 C(I)*C<I)+CC 

10 i 100 CONTINUE 
11: RETURN 

12: END 


1: SUBROUTINE MTRANS1A, AT, M,N> 

2 : C THIS SUBROUTINE CALCULATES A TRANSPOSE AND STORES RESULT IN AT 
3 : C M * ROUS OF A , N - COLUMNS OF A 
4: DIMENSION A(M,N> ,AT <N,H) 

S. DO 100 J*1,M 

6> DO 100 I«1 ,N 

7. AT(I,J) ■ A( J,I > 

8:100 CONTINUE 

9. RETURN 

10: END 


259 


1: SUBROUTINE VADD(A,B ,C,M> 

?;C THIS SUBROUTINE CALCULATES THE SUM OF VECTORS A AND B AND 
3 • C STORES THE RESULT IN VECTOR C 
4. DIMENSION A(M) ,8(M),C(M> 

5: DO 100 >1 ,M 

6 C(I>*A(I)+B(I> 

?: 100 CONTINUE 

S RETURN 

9: END 


1: SUBROUTINE VTRANS(A , AT ,M) 

2 • C THIS SUBROUTINE CALCULATES THE TRANSPOSE OF VECTOR A AND 
3 . C STORES THE RESULT IN MATRIX AT 
4; DIMENSION A(M) t AT(l,H> 

S: DO 100 I»1,M 

6. AT ( 1 , 1)*A < I > 

7. 100 CONTINUE 

8. - RETURN 

9: END 


1= SUBROUTINE MC0PY(A1 ,B1 ,M,N> 

2-C THIS SUBROUTINE COPIES A INTO B 
3 = C (BOTH ARE M ROUS BY N COLUMNS) 
4: DIMENSION A1(M,N) ,B1(H,N> 

S< DO S00 I»1 , M 

6. DO 200 J»1,N 

7. B1 ( I,J)»A1( I , J) 

8- 200 CONTINUE 

9. S00 CONTINUE 

10. RETURN 

11. END 


1. SUBROUTINE MPRINT(A ,M, N.NCHAR ,TEXT> 

2. C THIS SUBROUTINE PRINTS A MATRIX WITH TEXT AS A HEADING 
3C M - ROUS OF A , N - COLUMNS QF A 

4. DIMENSION A<M,N),TEXT(S0> 

S< IREM»MDD(NCHAR,6) 

6. IF(IREM.EQ.O) NUORDS-NCHAR/6 

7> IFdREM.NE . 0) NW0RD5-NCHAR/6 ♦ 1 

8. UR ITE(f>, 10 > (TEXT (I > ,1=1 ,NU0RDS) 

9.10 FORMAT (//IX ,S0A6//> 

10. DO 100 X»i ,H 

11. WRXTE(6,20)(A(X,J>,J»1,N> 

12.100 CONTINUE 

13<20 F0RMAT(5X,11G11 .S> 

14. URXTE (6,31) 

15.30 FORMAT (//) 

16. RETURN 

17 . END 
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APPENDIX E 


OP poor 


page r$ 
Quality 


NUMERICAL OUTPUTS FOR THE SIMULATION OF ADAPTIVE CONTROL 
SHUTTLE HARD DOCKING TO FOUR-PANEL SPACE STATION 


DURING 


r> c 

9 

• • 


o e o o • ■ 

o o o o • • 


i ■ a ■ a m 
j» t 

> -* A. i A » ■ 

> A i) U A ■ A 


• *A • • 
XX • •« * • 
*» » 9 X *•’ • A- 
• • &. m -A. a 


; x • ••►< 


» «*A « 


X ««*ASA o 

*> 0 . a. a 


AS«AM«IO 
► MA SANK 


o *» * a • 

r A N A MM K 

^M«5MNO 


XAM AM 
A M «3AA • 
• » O « • *K 

AE A AM 9**' 
AWN 3A.WM 


% 


o o » 


AAt ^oe * 

EA • •« Ayw AA |A 

• " A ***i ’•*•*»*»% 


-A A NAS • • 

X U *«*MA» O o OX 


% ttO • *00 03 


O A A SAW • • »AO 
5 •» • • «AAA •« 
EX *AEAOOOAt 

3NNNA ►« A 


AA kAAIAAOPAM 
AAm3nA««# *AJ 

• • • • ANN ^E 

* EM AMMAAAM « 
AM •3AAOOQX A 

• • o • » ••• mm mat 

A « A X ANN • A» 

AAA AM4 • « *30 
A AH 5 A AM lltfO 

• • • *ooo 

• A •«« I A»««0 • 


IOOOAOC OMOOA 
- -»9 • A • *# ~ 

• • 


c eeoac ac-c e c 

OOOdOOAOAO 96 

eeoo**o *09066 

? 000 *0909900 
OOOAON 6 • O O O 


90 *000009006 
•0900000 •OOO 

90NCC0C0— OOC 

R 0990990*099 
•A • • • • • N • • • 


> OOA-OOO 009 O OO 
•C *000009006 
•09000009000 
9CAOOOOO *006 


261 


o 




e o • *ooo«oo coc oc«»o • 
O 09*90.00 *0000000 *oo 
O OKKQOO-*0000000« o 
• ooo ooo*moo o oooo»> oo 
• *«n «n • • *9 ■ ••••• • m 
• • • • 


>000 *or o *000000000 
>ooooo»ooooooo « *00 

» O 0000900000000000 
>0 0 0 009 0 <u 000 00 — — O O 


o *0000000 *0 
00000000000 
00000000000 


« OO *00003 30000000** 

00 • *0 O OO O O 9 O O O O -*0 * 
o o «***o o oo o o o o o o o *o« 

0099000 *0000000909 


o o • *000000090 *0000 
0099000 *o J O *090000 

O 00000090909000009 
* *0*1 * • *9 • • *9 *9 • • • • 

• 0 



>oooeeoceee«cc 

1004003000040 
- _i 009 OO OO 000900 
>9000 *0900000 *04 
MK030«00CCCCO*0« 


% O 9 O 00900000* 

• • • • • *o • • • • • 


>000900090* 
1 ••*••••* • 


>0 
> o 
»9 


• 0**000 *900 *900000 


262 


df.tr in.fi iK.it in. if ito.it lot.i* i.ftoi i»*m i.osio i . 1**1 i.tioo 

l.ltll .*♦(*? .10100 .(((Al*0t 


ORIGINAL PAGE 13 

OF POOR QUALITY 



lOOOOOOOOOOOOOOO 


•I till • I • • • I 0 

©c© o c o c © c © o o o o © c. © c> o o eoo c c c 

• I • • • • • ■ • 

oKoe^ottaK 

••©©©©©•••VtOOOvt^OOOOOOOOOOOOOOOtO 


ooo ooooooooooooooooo ooo ooo 

>'«ro»«nc»r«trt»i>»»ttiit««npiwo» r r. c • 
i • •••••••••••••*••••••••••••••• 




8 

m 


o 




© 

© 

o 

© 

c 



o 


o 

© 






© 


ooooooo © 
OC 0»<l0*' c 
m e© oooo •© 


8 

3 ©© 


. eeooMcoc 

0000*009 

ec o© oec c 

I O© © ©M O © © 
OOCOItOCC 


o oo *ooooo 

S oorc oo© t 

OO —*© OO o © 


to# O OO OOO O © 
W * 000000 © 
M ITIOOOOOOO 
J — © © o © © c c 


263 





o ooooaoooooooooooooooM 

S ooooooocc ce oooooooooo 
0900000000000 00 00000AI 
C OOOOOOOOOOOOOOQOOOOOM* 
O OOOOOOOOOOOOOOOOOOOOAI 


o eooooooooeoooeeooeoco 

O 00900 30 30000000 900000 

© ooooooo 00000000000*00 

O OOO OOOOOOOO 9000000000 

o ooooooocoooooooooo«*eo 



© ooooooooooooaooooaooo 
o oooooooooooooooooioooo 
© 0000000000000000*0000 
O 00900000000 90000*0000 


© 00000000009000M000900 

© 00000 0 90 000000^000000 
O OOOOOOOOOOOOOO *000000 




o-ooo oeoooooccc 

• • • • 

» mom mm— mm— mmm • * 

© «0V%O'M*»»O «*0 ©*••*•>» 

M O O « 

* • * ** • • • • 

OOOOOO OOOOOOOOOOOOOO OOO 90 OOO 

mmmm—mmmm o«« ooo 
ooooo* • ©©•* «r-*©^w»«*©o© 

ooeo ©ooooooooococooc ooocc Q o o 

M « *M O * © 

t o 

M««|m 

• • • •* ■••••••••••••«*•••*•• • •* • • 

OO 0009000000000000009900 OOOO 

• • •!••»•••••••••»•••••• I • • • • 

** **99909000 

**•>**•*««* m mmmmm mm 01 *mmm mmmmmmmmm 
OM»nn«iito**OMOM*i« —mm o«to mm — mm mm m « 
oooo m —m — m mmm m mm— o—mmoommmm— mom 
m — m mm — m — m — m—m *««•* rm — mm — — m 

••*•••****»••••••• • •• 

o o oo o eo o ooo ooo ooo ooo oo oooeoo 
■ •••••••••••••■•■•• a ** •••••• 

oo**-*-*o©©o*999*o*.©©9©ooo**»oo m 
m m mmm — mm mm ommm o —mm mmmommm—mrnmm 

* OM0ooo«*omo"%MOtf»»»ooMnoo ooooooo**- 

m—m —m— m—m— m—m mm —m ——m — — — •—rn—mmm 

• t • *•••••••• • • 


264 


D»nr 


ORIGINAL PAGE IS 

OF POOR QUALITY 



a o 
o o 
o o 
© o 





o o 
o e 
o o 
o c 


o oo 
o o o 
o o o 
o o o 



© © © 
o o o 
o o o 




o 

o 

o 


© 

© 

o 

o 

e 


090 
CO© 
© © © 
o o o 
© © © 


9 9 
O © 
O C 
© © 
o © 


© © o o o c 


© © 
O 9 
O © 
o © 
o © 




oooooo— © coo©oo©ocooooeooo©oeco 
09000000000000000000000000000 © 
O O o O o o «•© ©oooooooooooooooooooo* 

000000«*0000000 000900000000009-* 


oo 
© © 
o o 
e O 
o © 


o o o © o o 
o © © © © © 
o o © © o o 
o o o o © © 
©09000 

o o oo o o 


ioooooooooo 

>©©ooo©c©o© 

t©00«009000 
l C ©0 — 00©©©© 


oooocoocooooccocmcc 

©000000090000000900 

oooooooooooooooo«©o 

©oo oooooooooooooooo 
ooooooocoooooeooa*oo 


e c 

9 9 
o e 


x c a. c 


OOOOO OOOO o ooooooooooooo o«*ooo 

090000000000000000000000900© 
© 00009000000000000000000«-»000 

© m © o o c ©oooooooooocooc o © © o c © c 
o««ooooooooaoooooooooo90ot»ooo 


o © 
o o 
o c 
© 9 


o e o o e e 

9 0 9 0 0 0 

e e © o e c 

O OOOOO 


90000000000000000000000«*©00004 
ooooooooooooo© ©ooooooooooooo o< 
ooooooooooooooooooooooo»»ooc © o < 

90000000000© 90000000000000090C 
•> © O o o o © © o O O O O ( 


>coooeoooo«*oooooo 


o © 
o o 
o © 
o o 


265 


leooo 
1 o o o o 
» 9 o e o 
> © o e © 


o • 0 

M % x fc x a. * •• 


w-*p»o ra.c 4 J 


.X4Z4X4X4X4S 1X43K»- 

i«»«»^ooo o ooc oooxxa 

, N NK»-^nMnii]<i «ir* r» ► ui 


0900000 
o o O O O O o 
0900000 
O O 900 O O 


000000009000000000 
o ooooooooooooooo o o 
oooooooooooooooooo 
000000090000000000 


100000004 
► O OOO OOO I 
>0000000* 
> oooooo o « 


o — O mO « • »«Vn«*4X4 X4Z 4 X4 X ft X I 
•t«*4S 4C & r<*-*X4 4«*<*nM«««4in>>00( 
4 X MM MM )• MOO. 33aM^*>MMHM4»>N«*»« 


>09 9 00000 OO 000000090000 000900 OOO 0 9 
>009000000000000000000000000000000 
>oo 9 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOM* 


9 9 000000 

O no* OM3IHO 




266 


#1 UF» « a k(TO«»OI 


ORIGINAL PAGE IS 
OF POOR QUALITY 


i4o«oooF>oro»o«o«o« o « o 01 •» « o 
]^o»o«oivo«o«oho«o<«o«o<iimm« 
I^0«0»>09 okio-oaio«omo« omo — 
I040wo»»0-040-c«0»0il0»«l «• 
ioro»o<wooocooo#ooo«09«4» 


o o o o o o oooooo e o © 

4.4.44. • • ••»••• * • 4 

o«offoco«etofeMi «« ••oo*o»o#.c 
o *o»o«ieMOi"oro«>««»HU*>o«o «c*»e 
o«to«io«o*ome«*QO«*«»*»o^o«o«kO 
o 4o»*o»o»o*c;o*2* K i;2' # 2f c - o 
o o«ono» ooo« t * » «o»owoooc o 


r««in»»>o 6 oaooooooz*K 

4rmKf)4U 


• • I • • • 4 

•*n«4^»o^««in*>«oA<*«^a:Lz&xit 
► ►►►►►EKxwwxwa.a. 33?aMN*-iNL ** 


* * • • • • • • • • i » * 4 

Mno»omoMo«o>> o«i o«onooo 0 <vo 
Nooeo«o<>a«o->oi«ono4o«oM-M 
OMOM OKOOO4O0O40 «000«9n«« 
--O4COO»>O-O-*OM0«COO4OOI«I4 

OMOQO«oeoitoi^0 4 0i*i3n-<n 

• • I • I I 

<r a x ar a x a r a r asazazaxax 3 a « 
4 444 KOOQO O0OOOOXZU 


O 000400 oooooo 000 

• » • • O • O at-00009 * i 

4 0 000 0000040MOC040000* 04*00 
«*0»«0<V000 40 404 OF>C«F>4^0xOir>0«< 
no»0 4000WO^OK090n044 04 090irt 
m o 1400 0 4 0 OOOOOO *► * 4W40-040 o 
O O 0*0 4 04 O 4*0 4 On O —ifl 9*4040004 

• 0 

MMtraoOMono «a»arDi«9axiXft s o. 
a n ca.z<«««a ca.c ax ^-aaa--*«m4* 4 } 
» »»*»*as«iiii4 4>ai9a dw9»^nn*»n 


55 ^ 5 , 


,oio5o ioio 
ronooo*90o 
*040«00000 


>44 4 00 o O 


00000 

00000 


o 

tacos 


c • ♦ 

C * O 4 


»k4 4U 


00 00000 00000000 o o 

• • 0 0 0 0 • 

o 4000 00 004c00«40* 4n«*ftOOO«* 
OMO* O OOO 0 40<*00444 «K *NO««On 
0»>0*04000r040*e*>«44cn0»0^ 
0400040^04 Of»0^«44N44-040«I 

0 4»0« o 000 OO 00*00 OOOOOO 00004 
• • ••• ••• •••••« ••••••••••• 

•4 0 o«*ooooo<40<40 M ?**? v 2 M 2 e • 

00 0000000 0 

»><VI4VI #00 44M 4 4*0 4* ooo«a title 
aracac«sacata««aaa — 00 44* 
****>>*a»M KMMMaaaaaa^MMaoa 


267 


00000* 0 



o m #* o ai cm 

o o o o o © 

• • • • l« 

•*o«o»-»re 

— C 9 C O f% *%t 9 

►oro»- «« 

• c«*o^a»-o 

«cro *»« v — 

• • • • 

IX A. SO. SNA- 

ooficcxrx 

ft. MM >->>0X0 
MM3 



X ft. X ft. Xft-ftl 
► ► MM ►►s *1 


i*r>-(Oo»-€ koi 
ft X ft. X - —ft ft ft — - 
IMKMM AXOl 3X > 


.XftSftSftxftzt.rLCftSftca.3 mm 
'm«*« 4 *>v>mocooooocoo:cxk 

4 M ft M ft MM MM — ft* M Ol « VJtn^KB^O 
MMMftMftMMMMOXO 


• • I I I 

ocontoo MM 

0«04 0 ^ 4 -* 
O — © 9 ©M © M 


OO-*O»O«O«OOOl 0 OK 9 — -*K««O 0 ‘O<VO 4 OOOftOK-Omo^ 09 O KO 4 
4 O9OKO0lOir0«iO*.OO9 — — 09 OOOOOMOOOMO— o AIO — 09 0 — M Ol 

© 0904 04000^ 34044K 9 — 9 090MOMOM0009000M 0«U O •» M — 

OlOPOO 0«C 0000 — 0 — M 4 CM — OOOOOOOCUOOOOOMOlTO — 09 — — 

— 044MO- 0 9© — O 9 09 *4 <¥ — 9 O 09 0MQ909 0«90V»0M09 040 O 9 M 9 


X ft X ft X 3 ft V 

oo oo oxzu 

• 4 0 40 M «4C 


> lOM O e — M99999M9M «9ftXftXftXftXfl 
. c ft r — — ft c ft = ft c—— a X a. — — m m •« 4 4« 
-MMMftXMM VMM K ft. ft. 300 MM M MM M M ft • 


O O O OO O 
ft ft I ft ft ft 

9. OOO— O M «9 
KOOOnOO 09 
•0090909 910 
OOOO 40*00 


OO o o o o o oooooooo ooooooooo ooo 

♦ ft ft • ft • ft • ••♦••♦• ft ft • ft • ft ft ft • ft ft ft 
O9O90M00 00O— OM9 940 MMO40O«O4ro«a0OMO40CO0MO — ©MOO 
O «• O «0 * 34 ftKOftOMOIOft ft 90 00009 09OM00OMO9090M0M OO 
0M09C— OOO 0*0 9 009 MM O 9 ft » O C C ft ft ft O O OO 0«C ft O — O 9 O 9 ft ft 
ft — O OO 400 OMOMOM 90M (ft Mm 90 — 0009090 — OmOHJOMO — Ok. At — 
O — 0 4)000 — O — O— O O 9 M o ft x0004040w0ft000000c«00004> n 




19*0009 — — MOOmO— 9 ftft M O — ft.Xft.Xftl 

> xa Xft xwtft ft ft x ft«««* ft ft. ft — — nftftv « 


kXftX fttftX ft.M ft 
>OOOOOOOXX IT 
Ift 9940 OK 44 U 


269 



ItTD.OJ T"| 0.00000 


o o o o o o oooooo ooooooooo oooo 

• I • • • I • • ft. ft. ft ft • • • • ♦ 

n a o o «o«c*wo«or» o*> e*o-oioo-o«o»o«Of o«oro-of 9 

• 0rO«OOe«O«C^<Ute«C<6OOOO0.0CC^6«O«O«Cl9O#^K^ 
tojo«cooeo»orf'm--«»«o«o»'o^o»>ooo-o«o*ONOoc«)o 
oooomO«o«oho«cai«m«oi"o>>o9o»okoco«o«oco^«kk 
9 o«o«oro-o«o*<KKoin4e mooo«i09 0»o«oao«ooo4^4 « 


:4C&rxxiz&z&xariLX&r4 3^^ 
• •dwm «««*>n»-ooceaooo coz » * 

IM»-t -mwfll* «4rftinK«< N» 

OXX ft.**a.XMXX3r a 

•># « 3 


00000 o oooooo ooooooooo OOO 

• • I • I • •••♦•■• ft ft • • • t • • • I I I 

»«009 0 0 0«0«0 90 » 9 O«O 4 OKO»OKO<*OinO«O 0 -OOO 09-0 
i0o<no»o^o*> oro-^m «».«o<"0 9 o%ooto^o»«o«ONOxc «Aim 
>«notnooo«o^oino» ooo»onoooAio«o0 9 0i> o<«oo>3 o 

> -o-oreoc^o-oe^-^i-* on o n o « o 00 «o*>o» o« ovi o«»» — 

> 9 on o^ono« ooiO«m«io«n o«oook09 0<uoko«oo>o j 00.4 

»*oO« 3^9-*0 90«0« «n**<t«OfllOMO *009090 9 3M O •• 09 O •• •* 
• • • t • • • • • I I 

w • <r^ « c - 9 nv> 4«o«rn«ofcsax4i^x4ikxax&iftxo k x2a9 
ca.txi--u.ra.ii c — — a. x x .--»9 900 on x aaoooaoooox zu 

- »-XXxa,X.MMMMMMa.tX 3 3 => M X X NN»>NX ^N--M99 9 # 04 O 

XX MM MXXMMX MX X MM Z O 


OO O O O O O OOOQQOOO OOOOOOOOO OOO 

♦ • ft ft ft ft ft ftftftftftftftft ft ft t ft ft ft ft • ft ftftft 

o-ooo<o«no-090N-«ftcnft 90xooooo~one909090«cn c 9 

0ft090tfl090ft0no999ft40M 00900 0 lAO 9 0«9ftOft>0 09l»0-99 > 
oo enoooneoono«9xo9X9Mooxooeoooo9ono-*exeoc «ft9 
oo oxoooo 090<9o-*xoxnnoAiooo90 non o — o o o 90 -o « 09 nx 

09 0X0X00090 oonnxooox 90000 OOOOOO COO 009 o o on c — 


• ft ft ft • ft ft 


271 


T»»N$»nsr w moo»l rtvtcn 






O O 




o 


o 


o © o o o o 

• I • • • I 

•W O rsi « o 


9 r « « Q O 


• 0 0 



9 

*> 

0 



O 



e 



o 

o 








0 0 0 



eocc eoeceeoc 

• iii 00110001 

i --w^oowNO-t ••••»• 

i m • o ooooo'VM-*«nn«a« 

— — — — — r 9 — 

• i * • i • ii« • i i i • i 


O 9000001 


>000090 O 00009000000000 


_ _ _ 00 ? 00 c_ _ # 88 

> ©©©** — — « arr corr *•««• 







• i i i i 


eoo oooooooo 

i • • 00000000 

9 oocc*>«»i>«»e 
ooooooownt^w^nc*** 

• it • • t 










ooooooo * 5555 





• II I I III II I I 


>*cr»«»*oo-o»r 

u • •••••••• 

> 0 III III 8 888 0 

88 
w 

O O 

at 8 8 

•* >*»i.8n«AV«0»*l»>M0OMOM»M4 O X 

8 «noaao0«rr0rr-0»0>>0>aiiM l 

m • • • • t • • •«< 080 t|N •••••• 

I t00«0080K00*>00 • • • t« 

-0 0 0 080 0 00 

« 

8 

K 

k8An«0«*Ae^c r8fc0t«0? 

0 »‘— •••■ I«»»I 8 •»•>«* ~~ 0 * •••••••• 

8 I ** 0*000 0 00 



272 


ooooooooooooooo 
**•**«*•« o <w c 0 CO««I« 


• ■tea t 

ooooooooooooooo 
onowMCtfiM* •• • ties 
r* o — ■’» • • 


i « • ■ • 

— m ■••nnKtMi 

oe oooooooc 


• • • 

i •• m m 

o c o 


V. — W, 


^•rrt n«t ail 

it o« ». 


»rt • 
i in 
i a o t 
* — ■* 


o o oooooooc 

• • • • I I I I t I 


- r» « • o m 




• • • • • • 


oo cooooeo 

• • •••■lit 

««l- 0 tlM». 0 «i« 

3 * « r> *x ©■ •*> o •* o 

• ► «••••• o • 


Or a> 4 » 


«*5S&«S£*rt- 


*• ►.••—so 

r m 

* 

0 c »«•*«» 

«s ns — 

— 

t- •«■ -*-aO J 
• »•••• 

• • I 

r* #**« •*•* 
M IWO OlAr 

• 

d 0 O 9 MM 

«W If---- 

1 • • 


«»•• o r r •• o iw 

teMMiHfiniatt -aO 

• • ■ • • 0 0 


c «ott 

V ► 1— mm mm 

• I • I 


o o e c c c c o o 

• • ••••••• 

mm mm mi ^ O 0 C W. t* © 


COO O 

• • • ■ 

mnitoo 


• • • • i • 

O O O O 3 O O O 


OO o o 


w. a" 1 ? •••“•• i»1^ •• 

• • I • • • ■ 

« M^MMatn.aaa.MnMtiMn 

o ooooooooooooooo 

~~ ~~~itti •••••••• 

i 0 M t O t t 

* — « — m- mi — m — n auri 


ooooooooooooooooooooooooooo 

Ooooco«*ooooooooccocococoooo 

000000000000000000000300000 
cooooc««ooooooooocooocooooco 
OOOOOOaaOOOOOOOOOOOOOOOOOOO o 


m. ooocttoooooooooc ccooooocc coo 

X 003000030000000003 90000033 3 

« oooooooocooooooocoooccooooo 

O O O O O ** O O OOOOOOOO OOO 30 30000 3 O 
•» OOOOMOOCOC oe ococcoooococooc 


eoaoooooooooooooo^ooooooo«o 
K oooooooocooooooooooooooooo o 
X 003000000000000000 0000300—0 
cot OOOOCCOCOOOOOOOOOOOCCC.C o 
o O 03000000000000 OOOOOOO*- O 

• ••»••••••••••••♦••••••• •• • • 

o 

m mm 

•* o 

« • 

• ooooooooooooooooooooooonoo o 

< 003000000000000090000000000 
* OOOOOOOOOOOOOOO oooooooo^ooo 
• 00000003000 300003300030000 
a — O C ooooooooooooococoooo^ooc 


273 


tOMO»«t* 00000* eoooo* ooooo* 

00000* ooooo* ooooo* ooooo* oooQo* ooooo* ooooo* OOOOO* ooono* ooooo* ooooo* 

ooooo* ? 0 «I 0 |?I* ooooo* ooooo* 


ooooooooooooooo 

000000000000900 







OOOOOOOOO 

00003000"* 

OOOOOOOOO 

00000030* 

OOCOOOOCAi 


oooooooo o 
0000000—0 
OOOOOOO** 9 
OOOOOOOOO 
0000000^*0 


000000*^00 
000000900 
OOOOOOOOO 
000000—00 
oooooo • o o 


O O O O 3 0 3 

00000 — 00 o 

OOOOOOOOO 
00000—000 
OOOOO • O O O 


000000000000900 
00000900000 oooe 

OOOOO *000000000 
OO OOOOOOOOO oooo 


000090000 

000090000 

OOOOftJOOOO 

OOOOOOOOO 

oooo • o o o o 


o 

o 

o 


000-00000 
OOOOOOOOO 
000*100900 
oootnooooo 
OOO *00000 


00000090000000000000000000000*% 
00000090000000000000CC ooocoocc 
— OOOOOO *003000000900000 0000009*9 
*9 OOOOOO— 000000000000000 0000000*9 
* 000000*90090000000000000000000*9 


M OOOOOOOOO 
r OOK 303303 
9 OOOOOOOOO 
a OOO 900900 
WOO *009000 


_1 

a oo o o «o oe e 

09 0909000 
9.0000 •OOO'* 



0-0000000 
OOOO 300 0 0 
OOOOOOOOO 
O *0990900 


OO *090009000000 -9 990009900*90000 

0090000000000000000000 oooooeeo 

90 — 3 99000 90009000 9 OO 9 9 3 3 3*3333 


000009000 
OOOOOOOOO 
*% 00090000 
>00000000 


I 9 900000000 9900090 3 90 OO 0^003 990 
I *00000000000000000000000000 0 CO 
< O 90 09000000000 9 90 0 000 9—009009 
1-0000000000000000000000 *000000 


274 


CO o o o c 

o o o o o o 

© o o o o o 

o o o o o o 

o c co c c 


o o 
o o 
o c 


e c 
o o 
e o 
o o 
o c 


c o 
o o 
o o 
o o 
c • 


o o 


o o o o 

o o o o 

o o o o o o 

CO CO © o 


CO oo oc cc 

oo oo oo oc 

oo oo oo oc 

OO OO OO OO 

oc oo oc oc 


o o o o o o 

o c o o o o 

o o oo o o 

o o o o © o 

o o o o o o 


o o 
o o 
o o 
o o 


o o 
o o 
o o 
o o 


oo 
oo 
o o 
c o 
o o 


oo 
o o 
oo 
oo 


o o 
o o 
o o 
o o 


o o 
o c 
o o 


o o o o o c 

o © c c CO 

o o o o o o 

CO o c c c 

O O O O o o 


o cooocoo* 
O ooooooo* 

c C OOOCOOO 

oo o o o o o o — 
c c o c o c o c * r - 


O O O O o o 

o o o o o o 

o o o c o c 

o o o o o o 

c c © c o c 


c c 
o c 
o o 
o o 
c c 


o e 
o o 
o e 
c o 
o c 


o o 
o o 
o o 
o o 
c o 


o o 
o o 
o o 
o o 
o o 


o o 
o o 
o o 
o c 


o © o © o c 

o o o o o c 

o o CO CO 

o o o o o o 

O © CO o c 


o oooooooo 
0000000*0 
e oooooocro 
ccoocoovo 
O OOOOOO — 9 


000000000000000000000000000000000000000000000 

ooooooococoooooocooocooocooooooooocoooocccccc 

000000000000000000 OOOOOOOOOOOOOOC© DCOOOC 0 0 0 3 0 
OCC OOOOOCOOOOOCCOOOCOCOOOOOOOOOOCCCOGOOOOC coc 
OOOOOOOOOO -OOOOOOOOOO OOOOOOOOOOOOOOOOOOOOOO© • 


C 00000*00 

o oooo©*-©o 

C 00000*00 
000000*00 
C 00000*00 


oooccoccocooeooocooccooeoooooooooocooeocoocoo 
ooeoooooooooooooooooooooooooooooeoooooooooooo 
eocooooooooooooooooooooocoooooooooeooocooocco 
oooooooooooooooooooooooooooooooooooooooooooo© 
oooceoo - o c o c o o o o o o o o o oo o ooo ooo o o o o o o o o c c c -coo 


0 00009000 
C 00009 OO© 

eoooo*-coo 
c ooooc ooo 
C 0000*000 


009000000000000000000000000000000000000090003 
oooooocoooooeoocoooooooooooooooooeeocoooooooo 
c c e oo ocooooococ oceoooooocooooooooc ooooc cc c coc 
o o.o ooooocooooooooeooooooocoooooocoooecococooo 
oooo •oooooooeooooooooooooooooooeeoooec - 00000 © 


o ©CO— ooc*c 
o 00000000 
c cooooooo 
o oootftoooO 
o ooo*»ccoo 


coocec c oooooooooooooooooocoooooooo 00 eoooocooc 

000000030000000000000000000009000030900000000 

eoocooooooooooeeoooooooooooooooccoocooooooeoc 

000000000000003000000000000000000000000000000 

o *000000000000000000000000000000000 *oooooccoo 


cccrococC 

OOOtAOOOOO 

o 00—00000 

OOO *00009 


000000000000000000000000000900000000000000000 
c c cc oc 000000 c oococooooococoooooocc 00 c oc c oc ccc 

000000000300000000000000000000000000000003000 

© ooooc cooooooooeooooooooooooooocoeooooooooooc 
000000000 00000000000000000000000 • o 00000000000 


ccoooocco 
O 33000000 
o c *000000 
o oooooooo 
00 • o c o o c* c 
• •—••••*• 


0000000 c ooooooocoooocoocooooooooccooooccooooo 
000000000000000009000000009000000009090000399 
OOOOC OCOOOOOC oocoooooooooeoooooooc OC coccoococ 

OOOOO O 90 O 90000 O 300000000000000000009000 3 90 300 
OOC OC OOOOOOOCOOCOOC OOOOOOOOOO • o o c cc e cooc ooco c 


0*0000000 
C *>0000009 
O 00000099 

c — oooccoc 

O •0000090 


00000000000000000000000009000000C009C0000090C 
OOOOO 90000000 OOOOOOOOOOOOOOOOOOOOOOO0OOOOOOOO 

occ ooooo oooocooocoooooooocoooooooc oc oc oooooco 
00000000000009000000000000 *9900 90009990990000 


— 00000009 

« oooooooo 
* 09000000 

* oooooooo 
• oooococo 



275 


ooooo 


C C O C CO 

O O O O o o 

o o c. o o o 

o o o o o o 

O © CO oc 


o o o o 


o — o o 

CM o o 

o • o o 

o — o o 

0*1 o c 


o o 
o o 
o o 
o o 
c o 


o o 
o o 
o o 
o o 
o o 


e o 
o o 
O o 
m o 
• o 


o o 


o o o o 

o o o o 

o o o o 

o o o o o o 

O O CO CO 


OOCOCCCC c 
OOOOOSOOC 
OOOOOO 3 O C 
ooooooooo 
ecccocccc 


o O O O 

CO CO 

O O O O 

O O O o 

o o o o 


o o 
o o 
o o 
o o 


o •» 
o • 
o o 
o — 


O O O C CO 

O O O O O o 

o o o O O O 

o o o o o o 


wo oo oo oo oo oo 
« o oo CO oo OO CO 
WO OO OO OO OO oo 

oo oo oo oo oo oo 
•o oo oo oo oo oo 


ooooooo oo 
oooocoocc 
ooooooooo 
ooccccoc c 
ooooooo oo 


3 3 


o o o o 

O O 3 0 

CO o c o o 

o o o o o o 

c C CO CO 


o o 
o o 
o o 
o o 
c c 


o o 
o m 
o o 
o — 


o o 
o o 
o o 


O o 
C c 
o o 
c o 


o o o o 

3 0 • O 

o O o 

3 0 O O 

CO • c 


o o 
o o 
o o 
o o 
c o 


O O 

o o 
c c 


o o 
o o 
o o 
o o 
e c 


o o o o o o 
O O 3 0 o O 
o o o o o o 
o o o o o o 
o o o o c c 


oooooooec 
©ooooooc. t: 

oococcocc. 
OOOOOO 33 2 
ooococccc 


O 000000000 * 0000000000-00000000000000000 ooooo— 
O OOOOOOOOO — OOOOOOOOOOmoOOOOOOOCCOC C 3000 OOOOOO 
0 00000000000009000090—00000000000000000000009 

o ccooococ o— ooooooccoo«ooocooooocoooocoooooo © • 

O 030000000—0 OOOOOOOOO *0300000000000000000000— 


ooooooooc 
ooooooo c c 
ooooococc 
oocococcc 
OOOOOO 3 3 -2 


OOOOOOC*>OOOOOOOOOOV>OOOOOOOOOCOOCOOOOOOOOCptOOO 
O OOOOOOO 3000000000*00000000000000000000000000 

0C.CO000C000OO0O00OOOOCC3000000O00O0000C0O— ooo 

OOOOOOO- — OOOO OOOOOOO OOOOOOOOOOOOOOO OOO 0000*000 
C OOOOOO — OOOOOOOCOO *0000000000000000003 ooon oo o 


ooooooo cc 
ooooooooo 
ooooooooc 
ooooooooc 
ooooooooc 


00333003300300*0000030030000300300039 - 003000 
C00«0000000000*0000000000000000000000»000000 
OC 0*0000000000 *0000000000000000000000 — 000000 
OOC* OCOOOOOOOOO OO COQOOOOOOOOOOCOOOOOOOOOOOOO 
o o c. — ooooooooooooooeoooooooooooooooooo»>oooooo 


OOOOO 3393 
OOOOCOOCC 

ooooecc cc 
oooocoocc 
ooooooooc 


■ooooccooooooooooocooooocoooooooo 

* 90000 0000001033000 03 9099000 300 3 0 0 
ftOOOOOOC OOO *00 0000 OQOOOOOOOOOOOOC 

roooooooooooooooooooooooooooooooo 
•cooooooooo— coooooooooooooooooooo 


0 0 3 0 0 1 


> © oo 
» o o o 
» e oo 
» oo o 
» o oo 


coococooc 

0090000 3 : 

ooccccocc 

000033933 

oooooooec 


© OC OOOOOO** OOOOOOOOOOOOOOOOOOCCC COOOOOCC OOO oco 
O 00909000 *0000000900000000000000— OOOOOOOOOOWO 
OOOOOOOOO — COOOOOOOOOOCOOOOOCOOOO — oooooooooowo 
OOOOOOOOOWOOOOOOOOOOOOOOOOOOOOOO — 000000000090 


000000009 

ooococccc 
ooooooooc 
oooooooec 
OO o O 9 0 OOC 


c coocowecooc ooococooooeoooooeooooooooooooocoo 

000393900009003000030090300 3 OK 309039000 900009 
o ooooo • o o o c OOOCOCOOOOOOOOOOCOOOOOCOOOOC o o o C 9 
O © © Q 9 O W 30000000000000000000001*000000000000 OOO 

e ooooowcoooooeooooooooooooooe — oooocc ooo o — c coo 


ccoccccoc 

000000099 

cooococ c c 
ooooooooc 

OOOOOOO OC 


0009 90000 

oooooooec 

00009000 O 

ooocoeocc 

OOOO 300 3 V 


9 0009000000 3 OO 3 OOOOOO 9000000000000**9090 990900 
• OOOOOOOOO OOOOO 000000000000000000090000000000 
9003990030303000000300303000300300— 0000000003 

— oooooooo ooooOo oooooooo — oocoeocooo • O O O O OOC ooo 


oooooooec 

OOOOOOOOC 

% oooooooc c 

O 900093003 

oeoecoc ec 


276 


UObOU* OOOUO 


original Pmi i§ 
Of POOR OuAtfTT 


ooooooooooooooo 


oocc cccroe»n»* » © »> *>©©.*. 
ooooo*’*««««©m©©©©«®*.^i 
ooooocinno««ngio«in« « c •« 
0000060>*rv«iPin»oain0» — m 
oecoc.*- — — — *w — — — • « - 


• it »i 


ooooooooooooooo 

ooooo7ocoo*«irrnioAjir«c 
0C00C^accci0-««0N9 m ►- — 
ooooo©*iH'o«©>’i*'9Mm« « e m 
C C COOOO-\-ipHf.l I r © © — Q 
DOOOO<*»x>*kn« m >9 © m 


I • I 


• • 


• t 


cc&cooooc o ooc oc 
o cc ooc»-®r.#'tc».©iro®»<-».c©*. 

O CS O © Q « ® © — tf>©*.M©M».M»r.»| 
C CCOOO* 9 /' 9 *^OKOCt»K* 9 IX 
c oooooo-C'M — • « o ivc f- o — t* © m 
c ccccc««c)xl^nn^rlnn«.«•r l 


• • • a • a a « 


OOOOOOOOOOOOOOO 

• •••••••aaa« III 

o oooooer>M»*> H ip|toN«ow«« 
o ooccori v'BxiMxaricxnMr 
O COOC 3 k« 

c ccocoir»ir>cB«ir**Mvii>> m • *» 
w COOOOkMV«»mm««9»««AI9 


• III II 1 



ooooooooooooooo 

O 9 9099«*o'*oomm«i© 0*0 «••# 
c cccccKKffK&^caimr^iMM^a 
C ccccca©r»»«r.*i« < i*©©©«-**® 
c ccocconr***offO-90-«K 

O OOCCCWKr.-^-l'if 

I I •••••••• 

ccccccceccccocc 
• ••••••••»••••• 

c ccc e.cc»9«9i.*^CKK»^ro 
o ooooo*^«*^o»Moo»i*«iirMn 

O OOOOC— ^—•©mMm*!© — © M V* m 
o oo9QOM»fvn«^*»i9A*nii>M» 
O OC OOC-««VIMW» --A.l'l ««<n 

• • I • > • • t I 



CCCCOC COOCOOODO 

c r © c c c ». © m m m © © m e n mo o nr 
o oooooJto «Mn»^^i«M>iv>*n^in 

C OOCCC*tB.«»*Otft©B‘Ml’)©«^©©M 

o oooo^»»r-*oa»«»tfliON--i»« 
O OC 0009»*“ m — m © © « *v © • m © 

• • • •••••••••• 

ooooooooooooooo 

O OOOOOrfMft ->m« 0 /t^O«*>-OOON 

e oocccm»«»M>«m»o«9«ii^« 
o ooooo**®M*»i*©*</>o*>«’t©©OM 
c ccccccoK-rtr»^ar o»«*9 

o OOOOO-* — — — r.*.©M:*MMM©©M 

» ■ • • ••••II 

ooooooooooooooo 

O OO O O OO© 

9 0000904)««l»/l4«f><*M»«Kn 

o ©©©©©**© ov*©*.«»r»©^o©©©M 

O 0 ? 00 d 0 pfW»«Mtfl 9 o«JI 09 ^« 


oo oo oo oo oo 


» o © o o o © o 


MM © * * *- ► e 

>o oo r> r* mc 


-• oo — — r<*> •< 

© • ©«* « © r o © « 

OO «»M flip — O OH 


oo oo oo oo oo 

• • » o — o m© m© 

• • • • • • • 

oo oo oo oc oo 

oirt mo M— mm oo 


oo oo oo oo oo oo 

©O «o ©m na. m« * » 

mm « r» «. • ©* © oo «« 

n« o vo Om v© ©i© 

OM Mm mm ►> M M*» IUN 

• •• I 1 I I • I I* 

oo oo oo oo oo oc 

v o * ip om c«t ©i « «io 


• I • •• • I • II II • I 

oooooooooooooooooooooooeooooooooo 

OlTMM©M©OM|Am»l<V®MM©«fy«M©©Min©«lMOMMOM 

»M©»MorvMrt©*MO*i»MOMV»i*»MOi»i«MMOM«»M* 

B B III • • VI B • B B B 

OOOOOOOC OOOOOOCOOOOOOOOOOOOOOC ooo 

© *• MO*MMMlAOOM®MMM|fl® © «V © m»V® OMOOmmO mm 

© M © © m©^mm©mOMMV*m»»*I« mm«*MM®Om© M © © M© 

• t B B B B BBB B B B B I* 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


• O • « OB 


B B 

© o 

MO 


© M 
• M 

on 


M * <*. 

© C © < 

• B 0 I 

IPOH'I 

* © © V 
r V 61 

O ©MH 


joocoocooccccc-coccccc c c cc 

I B • B • IBIBlBBBBIIIItllllll 
►m»M© C©r.«fU — ©••**©•.*.*.©©*,«« © 

• MB* •»» M M .* MM ©» B» MM I* ^1 M MO MV» ©» © ®M© 
.MM©.© m«® ©m© oo © © o © «i>nBvvv« 

• mmo aoo«*tv j* mo©m©oo«oo 

r©.®«*OMM^©|MM©l©\ M*. © ©©OMMMMMM 


BBB BBB B BBB BB B BBB 

>0000000000000000000000000000000 


B B B B BBB 


OOCCOOC cooooooooooooc ooc ocooocooc 


r««o4v©*iMO»»ov*o*M 

Umm9mmJWOT|*|4mI*«®* 

BB BBBBB B 


IOm©©o©«*©©©MM©k**V 

• m © II WmV V f mpPIBIBmmIBmm 

B B B BBBBB 


000000009000090000000000000900000 
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 
BBBmBOmV »VBVM©OnB®K|.tVB BOB^OOBBBMB 

► «*Omi*Omo©iO©««bb»«*m»i©©©>m© o*© © © omoon o 

M© m m© m©©W«©«©I©B*©OOCm|B©©M©I©©m© mm© m 
«©BBBB«B ©©©V*M©.rBB*MMMM«MO©©OBI*© * Al © * 
• ••••••••••a* ••*••••••••••••••••• 

BBB BBB BBBBB# BBBBBB BB BBB 


>©B*09 I^mOBIOv^mOOBIOVIIMO© «BmDB 
IBVBH^KniP BKBVBVPBtNBBiP M-N 1 
► ©©»©•*»©• m©©©«®m©*mo®©m©iv©* 
>a© a >dB«nV«VBMiA» JO©BI*«*lA«inM 


m • ©» b» © o © 
\ mb. ©or.©© 
I #. © © © © • BT. 
* * MO © MO© 


B B BBB 


c-y 


277 


HATNII AH 


OOCCCOOOC3 

OOOOOOOOO-* 

0000000009 

© 00000000-0 

ocooooooo^ 


O a M 

o o 

• I 

«*» c 

-1 9 


o o o 

ft ft • 

000090000040 
000000000^9 9 
0 0 9 OOOOOONi- 
000060000«M9 
ooooooooo-’* — — 


o o c 

• ■ 


OOOOOO 000***9<9 
OOOOOOOOOOM** 
OOO 000000^*109 
00000000 3 9 * >r 

eooooooeco9»> 


O QOOOOOOOOOOOOOOOOOOC 000000000900 OOO 
0000000000003000000000 9 00000009-0030** 
000000000000000000000000000000900009 
000000000000000000000000000000«00009 
OCOOOOOOCOOOOOOOO *0000000 00 00040000% 


OOO 

• • • 

OOOOOOOOO./'./'O 
OOOOOO COO — — 9 
OOOOOOOOO o V *U 

OOOOOO ooo^niu 
OOOOOOOOO-*** 


OOO 

I I • 

0000000009 4 *n 
OOOOOOOJOV'V^. 
OOOOOOOOO'U «- 

ooooeoooo<9-"**« 

O OOOCOOOO«V V* 


• * 


003000000000000000000003000090000000 
000000000000000000000000000090000 — oo 

ooooooooooooooc oooooo ooooooo*»ocoo*>-oo 
000000000000000000000000000090000«*«00 
OOOOOOOOOOOOOOO -00000000000090000<r00 


• ft < 

000300390— 
00000000099* 
OOOOOOOOO®'*'. 
00000000009 c 
00000000099*- 


00000000000000000000000000^00009*0000 
000000900000000000000030303030000003 
0000000000000000000000000000000**0000 
O 000000300000 330000<3000030.ft3000-*0000 
O O O O O O O O O OO O O *0000000000009*0000 *0000 


c c < 

• • I 

)OCOCOOC9M( 
► 03000 39 ©*** 
10000000 *** — 9 

I oooooo©——" 


O 0.0 0 0 00 000 0 00 0 9000000000**- 000 0*9 oooooo 
oooooooocoooooooocoocoooooooo^coocoo 
00000000000900000000000 o s* oooooooooo o 
o oooooooooooooo oooooooo o<«aooo«*o OO OOO 
O OOOOOOOOOO *000009000000 *0000 *000009 


OOO 

• ft * 


• * • 

OOOOOOOOO******* 
ecGOocccca ®*v 
000000000009 
OOOOOOOOO—OO 
O 00000000*0^*9 


09000000000000 OO OOOOOOOO 0009000000 OO 

ooooooooooooeooooooooo«*ioooo«woooooooo 

O 00000000000000090 a O ◦ 0« 3000*000900900 
OOOOOOOOO *000000000900 • O O O O *00000090 


OOO 

« • • 

OOOOOOCC c « st c 
OOO 3000909 /» — 

oooooecoe**n 
oooooooo afl**-* 

000000000*90 

• • • 


OOOOOOOOOOOOOOOOOOOOft-OOOOOOOOOOOOOOO 
0000003090000009000090000M0000009000 
00000000000000000000 — 00009*0000000000 
OOOOOOO *090000000000 *0300 *0000009000 


>00090000— **• — 
*00000093993 
tCOQOOOOOVIt*- 
»9000 9 090it**N 


0000000<30000000000«90000*»000000009000 

oooooooooooooooooo«*oooooooooooooooeo 

000300300000003000—30009300000000000 
OOOOO *000090000000 • O O O O *000000000000 


OOOOOOOOO** 

0000009000 

: ooooooooo® 

I OOOOOOOOO*" 
OOOOOO oc o«* 


OO 
ft ft 
19 — 
« 


278 


ISfckT.Ofc .«lkV'OJ ..JoO»J-Ok •.••1*1*0? •.102«1*0* •.moo-01 

eiJT.nk ..IklkO.O* ..MOJJ.O* ,|*90t«0» . *0010.0* ,10?*«*0* ,l»«fc«*0, 


ORIGINAL PAGE IS 
OF POOR QUALITY 


fs €* 9 * 
• o • 
- <0 « - 


o o o o 

III! 
Hi «9 » « 


o O O O 


• • • • 
• a 


o o o o 

till 


o o o o 
• •ft 
o 

f* r» <n ;/! 


• • • t t 


O 09000 




o o o o o o 
•> M 9 


O - 



9 9 0 9 9 9 


^9 JltfttM 

• t «l« o« 
9 O «• 


O O O O o O 


9 O ^ 


ft • 



OOOOOOO 909000000 OOOOOO 

■ §•«•••••••••••§•••«•• 

• m — « « ». xr. O 9 <W *» • 9 « O* l>« 4M • 

9900 4 IA 9 A*t^nt««(y« 4 xnO'M<*»> 

4 t - • 9 «*►•-• ••• «• MO X — • 9 0 -* — 


» I » • ••• I t • •• 


0000000000000000900000 

— » « -. * 9 9 > 3« 

^)»»» — o«oaxx^'Vi'xo<yo«%« 90 <^ — • 
9 Mf>o»mo«in«* 09 K 9 m««««n»kO 


• ••at t tat tt t 


oooocccoocoooocoococoo 

• t • ••••••• • t t I • ■ ■ • •• •• 

t»-4irMir*ir^t» ti»xKo«t tt« 
99 009-IO#Nlfl-<yo^- l n.'\(9 VIAA«4 
ntAlM»«C9«in o *» • 9 •> 

'*•«» 4tA«N3AN»lNMai««0«M«»« 

t tttt • • a • t • a a a 

OO OOOOOOOO 900000000000 

•o r» m •- m « «- *»» m <m«a m — © — «r««oc 
49AI/t4Kn.'U«ia4MOC94f«^^**VVI 
omcn««*o«4«niMoomfl«*4C on«<n 

a a a a a •• a a • ■ a ••• 

oocooooc oc eooooooooooo 
«vimr ownciwaMno«<w««ao«a«** 

9Ain^ani^a</irt« a*»nn 

ao«aa^<*i*tfto^4«aaAiiAoaoaa 

•* — IT — — f*©« a -> *U — Cir««« x *• — a» -• 

a • a a a a a • a a a a a a 

0009999990090000900090 

»>ooa«Qoa«{p<*«aoMAiaa^coa 

— — XX-«~0— X<X—M — 

• a • • a a • •• • a • a a 

OOOOOOOOO ooooooc OOOOOO 
0»^44«0^a0440 

va « a^ia oa-« *aa 

x M**#»X*«M09»»4n~9X«»9iVAt9X»»»»M 
u MX ««k «*m.n o » Aia 

M • ••••••••«•••••••••••• 

m • • •• a a a a a a • a • • 

at 

•- 

x 

K 


279 


BIBLIOGRAPHY 


1. Wang, S.J. and J.M. Cameron, "Dynamics and control of a large 
space antenna," Journal of Guidance, Control and Dynamics, Vol.7, 
No. 1, Jan. -Feb. 1984, pp. 69-76. 

2. Cameron, J.M., M. Hamidi, Y.H. Lin, and S.J. Wang, "Distributed 
control of large space antennas," paper presented at the 
Workshop on Applications of Distributed System Theory to the 
Control of Large Space Structures, July 14-16, 1982, Pasadena, 
Ca. , pp. 225-247. 

3. Landau, I.D., Adaptive Control: The Model Reference Approach, 
Marcel Dekker, Inc., New York, 1979. 

4. Kalman, R.E., "Design of self-optimizing control," Trans. ASME, 
Vol. 80, No. 2, 1958, pp. 468-478. 

5. Feldbaura, A. A., "Dual control theory,” collected in Optimal and 
Self-Optimizing, R. Oldenburger, Ed., MIT Press, Cambridge, 
Mass., 1966, pp. 458-496. 

6. Iserman, R. , U. Baur, W. Bamberger, P. Kneppo, and H. Sibert, 

"Comparison of six on-line identification and parameter 

estimation methods,” Automatics, Vol. 10, 1974, pp. 81-103. 

7. Ljung, L., "Analysis of recursive stochastic algorithms,” IEEE 
Trans, on Automatic Control, Vol. AC-22, 1977, pp. 551-575. 

8. Ljung, L., "On positive real transfer functions and the 
convergence of some recursive schemes,” IEEE Trans, on Automatic 
Control, Vol. AC-22, 1977, pp. 539-551. 


280 


9. Landau, I.D., "Elimination of the real positivity condition in 
the design of parallel MRAS," IEEE Trans, on Automatic Control, 
Vol. AC-23, No. 6, Dec. 1978, pp. 1015-1020. 

10. Martin-Sanchez , J.M., “A new solution to adaptive control,” 

Proc. IEEE, Vol. 64, 1976, pp. 1209-1218. 

11. Johnstone, R.M., D.G. Fisher, and S.L. Shah, "Hyperstable 
adaptive control — an indirect, input-output approach with 
explicit model identification," Joint Automatic Control 
Conference, 1979, pp. 487-494. 

12. Narendra, K.S., and L.S. Valavani, “Direct and indirect adaptive 
control," Automatica, Vol. 15, No. 6, Nov. 1979, pp. 653-654. 

13. Goodwin, G.C. , P.J. Ramadge and P.E. Caines, "Discrete time 
multivariable adaptive control," IEEE Trans, on Automatic 
Control, Vol. AC-25, No. 3, June 1980, pp. 449-456. 

14. Goodwin, G.C., P.J. Ramadge, and P.E. Caines, "Discrete time 
stochastic adaptive control," SIAM Journal of Control and 
Optimization, Vol. 19, 1981, pp. 282-285. 

15. Whitaker, H.P., P.V. Osburn, and A. Keezer, "New developments in 
the design of adaptive control systems," Institute of 
Aeronautical Sciences, Paper 61-39, 1961. 

16. Donalson, D.D., and C.T. Leondes, "A model referenced parameter 
tracking technique for adaptive control systems," IEEE Trans. 
Appl. Ind., Sept. 1963, pp. 241-262. 

17. WInsor, C.A. , "Model reference adaptive design,” NASA-CR-98453 , 
Nov. 1968. 


281 


18 


. Butchart, R.L. and B. Shackcloth, "Synthesis of model reference 
adaptive control systems by Lyapunov's second method," Proc. 
1965 IFAC Symp. Adaptive Control (Teddington, England), ISA, 
1966, pp. 145-152. 

19. Parks, P.C. "Lyapunov redesign of model reference adaptive 
control systems," IEEE Trans, on Automatic Control., Vol. AC-11, 
July 1966, pp. 362-367. 

20. Phillipson, P.H., "Design methods for model reference adaptive 
systems," Proc. Inst. Mechanical Engineers, Vol. 183, Part I, 
-1968-1969, pp. 695-706. 

21. Gilbart, J.W. , R.V. Monopoli and C.F. Price, "Improved 
convergence and increased flexibility in the design of model 
reference adaptive control systems," Proc. 9th IEEE Symp. 
Adaptive Processes Design and Control, Dec. 1970, pp. IV. 3.1- 
3.10. 

22. Landau, I.D. and B. Courtiol, "Adaptive model following systems 
for flight control and simulation,” AIAA 10th Aerospace Sciences 
Meeting, No. 72-95, January 1972. 

23. Popov, V.M. , "The solution of a new stability problem for 
controlled systems," Automation and Remote Control, 24, 1963, 
pp. 1-23. 

24. Landau, I.D.,. "Synthesis of discrete model reference adaptive 
control," IEEE Trans, on Automatic Control, Vol. AC-16, 1971, 
pp. 507-508. 


282 


25. Landau, I.D. and H.M. Silveira, "A stability theorem with 

applications to adaptive control system," Proc. 3rd IFAC Symp. 
on Sensitivity, Adaptivity, and Optimality, Ischia, Italy, June 
1973, pp. 282-286. 

26. Bethoux, B. and B. Courtiol, "A hyperstable discrete model 

reference adaptive control system," Proc. 3rd IFAC Symp. on 
Sensitivity, Adaptivity, and Optimality, Ischia, Italy, June 
1973, pp. 287-289. 

27. Anderson, B.D.O., "A simplified view of hyperstability," IEEE 
Trans, on Automatic Control, Vol. AC-13, 1968, pp. 292-294. 

28. Monopoli, R.V., "Model reference adaptive control with an 

augmented error signal," IEEE Trans, on Automatic Control, Vol. 
AC-19, Oct. 1974, pp. 474-485. 

29. Morse, A.S., A. Feuer, and B.R. Barmish, "An unstable dynamical 
system associated with model reference adaptive control,” IEEE 
Trans, on Automatic Control, Vol. AC-23, June 1978, pp. 499-500. 

30. Narendra, K.S. and L.S. Valavani, "Stable adaptive controller 

design - direct control," IEEE Trans, on Automatic Control, Vol. 
AC-23, Aug. 1978, pp. 570-583. 

31. Narendra, K.S., Y.H. Lin and L.S. Valavani, "Stable adaptive 
controller design - part II: proof of stability,” IEEE Trans, on 
Automatic Control, Vol. AC-25, No. 3, June 1980, pp. 440-448. 

32. Morse, S. and A. Feuer, "Adaptive control of single-input 
single-output linear systems,” IEEE Trans, on Automatic Control, 
Vol. AC-23, August 1978, pp. 557-570. 


283 


33. Morse, A.S., "Global stability of parameter-adaptive control 
systems," IEEE Trans, on Automatic Control, Vol. AC-25, No. 3, 
June 1980, pp. 433-439. 

34. lonescu, T. and R.V. Monopoli, "Discrete model reference 
adaptive control with an augmented error signal," Automatics, 
Vol. 13, 1977, pp. 507-517. 

35. Narendra, K.S. and Y.H. Lin, "Stable discrete adaptive control," 
IEEE Trans, on Automatic Control, Vol. AC-25, No. 3, June 1980, 
pp. 456-461. 

36. Suzuki, T. and S. Takashima, "A hyperstable scheme for discrete 
model reference adaptive control system,” International Journal 
of Control, Vol. 28, 1978, pp. 245-252. 

37. Johnstone, R.M., S.L. Shah, and D.G. Fisher, "An extension of 
hyperstable adaptive control to nonminimum phase systems,” 
International Journal of Control, Vol. 31, 1980, pp. 539-545. 

38. Narendra, K.S. and L.S. Valavani, "Direct and Indirect model 
reference adaptive control," Automatics, Vol. 15, 1979, pp. 653- 
664. 

39. Narendra, K.S. and L.S. Valavani, "A comparison of Lyapunov and 
hyperstability approaches to adaptive control of continuous 
systems," IEEE Trans, on Automatic Control, Vol. AC-25, No. 2, 
April 1980, pp. 243-247. 

40. Monopoli, R.V. and C.C. Hsing, "Parameter adaptive control of 
multivariable systems," International Journal of Control, 22, 
No. 3, 1975, pp. 313-327. 


284 


41. Sobel, K., H. Kaufman, and L. Mabius, "Model reference output 
adaptive control systems without parameter identification," 18th 
IEEE Conference on Decision and Control, Ft. Lauderdale, FL, 
Dec. 1979, pp. 347-351. 

42. Sobel, K. , H. Kaufman, and L. Mabius, "Implicit adaptive control 
systems for a class of multi-input multi-output systems,” IEEE 
Trans, on Aerospace and Electronics Systems, Sept. 1982, pp. 
576-590. 

43. Broussard, J.R. and M.J. O'Brien, "Feedforward control to track 
the output of a forced model," Proc. 17th IEEE Conference on 
Decision and Control, Jan. 1979, pp. 1149-1154. 

44. Balas, M.J. and C.R. Johnson, “Toward adaptive control of large 
structures in space," Applications of Adaptive Control, Academic 
Press, 1980, pp. 313-344. 

45. Johnson, C.R., "Adaptive modal control of large flexible 
spacecraft," AIAA Journal of Guidance and Control, Vol. 3, No. 
4, July-August 1980, pp. 369-375. . 

46. Rohrs, C.E., L. Valavanl, M. Athens, and G. Stein, "Robustness 
of adaptive control algorithms in the presence of unmodeled 
dynamics," Proc. 21st IEEE Conference on Decision and Control, 
Dec. 1982, pp. 3-11. 

47. Bar-Kana, I., H. Kaufman, and M. Balas, "Model reference 
adaptive control of large structural systems,” AIAA Journal of 
Guidance and Control, Vol. 6, No. 2, March- April 1983, 

pp. 112-118. 


285 


48. Workshop Briefing Charts, Study Task No. 9 - Berthing, NASA 

Space Station Task Force Concept Development Group, Dec. 5-9, 
1983. 

49. Craig, M. , Study Task No. 11 - Configuration and Controllability 

Studies, NASA Space Station Task Force Concept Development 

Group, Dec. 16, 1983. 

50. Clough, R.W. and J. Penzien, Dynamics of Structures, McGraw- 

Hill, N.Y., N.Y., 1975. 

51. Juang, J.N. and S.J. Wang, "An Investigation of quasi-inertial 
attitude control for a solar power satellite,” Space Solar Power 
Review, Vol. 3, 1982, pp. 337-352. 

52. Hahn, W., Theory and Application of Lyapunov's Direct Method, 
Prentice-Hall, Inc., New Jersey, 1963. 

53. Lyapunov, M.A., "Investigation of a singular case of the problem 
of stability of motion," Mat. Sbornik 17, 1893, pp. 252-333 

54. Bar-Kana, I., "Direct multivariable model reference adaptive 
control with application to large structural systems," Ph.D. 
Dissertation, Rensselaer Polytechnic Institute, Troy, New York, 
May, 1983. 

55 . Landau, I.D., "Design of discrete model reference adaptive 
systems using the positivity concepts," in Proc. 3rd Symp. on 
Sensitivity, Adaptivity and Optimality, Ischia, Italy, 1973, 

pp. 307-314. 


286 


NASA— JPL-ComL, LA, Caltf. 


TECHNICAL REPORT STANDARD TITLE PAGE 


!* Report No> JPL 85-57 


2. Government Accession No. 3. Recipient's Catalog No. 


4. Title and Subtitle 

Dynamic Modeling and Adaptive Control for 
Space Stations 


5. Report Date 

July 15, 1985 


6. Performing Organization Code 


7. Author (s) 


Che-Hang Charles Ih and Shyh John Wang 


9. Performing Organization Name and Address 

JET PROPULSION LABORATORY 
California Institute of Technology 
4800 Oak Grove Drive 
Pasadena, California 91109 


12. Sponsoring Agency Name and Address 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Washington, D.C. 20546 


15. Supplementary Notes 


8. Performing Organization Report No. 


10. Work Unit No. 


11. Contract or Grant No. 
NAS 7 -9 18 


13. Type of Report and Period Covered 

JPL Publication 


14. Sponsoring Agency Code 


16. Abstract of all large space structural systems, space stations present a unique 
challenge and requirement to advanced control technology. Their operations require 
control system stability over an extremely broad range of parameter changes and high 
level of disturbances. During shuttle docking the system mass may suddenly increase by 
more than 100% and during station assembly the mass may vary even more drastically. 

These coupled with the inherent dynamic model uncertainties associated with large space 
structural systems require highly sophisticated control systems that can grow as the 
stations evolve and cope with the uncertainties and time-varying elements to maintain 
the stability and pointing of the space stations. 

This report first deals with the aspects of space station operational properties 
including configurations, dynamic models, shuuttle docking contact dynamics, solar panel 
interaction and load reduction to yield a set of system models and conditions. A model 
reference adaptive control algorithm along with the inner-loop plant augmentation design 
for controlling the space stations under severe operational conditions of shuttle 
docking, excessive model parameter errors, and model truncation are then investigated. 
The instability problem caused by the zero-frequency rigid body modes and a proposed 
solution using plant augmentation are addressed. Two sets of sufficient conditions 
which guarantee the globablly asymptotic stability for the space station systems are 
obtained. 

The performance of this adaptive control system on space stations is analyzed through 
extensive simulations. Asymptotic stability, high rate of convergence, and robustness of 
the system are observed under the above-mentioned severe conditions and constraints 
induced by control hardware saturation. 


17. Key Words (Selected by Author(s)) 

Astronautics (General) 

Engineering (General) 

Computer Programming and Software 
Systems Analysis 


18. Distribution Statement 


Unlimited/Unclassified 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21 . No. of Pages 

22. Price 

Unclassi fied 

Unclassi fied 




** 


























